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Abstract 

The scaling theory of irreversible aggregation is discussed in some detail. First, 
we review the general theory in the simplest case of binary reactions. We then 
extend consideration to ternary reactions, multispecies aggregation, inhomogeneous 
situations with arbitrary size dependent diffusion constants as well as arbitrary 
production terms. A precise formulation of the scaling hypothesis is given as well 
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as a general theory of crossover phenomena. The consequences of this definition are 
described at length. The specific issues arising in the case in which an infinite cluster 
forms at finite times (the so-called gelling case) are discussed, in order to address 
discrepancies between theory and recent numerical work. Finally, a large number of 
exactly solved models are reviewed extensively with a view to pointing out precisely 
in which sense the scaling hypothesis holds in these various models. It is shown that 
the specific definition given here will give good results for almost all cases. On the 
other hand, we show that it is usually possible to find counterexamples to stronger 
formulations of the scaling hypothesis. 

Key words: aggregation, scaling, exactly solved models 
PACS: 05.70.Ln, 82.70, 82.35.+t 



1 Introduction 

1.1 What is this all about 

What do a glass of fresh milk, a planetary system in formation and Los Ange- 
les on a smoggy day all have in common? At first sight, not much, even from 
the slightly skewed viewpoint of a physicist. The various phenomena involved 
in these systems have quite different physics and the interesting questions are 
wholly unrelated. Nevertheless, there is one aspect which these three (and a 
large number of other) systems share: in all of them one finds some process of 
irreversible aggregation, the characterization of which is of some importance. 
To be specific, in the case of milk, we are dealing with small globules driven 
by buoyancy which coalesce irreversibly to form cream; in a planetary system, 
on the other hand, one has planetesimals colliding inelastically to form ever 
larger blocks, eventually leading to planet formation. Finally, in aeorosols such 
as those found in smog, we have airborne particles so small that Brownian mo- 
tion is (frequently) an important transport mechanism. There the aggregation 
mechanism is due to the presence of van der Waals forces between the particles, 
which are considerably stronger than the effects of thermal agitation. 

Let us therefore review what these systems have in common: in all of them, 
we have identifiable aggregates of particles playing an important role. These 
aggregates grow by sticking to each other. Specifically, they are usually not 
limited to growing by accretion of the very smallest clusters, which is another 
topic altogether. One has, therefore, a situation in which there is originally 
quite a large number of small particles, and these then coalesce as time goes 
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on. No allowance is made for the possibility of break-up. One is thus not 
concerned with a problem of approach to equilibrium, but with an entirely 
dynamical problem 1 . 

Aggregation phenomena are, as hinted at above, quite common in nature. 
Apart from the examples mentioned above, we may mention astrophysics (for 
further discussion of the subject with references see e.g. [55]), cloud physics 
(see [86]) and polymer chemistry, where a great deal of the theoretical work 
underpinning the gelation transition actually originated (see [95]). Even to 
attempt to give a fair idea of the variety of applications of aggregation and 
related concepts, however, is beyond both my abilities and the scope of this 
work. The reader specifically interested in the applications to aerosol physics 
is referred to the review by Drake [23] which addresses a large number of 
such issues in detail. Friedlander [32] and Hidy [37] are also important general 
references on the subject. 

What, then, are the issues we want to address when dealing with such a sys- 
tem? Overall, the questions most frequently asked fall into two categories: 
the first concerns the sizes and the second the morphology of the aggregates. 
Let us deal shortly with the latter first. It is, of course, quite obvious that if 
two liquid droplets coalesce, they will, under the influence of surface tension, 
rapidly relax to a spherical shape. Slightly more complex liquid systems may 
behave in a less obvious way (think for example of a water drop on an inclined 
glass plane, where surface tension and gravity act together), but the geome- 
try still remains fairly straightforward. Matters change dramatically, however, 
when solid particles stick together. In this case, depending on the degree of 
rigidity of the bonds, the aggregates may or may not rearrange. The simplest 
case, because it is quite well-defined, is the one in which there is no rearrange- 
ment whatever. Then, as is well-known, the aggregates grow fractal shapes 
and their morphology can only be characterized at the statistical level. The 
main tool used in order to characterize different morphologies is the so-called 
mass-radius relationship 

R(m) = const. -m 1/D f, (1.1) 

for large values of the mass m. Here R(m) is some average measure of the 
radius of the aggregate as a function of its mass m and Df is an exponent 
known as the fractal dimension. A considerable amount of work has gone into 
the study of such systems (for a review and for an introduction to the extensive 
literature, see [47,78]). The upshot is that a fractal dimension can indeed be 
defined, but that it depends, although not very sensitively, on some of the 

1 It is argued by some that this is unrealistic, since on fundamental grounds back- 
ward reactions always exist. To this the short answer is that I shall only be interested 
in time scales for which the backward processes can indeed be neglected. 



4 



details of the transport mechanism. In particular, ballistic transport yields 
more compact structures (and hence larger values of the fractal dimension 
Df) than diffusion, and low sticking probability has a similar effect. 

A short summary of the above discussion could therefore be stated as follows: 
morphology is extremely important, but it is strongly dependent on the spe- 
cific mechanisms involved in each particular system (many features of general 
interest such as, say, the appearance of complex self-similar structures in the 
case where no restructuration is present, occur quite generally, but detailed 
features, such as the actual values of the fractal dimension, do depend on 
specifics). 

What will principally concern us in this review is the study of the masses of the 
aggregates. Specifically, since there are always a large variety of masses at one 
given time, we need the so-called cluster size distribution as a function of time. 
(Note that in the following, by size, we shall always mean mass and never some 
measure of geometric size such as radius or volume.) By this we mean that we 
study the function c(m, t) giving the concentration of aggregates of mass m at 
time t. Roughly speaking, one uses the following strategy: in many physical 
situations, it is possible to express the rate of aggregation as the encounter 
probability between aggregates of masses m and m', multiplied by a reaction 
rate K(m, m!) depending in quite a general way on the two masses m and m' . 
The possibility of describing the aggregation process in this very simplified 
manner depends on several assumptions, which will be shortly discussed later. 
The most important is a mean-field hypothesis, which states that no spatial 
correlations develop during the course of the reaction. From the theoretical 
point of view, this is an extremely severe assumption, but in practice it is 
often quite well satisfied. 

If this can be done, it means that we have summarized all the information con- 
cerning the complicated and specific physics of the various individual systems 
in the matrix K(m,m') characterizing the mass-dependence of the reaction 
rates. It should be emphasized that the determination of the rates K(m, m!) is 
in general by no means straightforward. However, as it depends intimately on 
the specific physics of the system involved, we shall not consider it here, but 
rather always assume the rates K(m, m') to be externally given. One is then 
led to a system of infinitely many coupled non-linear rate equations. These 
still present a formidable problem, which usually cannot be solved by analytic 
means. As for numerical studies, while they are certainly very valuable, they 
are necessarily limited in their ability to deal with either large cluster sizes 
or large times. In this review, we extensively describe an approach due origi- 
nally to Friedlander, and then considerably expanded upon by Ernst and van 
Dongen, namely the scaling approach. The crucial idea consists precisely in 
limiting oneself to the case of large masses and large times. One may then, 
as was shown, take only the crudest features of K(m, m') into account and 



5 



arrive at quite precise predictions for the overall features of the cluster size 
distribution. 

Let us illustrate this by a simple and classical example: consider particles 
moving with a mass-dependent diffusion constant D{m) having a radius R{m) 
growing as in (1.1). It is then reasonable to assume that the reaction rates are 
given by 

K(m, m') = {D{m) + D(m'))(R(m) + R(m')) d -\ (1.2) 



where d is the dimension of space. This follows from the elementary fact that 
the rate at which a pointlike particle diffusing with diffusion constant D is 
captured by a spherical trap of radius R is of the order DR d ~ 2 . In order to 
proceed, we need to know the asymptotic behaviour of the diffusion constants 
D{m). Let us assume, in keeping with tradition, a fully general dependence of 
the type 

D{m) = const. ■ m~ 7 . (1.3) 



The scaling theory referred to above then states the following facts concerning 
the cluster size distribution c(m, t), that is, the number of aaggregates of mass 
m per unit volume present in the system at time t: 

(1) The so-called typical cluster size s(t) grows as t z , where the exponent z 
is given by 

D ' (1.4) 



Df + T-id-2) 



This means, as we shall see in far greater detail, that if we consider 
an appropriate rescaling of the cluster size distribution on this scale, it 
collapses onto a single time-independent distribution. 

(2) From the above follows also that 

oo 

M p (t) := J m p c(m,t)dm = const. ■ &~ 1)z , (1.5) 
o 

at least for sufficiently large values of p. Further work shows that if 7 > 0, 
which is certainly the most common case, (1.5) actually holds for all 
values of p. When 7 < 0, on the other hand, the critical value for p is not 
known, but is known to lie between and (d — 2)/Df — 7. 

(3) Finally, the behaviour of the aggregates that are much smaller than s(t) 
can be evaluated. Thus it can be shown that such aggregates decay as 
exp(— const. ■ t( d-2 )/ D /) if 7 < 0, whereas they decay as a power law 
otherwise. 
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The way in which these and similar conclusions are drawn will be the topic of 
this paper. 

At this stage, the better read may be experiencing some doubts: all of the 
above example, in fact, was essentially known in the sixties, see for example 
[96]. So why am I talking about such things, and what is new? The answer 
comes in two parts. 

First, I have attempted to show how scaling theories can be constructed in a 
nearly automatic manner for almost any type of problem involving irreversible 
aggregation. By this I mean that I have treated a large number of extensions 
(source terms, aggregates with variable compositions, inhomogeneous situa- 
tions and, finally, the effect of three-body interactions). In all of these cases I 
have shown a perfectly routine way of obtaining a scaling theory, which yields 
sensible results. The point of this excercise, which may be said to be the main 
purpose of this review altogether, is to show how to compute scaling functions 
systematically in almost any reasonable situation by simply turning a crank. 
In this respect, the case of three-body interactions is particularly instructive: 
indeed, I was just busy working out just how the crank worked in this particu- 
lar case, when suddenly it failed in quite an unexpected way. What happened 
is that an exact result, quite easy to derive, conflicted with the equally evi- 
dent scaling theory I had devised. Appendix B presents my present view of 
what went wrong, and should be borne in mind whenever constructing any 
kind of scaling theory along these lines. Another rather tricky issue concerns 
gelation: in this case, all the original literature had predicted certain values 
for various exponents, which were clearly refuted in the remarkable numerical 
simulations of Lee [56]. An attempt is made to show that the previous argu- 
ments [59,116,117] were indeed inconclusive 2 , to show in detail why and to 
present some indication of how a theory of Lee's work might go. 

As the title of this work indicates, however, there is yet another side to this 
review, which concerns exactly solved models. Indeed, there exists a consid- 
erable body of knowledge concerning those few models for which the rate 
equations can be solved analytically. This is extremely useful in order to test 
the scaling approach. The problem is that the exact nature of the error made 
when using the scaling ansatz is not known a priori On the other hand, when 
one is dealing with exactly solved models, one has a full overview of the so- 
lution's behaviour, so that a detailed comparison with the scaling predictions 
becomes possible. To this end, however, it is necessary to investigate the ex- 
act solutions in the case of arbitrary initial conditions, which has not, to my 
knowledge, been done earlier to any great extent. Here I present such exact 
solutions and show that scaling, in a very specific sense of the word, indeed 
invariably holds, whereas many of the consequences which are ordinarily be- 



2 or, in plain English, wrong 
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lieved to follow directly from it, actually fail. This may be viewed as a list 
of cautionary counterexamples, which show the need for care in stating what 
follows from scaling and what does not (thus, I have been characteristically 
careless in asserting, for example, the behaviour of the small clusters stated 
in the example above). The whole issue is rather complex, but a simple sum- 
mary 3 may be given as follows: whenever a given property of the cluster size 
distribution involves a finite fraction of the total mass, one may be reasonably 
confident to obtain it correctly via scaling. Whenever it does not, however, 
counterexamples are possible, and indeed, in generic cases, rather likely. 



1.2 Outline of the Paper 

Let me here briefly summarize the layout of this paper. In Section 2 I state the 
known results on the properties of existence, uniqueness and mass conservation 
for the solutions of the rate equations. These are rather few, but quite instruc- 
tive. Indeed, they are probably also in some sense optimal, according to the 
results suggested by the scaling approach. In Section 3 we develop the scaling 
approach in several stages. First, we show how scale invariant solutions exist 
whenever the rate kernel K(m, m') is homogeneous in its arguments. This in 
itself, however, says nothing about the general solutions of the equations. To 
this end, we require the scaling hypothesis, which states that all solutions which 
start from an specific class of initial conditions eventually approach this scale 
invariant solution at large times. We then go on to discuss the precise sense in 
which such convergence should be assumed. This question may seem to be en- 
tirely technical and devoid of interest, but in fact it lies at the very heart of the 
issues concerning the validity of scaling. I shall suggest that the appropriate 
choice of mode of convergence is weak convergence, where suitably normal- 
ized measures are defined in terms of the cluster size distribution c(m,t). To 
express this in non-technical terms, I will say that a cluster size distribution 
c(m, t) converges weakly to a given distribution <&(#) with respect to a typical 
size s(t) if the expectation value of any function which varies smoothly on 
the scale s(t) tends to the corresponding expectation value as t — > oo. If this 
specific type of convergence is used, it can be shown, as will be done exten- 
sively in Section 4, that scaling obtains under all reasonable circumstances. 
In other words, all the issues that have been raised concerning so-called "vi- 
olations of scaling" are seen to arise from the fact that the correct questions 
are not being asked, or equivalently, that conclusions are being drawn from 
the scaling hypothesis, which are not legitimate once scaling is defined in this 
precise manner 4 . I then proceed to derive a general integral equation for the 



3 valid at least for the non-gelling case 

4 Before scaling is defined in some way or other, of course, the whole debate remains 
void 
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scaling function <3>(x). This is simply a mild generalization of the Ernst-van 
Dongen equation and allows in the usual way to derive the standard cases I, 
II and III 5 for non-gelling aggregation. A similar equation is also obtained 
for the gelling case, and some very partial results for the behaviour of the 
solution at the origin are obtained. These, while still quite incomplete, clearly 
show that the original claims of [58] concerning the so-called r-exponent in 
the gelling case were incorrect. The resulting claims are now more nearly in 
agreement with the numerical findings of Lee [56] . The theory further suggests 
a new scaling relation which should, perhaps, be verified numerically. After 
these general remarks, I proceed to show how the scaling approach can always 
be obtained in an essentially mechanical fashion. In particular, extensions to 
monomer production, spatially inhomogeneous systems, aggregation with ag- 
gregates of many species multicomponent aggregation and, finally, three-body 
aggregation, are all worked out in some detail. A theory of crossover, for the 
case in which one aggregation mechanism is supplanted by another at large 
sizes, is also presented. In Section 4, on the other hand, we run through most 
of what was worked out in Section 3, but now using exactly solved models. 
Our main workhorse, of course, is the constant kernel, for which a great deal 
is known, and which can be solved exactly in almost any of the extensions one 
cares to think about. However, we also consider the sum and product kernels, 
and present some results on the general bilinear kernel, in order to illustrate 
the concepts of crossover previously developed. We also discuss two more re- 
cent exactly solved models, both of them variations on the constant kernel: 
one is a set of reaction rates K(k, I) which depend on the parity of the sizes 
k and I of the reactants, but on nothing else; the other is also only defined 
for discrete values of the masses and is given by 2 — q k — q l , where q is an 
arbitrary number between zero and one. Both of these can be solved exactly 
to a large extent, and both display strange effects which have occasionally 
been argued to be "violations of scaling" . Again, a discussion of these models' 
behaviour shows the contrary, but also sheds considerable light on the limits 
of our definition of scaling, which has a very hard time making any reasonable 
statements about clusters of fixed size. But then, as these models in fact show, 
the behaviour of fixed size clusters can be completely unexpected. Finally, in 
section 5, I discuss the issue of models for which the mean-field theory does 
not apply. In the cases I discuss, spatial correlations are built up by the inter- 
play between transport and reaction. Whenever the reaction is slow, these are 
negligible, but as soon as transport becomes the rate-limiting factor, spatial 
correlations invalidate the use of the rate equations. In this case, no straight- 
forward application of the formalism developed in the previous sections is, 
of course, possible. Scaling concepts, however, remain extremely useful as we 
shall see, as they allow to analyze the information one has in a systematic 
way. 



5 These will be described in due time. 
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1.3 Various Topics of Interest not Treated Here 

There are, as always, a considerable number of interesting subjects which 
could, and perhaps ought to, have been included in this review, but which 
have been left out, for reasons either of space, of limited competence on my 
part, or else finally because, although undeniably beautiful, certain subjects 
do not fall readily within the scope of this paper as it has been described in the 
previous subsections. A necessarily incomplete list might include the following 

(1) The rate equations are not, of course, primary: rather they are derived, 
after some approximations, from a master equation for some appropri- 
ate Markov process. Using then the so-called f2-expansion devised by van 
Kampen, one may derive the rate equations. This leads to the discussion 
of corrections due to particle number fluctuations, for a thorough discus- 
sion of which the reader is referred to [70,104,105,109,11]. A considerable 
literature has also arisen in the literature on probability concerning such 
models and their connection to the Smolchowski equations. For a re- 
view, see [2]. An altogether remarkable claim of convergence of the above 
stochastic process to the rate equations treated in this paper is made in 

[si]. ; 

(2) It is possible to go into much greater detail concerning the exact solution 
of the full bilinear kernel. I have essentially limited myself to analyzing 
various types of crossover and showing that scaling behaviour arises for 
arbitrary initial conditions. Much more explicit forms can in fact be ob- 
tained, and a very interesting combinatorial interpretation of them can be 
given. The reader is referred to [90,91,92] for further discussion of these 
issues. 

(3) The list of extensions to the aggregation equations is of necessity incom- 
plete. Thus one findstreatments of aggregation with replication or with 
annihilation [52,6] as well as much more. The elements presented in this 
paper could also be multiplied in various ways. This has not been at- 
tempted. Rather, I have tried to emphasize the applicability of a general 
method, which could then be carried over to arbitrary cases as desired. 
The systems studied were chosen on the basis of a wholly subjective feel- 
ing concerning both their practical importance and the relevance of the 
conceptual issues attached to them. 

(4) As will be seen, there exist several cases in which the polydispersity expo- 
nent r cannot be determined via a straightforward study of the integral 
equation for the scaling function. For these a very powerful method was 
developed by Cueille and Sire [18] in order to obtain rigorous bounds 
on the r exponent. The method turns out to be so powerful that these 
bounds are often tantamount to exact evaluations: the upper and the 
lower bounds coincide to, say, four or five decimals. They also give a sys- 
tematic approach to make ever better bounds. This should therefore be 
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viewed as a solution of the problem of evaluating r for those particular 
cases. However, the approach is quite subtle, so that I must refer the 
reader to the original paper. 
(5) My treatment of non mean-field models is completely insufficient. This, 
however, may perhaps be excused in the light of the fact that my aim here 
is to present the scaling theory of aggregation, not all the beautiful results 
in some way connected with irreversible processes and models for chemical 
reactions, such a task being quite beyond my abilities. I have therefore 
strictly limited myself to such models as were strictly aggregation models, 
and have only shown how the scaling theory can be brought to bear in 
such case. Even so, I have not, of course, been able to do the subject 
justice. The references given in this Section may get the reader started. 



2 Irreversible Aggregation: The Mean-field Approach 

2. 1 Generalities 

As already stated in the Introduction, the phenomenon of irreversible aggre- 
gation is described as follows: Aggregates of mass m and m! react to form 
aggregates of mass m + m! . In mean-field theory, one writes down rate equa- 
tions for this process, asssuming that 

A m + A m r - > A m + m i, (2-1) 
K(m,m ) 

where the K(m, m!) are the rates at which the aggregation process takes place. 
The index m may range over a discrete or a continuous range of values. If 
we then denote by c(m, t) the concentration of aggregates A m at time t, we 
obtain the following rate equations to describe the dynamics of the cluster size 
distribution 

d t c(m, t) = \J dm, dm 2 K(m u m 2 ) c(m u t)c(m 2 , t) x 

x [5{mi + m 2 — m) — 5{mi — m) — <5(m 2 — m)] . (2.2) 

Here, as in the rest of this paper, I choose the convention that m runs over 
continuous values, since the modifications to obtain the discrete case can al- 
ways be viewed as specializations to singular measures. However, since in many 
cases it is much more convenient to solve the case in which only discrete multi- 
ples of a certain mass occur, let me introduce some notations which will recur 
throughout this paper: if we have discrete initial conditions, that is, if there 
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is an mo such that 



oo 



c(m, 0) = Ck{0)5(m - km Q ), 
k=i 



(2.3) 



then the solution c(m, t) at all times has the form 



c(m, t) — Ck(t)5(m — km ). 
k=i 



(2.4) 



Here the r.h.s of (2.4) defines the functions Ck(t). A very important special 
case is that of so-called monodisperse initial conditions, for which Cfc(0) = Sk,i- 

The derivation of (2.2) involves certain hypotheses, which I shortly discuss: 
First, the effect of cluster morphology on the rates of aggregation must ei- 
ther be neglected or be taken into account using some appropriate averaging 
techniques. Second, we need to assume that no spatial correlations between 
the clusters build up as the reaction proceeds. It is well-known that such an 
assumption is by no means trivial. In fact, one may argue that most realistic 
systems will, for sufficiently large times, display such correlations. However, it 
can be shown that these can indeed be neglected when the transport mecha- 
nisms responsible for bringing the clusters to react are significantly faster than 
the reaction step. Since this is a fairly common situation in practice, we will 
restrict ourselves to this mean-field situation except in Section 5. 

Equations (2.2) can, of course, be generalized in many ways. One may include, 
among others, the following additional effects: 

(1) Reactions of higher order than binary. This is certainly important when 
the volume fraction of the aggregates grows in time, which occurs when- 
ever a model displays fractal growth. In this situation, however, it is quite 
unlikely that the hypotheses necessary for mean-field theory to hold re- 
main valid when these terms become dominant. However, their effect as 
corrections can be estimated within this framework, which allows to de- 
cide whether such corrections are likely to be important or not. 

(2) The parameter m can be generalized to be a multicomponent vector. 
This describes, for example, situations in which the aggregates consist of 
various substances. In this case the scaling theory must be generalized 
somewhat to incorporate various different mass scales which describe dif- 
ferent aspects of the cluster size distribution, namely the average as well 
as the spread in composition, which both scale differently, due to the 
central limit theorem. 

(3) Spatially inhomogeneous situations can be considered. In this case, diffu- 
sion terms must be explicitly added to the equation. The diffusion con- 
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stants involved have, of course, in general a non-trivial mass-dependence, 
which leads to quite complicated equations. 
(4) Finally, monomer production can be introduced. One may either consider 
the case in which this occurs homogeneously throughout the system, or 
on the contrary, we may assume that monomer is being produced at a 
given point and diffuses as it aggregates. 

As we shall see, the scaling approach can be extended to give a straightforward 
description of these various systems. However, before we proceed, I would like 
to review shortly some of the rigorous results known on the existence and 
uniqueness properties of solutions of (2.2). 

2.2 Existence and Uniqueness Results 

In the following, we shall treat the discrete case of (2.2) first: results have usu- 
ally first been obtained for it, since the general case involves some additional 
difficulties. These equations represent an extremely challenging problem from 
the point of view of pure mathematics. The reason is that they consist of an 
infinite number of coupled nonlinear ordinary differential equations. There do 
not exist standard results for such cases. 

In this case, existing results fall into two classes: First, those which make an 
assumption of the type: 



Here, and throughout the rest of this paper, reaction rates will always tacitly 
be assumed to be positive. From such an assumption it can be shown [112] 
that the c(m, t) decay exponentially in m for all t for which the solution exists. 
From this one shows that the quantity 



remains constant as long as the solution is defined. This in turn allows to 
prove existence and uniqueness of the solution for all times. Note that it is 
formally trivial to prove the constancy of Mi(t), which physically corresponds 
to the conservation of the total mass contained within the aggregating system. 
However, from a rigorous point of view, it is necessary that 




(2.5) 




(2.6) 



m=l 



■oo 



/ 



mK(m,m')c(m,t)c(m' ,t)dmdm < oo 



(2.7) 



o 
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for all times. 

These results are in a sense optimal: they yield both existence and uniqueness, 
as well as mass conservation. Further, it is expected that mass-conservation 
will fail whenever the hypothesis (2.5) does. However, this hypothesis is of- 
ten too restrictive. Thus, a reasonable model for aggregation of branched ag- 
gregates (without loops) assumes that the number of reactive sites on each 
aggregate grows as m, from which follows 

K(m,m') =Cmm', (2.8) 

which obviously does not satisfy (2.5). In this case, it was first shown [76,77], 
that no solution satisfying constancy of Mi(t) could exist beyond a certain 
time. This negative result was later put into proper perspective by the discov- 
ery [115,58] of an exact solution for all times, for which Mi(t) decays after a 
finite time. This phenomenon, which is found to be linked to the appearance 
of an infinite cluster at finite time, is known as gelation. It is of considerable 
importance and we shall discuss it extensively further on in this paper. 

This led to a different set of existence results [58] , which did not aim to prove 
either constancy of Mi(t) or any decay condition leading to (2.7). It could 
then be shown in [58] that, if 

K(m,m') < r{m)r{m') r{m) = o(m), (2.9) 

then a solution of (2.2) with finite and monotonically decreasing Mi(t) exists 
for all times, but no statement was made either concerning the conservation 
of mass or uniqueness of the solution. The latter is certainly unfortunate, 
but is intimately related to the technique of proof used. This issue has now 
essentially been solved in [81], where many other striking results are shown. 
The result claimed is that under the hypothesis (2.9) uniqueness holds. The 
former is, on the other hand, only to be expected. In fact, as we shall see, the 
scaling theory strongly suggests that the bound (2.5) is, in fact, exactly the 
one that separates gelation at finite time from regular behaviour at all times. 
As a further indication of this fact, the model 

K(m,m') =m a 5 m , m , (2.10) 

can be rigorously shown to violate mass conservation if a > 1 [60], whereas 
it satisfies (2.5) — and hence mass conservation — if a < 1. The marginal case 
a — 1 is discussed in greater detail in [10] . 

Let me shortly describe the technique involved, since it turns out to be of 
fairly general applicability: one first defines a finite, mass-conserving system, 
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for which standard results can be invoked to guarantee existence of the solution 
for all times. A compactness argument is then used to show that the solutions 
to the finite systems have a convergent subsequence. It can then be shown that 
the limit of this subsequence satisfies a weak form of (2.2), and finally that it 
satisfies (2.2) in a strong sense as well. The mode of convergence, however, is 
not strong enough to conclude that the mass contained in the limiting solution 
is the same as that in the finite approximants, which is what gives to this 
approach the requisite generality. The extension to the case of continuously 
variable m presents some technical difficulties. These were adressed by Bak 
[3] in the case of a reaction kernel bounded above and below by a constant. 
The general case was treated by Ball and Carr in [5]. 

Another extremely important issue from the point of view of applications 
concerns the extension to non-uniform systems This involves replacing (2.2) 
by a reaction diffusion equation. This has been treated by Slemrod [88] , along 
lines similar to those sketched above. Here, however, it should be emphasized 
that even the "standard" results on the finite approximants are highly non- 
trivial in this case. While the results obtained are remarkable indeed, I believe 
that they are probably not yet optimal: indeed, the strongest results obtained 
concern diffusion constants bounded from below and above by a constant. In 
a realistic situation, however, these always decrease with mass; an extension 
of these results to more general cases would certainly be very desirable. 

Finally, it should be emphasized that the issues adressed here are not al- 
together idle: existence may indeed fail under some circumstances. There is 
considerable evidence, for example, that if K(mo,m)/m — > oo as m — > oo for 
fixed m , then either no solutions exist or these solutions have the peculiar 
property of violating mass conservation from the very start. In this case it 
is an open problem whether the solution exists at all for a sufficiently wide 
set of initial conditions. For more details on these fascinating questions, see 
[107,16,55,56]. 



3 The Scaling Hypothesis 

3. 1 Basic Concepts: Moments and Typical Size 

The crucial observation underlying the scaling approach is the following: at 
large times the value of most quantities of interest — namely those that result 
from an average over the whole cluster size distribution — can be computed 
using a single "typical size" which grows indefinitely with time. Under these 
circumstances, we expect that only the coarsest features of the K(m, m!) will 
be relevant. Let us therefore assume that K(m,m') is asymptotically homoge- 
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neous of degree A, that is, that there exist a A and a function K(m,m f ) such 
that 



K(am, am') = lim s x K(sm, sm!) 



(3-1) 



The function K(m,m') is then clearly homogeneous of order A. To describe 
K(m, m') in greater detail, note that it can be written as 



m 



K(m, m') = m k — , 
\ m I 



(3.2) 



where k(z) is a function subject to the symmetry condition 
k{z) = z x k{l/z), 



(3.3) 



but which is otherwise quite arbitrary. It can therefore be chosen arbitrarily, 
say, between and 1. We now define the exponent \i and the prefactor T which 
describe the (power law) behaviour of k(z) near the origin. More specifically: 



k(z) = rv[i + o(i)] 



(3.4) 



One additionally defines v as the exponent describing the behaviour of k(z) 
at infinity. It follows from (3.3) that 

v — A — /x. (3.5) 



It follows from these definitions that there exist reaction rates corresponding 
to arbitrary values 6 of fj, and u, as evidenced by 



K(m, m) 



mrm 



mm 



(m < m') 
(m > m') 



(3.6) 



At this stage, we should rephrase the existence and uniqueness results de- 
scribed in the previous subsection. Condition (2.5) is equivalent to A < 1 
and i/ < 1. It is indeed presumably optimal for the existence of a global 
mass-conserving solution of (2.2). On the other hand, condition (2.9) which 
guarantees existence only, is equivalent to the conditions A < 2 and v < 1. 

6 This issue has caused some confusion in the literature, as a lot of work has gone 
into the study of kernels of the form m a m' X ~ a + m x ~ a m' a , which do not generate 
the whole range of possible jj, and v values, as opposed to (3.6) 
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However, a straightforward extension of the proof given in [58] shows that 
the former assumption can be dispensed with. The only essential condition is, 
therefore, v < 1. 7 

Let us temporarily assume that K{m,m') is exactly homogeneous. From this 
assumption follows [36] that, if c(m,t) is an arbitrary solution of (2.2), then 
so is 

{T a>b c) (m, t) = a x+l b- l c{am, bt) (3.7) 



for arbitrary positive a and b. The existence of such a symmetry group moti- 
vates the search for solutions which remain invariant under it. In other words, 
we ask whether solutions exist, which, as time evolves, merely transform into 
similar solutions generated from a constant cluster size distribution by a vari- 
able group transformation of the form (3.7). In the so-called regular case, that 
is, when mass conservation holds, see (2.5) , this can only occur if the group 
transformation belongs to the subgroup of (3.7) which maintains the total 
mass in the system invariant, that is T aa x-i. From this follows that such a 
solution necessarily has the form 

c(m, t) = Wt- 2,{l - x) § (m/t 1/(1 - A) ) , (3.8) 



where W is a constant factor which will be of use later 8 . The gelling case, 
as we shall see below, requires more careful considerations. Quite generally, 
however, we shall say that a solution is of scaling form if it can be written in 
the form 



c(m,t) = Ws(t)- 2 ^^j, (3.9) 



where s(t) is an arbitrary (growing) function of time, which is called the typical 
size. From this follows, as we show in greater detail later, that if we define the 
moments of the distribution as follows 

oo 

M p (t) = J m p c{m,t)dm, (3.10) 
o 



The exponents A and \i above can be defined in a rather rough manner, except if 
we consider the limiting cases A = 1 and v = 1. For these it is essential that there 
be no logarithmic corrections to the power law behaviour, see e.g. [10] for details. 
8 Specialists will recognize here the separation constant w of the papers of Ernst 
and van Dongen 
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one finds as an estimate for the order of magnitude of s(t) 
M v+ i(t) 

(3 ' n) 



at least when p is larger than a given p , the value of which depends on the 
details of the system under consideration. In fact, in many contexts, the typical 
size is defined as some such moment ratio: such are, for example, the weight 
average, defined as M 2 (t)/M 1 (t) and the ^-average, given by M 3 (t)/M 2 (t). 
The constant W will allow us to maintain any such definitions, while keeping 
simple normalizations for the function 

We shall show in the following sections that, if such a solution exists, the func- 
tion <&(#) must satisfy a certain integral equation. Conversely, it is straight- 
forward to see that any solution to the integral equation also defines via (3.8) 
a self-similar solution of (2.2). 

The scaling hypothesis now reduces to the following statement: every solution 
starting from a "narrow" cluster size distribution (that is, one that does not 
have already some kind of power law distribution at large masses) evolves for 
large times into a self-similar solution of the type we have just described 9 . We 
are therefore assuming that these specific solutions, which can be described 
quite accurately in a broad variety of cases, are in fact the only relevant 
solutions at large times, as they are approached from essentially any initial 
conditions. It should be emphasized that there is no proof of this so-called 
"scaling hypothesis", nor is there any idea of how such a proof might go. 
However, the numerical evidence for it is good, and we shall see that it is 
satisfied in all exactly solved models. 



3. 2 Different Kinds of Scaling 



In the following, we define more accurately in which sense we expect solu- 
tions of (2.2) to tend towards a scaling solution. This issue is of fundamental 
importance, as it will accurately determine the kind of conclusions one may 
legitimately draw from scaling and those one may not. 

To justify the definition I shall make, let me first state the nature of the 
object we are interested in: this is the cluster size distribution c(m,t). The 
information it gives us is the following: how many clusters are there at time t 
in the mass range going from m to m+dm. As such, the cluster size distribution 

9 In fact, we shall see that, quite often, even initial conditions with power-law tails 
will approach the same scaling limit as narrow ones. The above is merely the most 
conservative formulation of the hypothesis. 



18 



is essentially a time- dependent measure on the set of all masses. We therefore 
need a notion of convergence adapted to measures. 

It turns out, as we shall justify in the following, that the appropriate type of 
convergence is weak convergence. Specifically, we shall say that a cluster size 
distribution c(m,t) tends to a scaling form if there is a function s(t) — > oo as 
t — > oo and such that 

oo oo 

J mc(m,t)f[m/s(t)]dm t — > W J x$(x)f(x)dx (3.12) 





for every continuous bounded function f(x) on the positive real axis. Note that 
the use of x$(x) in the definition is purely conventional: it serves to make 
notation consistent with the previous subsection, as well as with ordinary 
usage. 

We shall argue in the following that (3.12) is the most appropriate definition 
for the appproach to scaling. On the one hand we shall find that the various 
exactly solved models which we discuss in Section 4 yield counterexamples to 
many attempts to make sharper statements, whereas (3.12) is, as we shall see 
in some detail, fulfilled in the vast majority of cases. 

Let me here make a technical — but quite important — remark: in order to 
prove weak convergence of measures, it is necessary and sufficient, by a stan- 
dard result in probability theory [114], to show that the corresponding Fourier 
transforms converge. Since all the measures we consider are concentrated on 
the positive real axis 10 , this is equivalent to convergence of the correspond- 
ing Laplace transforms, or generating functions in the discrete case. In other 
words, if we define 

oo 

G((;t) = J c(m,t)e< m dm (3.13) 
o 

we obtain the following two equivalent formal expressions for the scaling func- 
tion <&(x): 



W ( $(x) (e~ px - l) = lim {s(t)G[-p/s(t),t}} 



= W-'C- 1 {lim [G p (-p/( 8 (t),t)]} (*), (3-14) 



10 Well, almost all. When we look into the distribution of compositions in multi- 
component aggregation, we shall in fact use Fourier transforms 
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where C~ l denotes the inverse Laplace transform. Since generating functions 
such as G((,t) are ubiquitous in exact solutions, where they arise in a very 
natural manner, we may already argue from an aesthetic viewpoint in favour 
of the above defiinition of convergence to scaling. 

Finally note that there is one element of arbirariness in the definition (3.12) 
that I have not yet addressed: it concerns the choice of the measure we decide 
to consider. In (3.12) we chose mc(m,t), partly because it is automatically 
normalized. It is, however, quite possible to use, for example, m n c(m, t)/M n (t) 
in its stead. We shall say that the n'th moment of the cluster size distribution 
approaches a scaling form if there is a function s(t) as above such that 

oo oo 

Mnity 1 J m n c(m,t)f[m/s(t)]dm^W J x n $(x)f(x)dx (3.15) 



for every continuous bounded function f(x) on the positive real axis. Now, it is 
straightforward to check that, if the n'th moment of a cluster size distribution 
tends to a scaling form, then so do all higher ones. However, the opposite 
need not be the case. Consider for example the exact solution to (2.2) for 
K(m,m f ) = mm' and an initial condition c(m, 0) = 5(m — 1), which is given 
by 

oo 7 fc-2 (f p -t\k 



and is valid for t < 1. In this case, one readily verifies that, choosing s(t) = 
(1 — t)~ 2 , the second moment of (3.16) approaches a scaling form with $(x) 
given by x~ 5 ^ 2 e~ x , whereas the first does not tend to a scaling form at all. 
Generally speaking, the scaling function to which the n'th moment tends must 
satisfy 

oo 

J dxx n <S>(x) < oo (3.17) 



so that a singularity of the type x~ 5 ^ 2 could never arise as a limit of first 
moments. 



3.3 Some Consequences that Follow from Scaling, and More that Don't 



We now need to understand which consequences follow from the scaling hy- 
pothesis. In many systems, it is possible to measure some of the moments 
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M p (t) of the cluster size distribution c(m,t), as defined in (3.10). In fact, it is 
often, at least in principle, possible to reconstruct the full distribution from a 
knowledge of its moments. We therefore concentrate first on these. From the 
definition of convergence to a scaling form (3.12) follows 



lim 



s(t) _(p-1) M p (t)] =W [ x p $(x)dx = WM p , 



(3.18) 



whenever the integral on the r.h.s converges. Here the last equality defines 
the scaled moments M p which play a considerable role throughout the theory. 
As we shall find out later, <&(#) has no singularities for x ^ and decays 
exponentially as x — > oo. 

The only issue is therefore the behaviour at zero. In many cases, this is given 
by a power law as follows 



$(#) ~ x 



(s->0). 



(3.19) 



The meaning of such a behaviour at the origin is the following: it states that 
the concentrations averaged in the low end of the cluster distribution go as 
m~ T . More formally, 



lim 



S(tf 



1 



2es(t) 



es(t) 



c(m' , t)dm' 



0(e- T ) 



(3.20) 



as e — > 0. Here the l.h.s. is an averaged concentration rescaled by s(t) 2 so as 
to cancel the time dependence. 



It now follows that a necessary condition for (3.18) to hold is 

P > T — 1. 



(3.21) 



(3.18) and (3.21) therefore confirm the claim made in (3.11) that s(t) can 
be defined as the ratio of successive moments, if the order of the moments is 
sufficiently large. 

Let us now introduce some notation: denote by z the exponent with which 
s(t) grows as a function of t and a p the one with which M p (t) grows 



s(t)~t z 
M p (t)~t a ? 



(3.22) 
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From these definitions we obtain the following relationships 
ap = (p-l)z (p>r-l). 



(3.23) 



We may further define w as the exponent with which c(m, t) decays when s(t) 
becomes much larger than m. Formally once more 



lim 



2es(t) 



es(t) 



c(mf, t)dm 



es(t) 



o(r w ) 



(3.24) 



From (3.19) and (3.24) one readily shows the following celebrated identity 
[110] 



(2 — t)z = w. 



(3.25) 



This relation relates large time behaviour of small clusters to the shape of 
the cluster size distribution at the low end of the distribution for fixed large 
times. Since these two quantities are not related in a straightforward way, this 
relation often gives very useful and non-trivial information. 

Note, however, that the above definitions of r and w seem quite artificial. It 
appears much more natural to ask about the large time behaviour of c(m, t) at 
fixed m and the small m behaviour at fixed t. Here, however, it must clearly be 
stated that these asymptotic behaviours cannotbe determined from the scaling 
hypothesis in the form that we have given it up to now. Similarly, moments 
M p (t) of order lower than r — 1 cannot be described by scaling theory as it 
stands. In fact, their behaviour is linked to that of the c(m,t) at fixed m for 
large t, which lies beyond the reach of the scaling approach. We may define 
some additional exponents: for example, define r' and w' as follows 



c(m, t) ~ m T ' (1 <C m « 

c(m,t)~t- w ' (*->oo). (3.26) 

At this stage, however, I should again remark that exactly solved models are 
known (see section 4) for which the w' exponent depends on m, whereas the 
t' exponent does not exist. These, however, satisfy ordinary scaling and have 
well defined values of the exponents r, z and w. 

Nevertheless, since such systems are not very frequent (they are not, however, 
pathological), it makes sense to say that a system satisfies the strong scaling 
hypothesis if r' and w' both exist and satisfy 

T ' = T w' = w. (3.27) 
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Under these circumstances, we can say how the moments M p (t) of order p < po 
behave if they are finite initially: the a p are given by 



ol. 



(p — \)z {p > t — 1) 

( P <T-1). 



(3.28) 



w 



Note, however, that if strong scaling is violated, it may well happen that 
a p 7^ {p~ I)- 2 5 even when p > r — 1, as we shall later show in explicit examples. 
The only generally true statement is that a p = (p — l)z if p > 1. Since the 
evaluation of moments is of central importance in applications of the theory 
to real systems, this result is basic. It shows that the large time behaviour of 
moments can always be obtained to a knowledge of the exponents r and z for 
given forms of the reaction kernel K(m,m'). It is this task to which we now 
turn. 



3.4 The Scaling Equation and Its Derivation 



The basic kinetic equations (2.2) together with the definition (3.8) of the fact 
that the first moment of c(m,t) approaches a scaling limit, given in (3.12) 
lead to the following condition on the scaling function $(x): for all continuous 
functions on the positive real axis which vanish sufficiently fast at infinity, one 
has, as shown in Appendix A: 



dxdyK(x,y)<Z>(x)®(y)x[f(x + y) - f(x)\ = / dx x 2 f(x)$(x). (3.29) 



Here K(x,y) is the continuous function defined by 



K(x,y) = lim s x K(sx 1 sy). 



(3.30) 



As part of the derivation of this equation, one obtains a condition on s(t), 
namely 



s(t) = Ws(t)\ 



(3.31) 



the solution of which is in agreement with the statements obtained on quite 
general grounds in (3.8). This therefore confirms that 

z = T ^-. (3.32) 
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when gelation does not occur. The derivation of (3.29) is somewhat technical 
and is given in appendix A. We present it in some detail, though, as the scheme 
is very straightforward and can be extended in a mechanical way to all the 
other cases we shall be considering in this paper. 

A particularly useful form of (3.29) is the one in which f(x) is taken as an 
arbitrary exponential: 



oo 

I(p) ; = J dx dyK(x,y)<S>(x)<S>(y)xe- px [l - e~ m ] (3.33) 
o 

oo 

I(p)=p J dxx 2 e- px ^(x). (3.34) 



Here I{p) is defined for future reference in (3.33) A very similar equation is 
central to the work of van Dongen [108]. It is obtained from (3.29) by setting 
f(x) = Q(a — x) for arbitrary a and yields the relation 



a oo 



a 2 $(a) = J dx J dyxK{x,y)<5>{x)<5>{y) (3.35) 



a—x 



I present (3.35) because of its importance in previous work, but shall not make 
much use of it in the sequel. From (3.29) also follows immediately the following 
central expression for the scaled moments M p . 

oo 

pM p+1 = J dxdyxK(x,y)$(x)$(y)[(x + y) p - x p ] (3.36) 
o 



When K(x,y) is of the form x a y x ~ a , then (3.36) becomes a nonlinear recur- 
sion between the scaled moments whic is frequently useful, see for example 
subsection 3.5.2. Also, as we shall later show, the scaled moments enter in the 
determination of r for certain cases (see (3.48) below), and (3.36) then turns 
out to be very useful to get interesting results from these connections. 

Finally, let me make some remarks on uniqueness and normalization of the 
solutions of (3.29) as well as all its special variants. It is clear that (3.29) inher- 
its the symmetry of the equations (2.2), as stated in (3.7). Namely, whenever 
$(x) is a solution, so is 

S 6 (x) = 6 1+A $ (£) (3.37) 
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for all b > 0. This ambiguity in the definition of arises from the absence 
of scale in the equations (2.2) at least when aggregates are large enough and 
differs from the ambiguity which simply results from the choice of the scale of 
s(t), which we have absorbed in the prefactor W in (3.12). It should generally 
be checked that final results are not affected by the transformation (3.37). 
To fix the normalization, we use the fact that, since x<&(x) is defined as the 
weak limit of a probability distribution, it must be normalized (gelation is not 
assumed to occur in this case). This means 

oo 

J x$(x)dx = 1. (3.38) 
o 



This condition then determines b in 3.37) and hence <$>(x) uniquely. Finally, 
the absolute scale of <&(#) is related to that of K(x, y), which is itself fixed by 
the choice of time scale. This can be stated as follows. If <&(#) is a solution with 
rate K(x,y), then 6$(x) is a solution to the equation with rates K(x ) y)/b. 

Note that a slight difficulty arises when A = 1, since the transformation de- 
fined by (3.37) then leaves / x<$>(x)dx invariant. In this case we must keep the 
constant W\ explicitly in (A. 5) and set the normalization by allowing it to 
vary. 

The above relations were all derived under the hypothesis that the first mo- 
ment of the cluster size distribution tends to a scaling form. This is not always 
the case. In fact, it is easy to verify that, when r > 2, it is impossible to have 
such convergence. Therefore, in the gelling case, which always has r > 2, we 
must consider convergence of higher order moments. In this case one finds 



J dx [(3 - r)f(x) + xf'(x)] x 2 §(x) = j dxdy K(x, y) x 

x$(x)$(y)x [{x + y)f(x + y) - xf(x)} . (3.39) 

Here we have used the strong scaling assumption to show that 

M 2 (t) V J s(t) K ! 



The appearance of the exponent r in this equation may at first appear discon- 
certing, but there is no real problem: (3.39) can be solved for any r. It then 
remains to see for which value of r the condition (3.19) obtains. We now show 
how these equations can be used to determine r from the knowledge of the 
exponents A and ji characterizing the reaction kernel K(x,y). 
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3. 5 The Solutions of the Scaling Equation 



3.5.1 The Non-gelling Case: Small-x Behaviour 

A very useful mathematical concept in the determination of the behaviour of 
$(x) near the origin is that of regular behaviour. We say that a function <J>(x) 
has regular behaviour at the origin if the limit 

A (a) = limiM (3.41) 



exists. If this is the case, it is easy to see that a number r must exist such that 
h(a) = a~ T . (3.42) 

This is shown by proving that h(a) must satisfy the functional equation 

h(ab) = h(a)h(b), (3.43) 



for which power laws are the only solutions 
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We now assume that has regular behaviour near the origin. That is, 

we assume that it has a power law behaviour, possibly dressed by arbitrary 
logarithmic corrections, but no behaviour of the type cos(61na;). We shall find 
that this assumption is self-consistent for p > 0, but that it is not so when 
p, < 0. In the former case it allows us to derive more or less explicit expressions 
for z and r. In the latter, we shall see that $(x) has a stretched exponential 
behaviour near the origin, which can be determined by other techniques. It 
should be realized, however, that self-consistency does not prove that the 
function <3>(x) does display regular behaviour at the origin. In fact, there is 
substantial numerical evidence [56] that when p > the function &(x) has 
non-regular behaviour at the origin, with (possibly) undamped logarithmic 
oscillations. 

Since we are interested in small values of x, we naturally consider the large-p 
behaviour of the two sides of (3.34). We show in appendix C, see (C.6) and 
(C.7), that 

I(p) = B ip -^^{l/pf + B 2 p-^(l/p)k(l/p), (3.44) 



where B\ and B 2 are definite integrals defined by (C.7) with the following 
peculiarity: of the two integrals involved, one will diverge and the other one 

11 1 disregard non-measurable functions throughout this paper. 
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converge for any choice of Only the term with the convergent integral 

should be retained in (3.44). 



From this follows, by matching I{p) with the r.h.s. of (3.34) which is of order 
p~ 2 $(l/p), the following well-known three cases (the numbers have become 
standard terminology): 

• Case I: /i > 0. For these values of the parameters one finds 

r=l + A w=l. (3.45) 

In this case, however, there exists some doubt on the very existence of a 
well-defined solution to (3.29). The issue is the following: it is readily seen 
that the function x~( x+1 ^ solves (3.29) quite generally. However, it clearly 
never fulfils the essential normalization condition (3.38). This means that 
the condition of strong decay as x — » oo is essential in determining the 
solution. This has cast some doubt on the existence of such a solution. On 
the other hand, recent numerical work by Lee (see [56]) strongly suggests 
that such a solution does exist 12 , though it may not behave regularly at 
the origin: specifically, Lee observes oscillations in 3>(x), which might well 
be of the type cos In x. These do not conform to the regularity hypothesis. 
On the other hand, if $(x) really does behave regularly at the origin, we 
may state that 

= B lX - T , (3.46) 

where Bi is given by (C.7), which is a definite integral that can be explicitly 
evaluated for every single kernel. 

Finally note that the value for w is the lowest possible for the monomer 
decay exponent in a system with binary reactions, since c\{t) can always 
decay at least via a reaction with itself. 

• Case II: /i = 0. In this case, the exponent r is non-universal, but satisfies 
the inequality 

t < 1 + A. (3.47) 

If $(x) behaves regularly at the origin, then a comparison of the prefactors 
of the r.h.s. and l.h.s. of (3.29) yields the well-known relation [108] 

oo 

t = 2-T j x x <$>(x)dx = 2- r¥ A , (3.48) 
o 

where T is defined by (3.4) and M p is defined as in (3.18). Note that this 
result is invariant under the symmetry transformation (3.37) as well as under 
a common scale change of $ and K (due to the presence of V). 



12 contradicted, however, by earlier work, see in particular [53] 
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Of course, (3.48) does not straightforwardly determine r, since it requires 
in principle the knowledge of the full scaling function <&(x), from which the 
exponent r could have been derived without recourse to (3.48). However, 
as has been shown in [108,18], (3.48) can be used to determine bounds on r 
for any specific kernel independently of $(#). In particular, Cueille and Sire 
[18] have developed a large number of very elegant methods which allow to 
obtain extremely sharp rigorous upper and lower bounds on r, to the extent 
that one may view these as the equivalent of an exact determination. In some 
specific cases, r can also be displayed as the solution to a transcendental 
equation [102]. 

Case III: fj, < 0. In this case, the assumption of regular behaviour of $(x) 
near the origin leads to inconsistencies. We may therefore assume a decay 
of $(x) near the origin that is faster than every power. One may then 
approximate I{p) as shown in appendix C, with the result 

oo 

I{p) = TMx-^J x 1+ ^(x)e~ px dx. (3.49) 
o 

From this and (3.34) follows, by inverting Laplace transforms, that within 
our approximations 



x $(x) 



rMA^x^^fi). (3.50) 



From this one finds the approximate form for <$>(x) near the origin in this 

case: 



x 2 &(x) = const. ■ exp 



x -\iA 

H 



(3.51) 



Here only the exponential behaviour with the constants indicated are reli- 
able. The power law x~ 2 does not follow at the degree of accuracy we have 
used here. For a detailed treatment of corrections to scaling for both Case 
II and Case III, see [108]. It is an entertaining excercise for the reader to 
verify that (3.51) is indeed invariant under the symmetry transformations 
(3.37) as well as under a simultaneous change of scale in K and <&(#). 

Yet another dificulty arises when A = 1 and /i > 0: in this case, since r is 
predicted by (3.45) to be two, the normalization condition (3.38) for <3>(x) is 
not satisfied any more. Under these circumstances the easiest solution is pre- 
sumably to view this as an instance where it is more appropriate to look at the 
convergence of the second moment to a scaling form, so we defer consideration 
of this issue to the next subsection, in which we consider also the gelling case, 
which cannot be treated in any other way. For a different approach, however, 
see [108]. On the other hand, if // < 0, we find ourselves either in Case II 
or in Case III. Due to inequality (3.47) in Case II and (3.51) in Case III, no 
convergence problems arise even when A = 1. 
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3.5.2 The Non-gelling Case: Large-x Behaviour 



The large-x behaviour in the non-gelling case has been treated by van Dongen 
and Ernst in [106,108]. We present here a slightly modified version of their 
approach. Let us first show that whenever K{m,m') is less than to + to', the 
scaling function $(#) decays exponentially for large values of x. This can be 
done as follows: consider the moments M p of the scaling function. These satisfy 
the relation (3.36), which implies, using the bound on the kernel K(m,m') 



M p+1 <^-J2 r \M k M p . k+1 . (3.52) 



2p 



k=i 




If we now denote by the moments which satisfy (3.52) as an equality, it 
is easy to show inductively that 

M P <M ( P °\ (3.53) 



for all p. On the other hand, it is seen that the M p °^ can be evaluated analyti- 
cally, since they are the moments of the scaling function of the kernel to + to', 
which is known exactly. One therefore obtains 

M p < Hf = const. ■ p~ 3/2 p\R p (3.54) 

for all p and some fixed value of R, which strongly depends on various nor- 
malizations. (3.54) now immediately leads to the desired conclusion by using 
standard results on the moment problem. Similar lower bounds can also be 
shown assuming correponding lower bounds on K(m,m'). 

Assuming some such lower bound, it follows that the integral equation (3.34) 
also holds for negative values of p, up to a critical value — p c beyond which all 
integrals diverge. Introducing 

h{x) = <$>(x)e pcX (3.55) 
one obtains the following equation from (3.34) for h(x) 



oo oo 

dxx 2 h(x)e~ ip+pc)x = -p J alx dy xK{x,y)h{x)h{y)e- {p+pcKx+y) x 

x [1 - e py ] (3.56) 
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If we now assume that h(x) behaves regularly at infinity with exponent 
that is, 

.. h(ax) .0 , 
hm -V^ = a (3.57) 
x ^°° h[x) 



then one obtains, using altogether the same techniques as for the small- a; 
behaviour, that if v < 1 

■& = A. (3.58) 



If, on the other hand, v — 1, then the integral on the l.h.s of (3.56) diverges. 
A treatment entirely similar to the one corresponding to the indeterminate 
Case II for the small x behaviour shows that the exponent d in this case also 
cannot be determined uniquely using only the type of scaling considerations 
discussed above. Details on this subject are to be found in [106]. 



3.5.3 The Dominant Singularity Hypothesis 

The following scenario for the explanation of scaling may, at this stage, quite 
possibly have occurred to the thoughtful reader: if we consider the generating 
function (or the Laplace transform) of the cluster size distribution function 
G((,t) defined by 

oo 

G(£,t) = J c(m,t)e< m dm, (3.59) 
o 



one may assume that the appearance of a unique divergent size s(t) in the clus- 
ter size distribution c(m, t) might translate itself in the appearance of a single, 
dominant singularity of G((,t) on the positive real axis. Certainly, standard 
theorems on Laplace transforms state that the nearest singularity does indeed 
lie on the positive real axis, so that such a hypothesis may seem plausible. 
This assumption might be further strengthened by the careful study of a large 
number of exactly solved models, since indeed the mechanism explaining the 
growth of the typical size and the appearance of scaling is invariably of this 
nature. 

Let us now use our understanding of the nature of the scaling function <&(#) 
in order to assess this possibility. Put in equations, our assumption states: 

G((,t) = const. ■ [(c — (s(t)]® +1 + higher order terms, (3.60) 



where ( c is some appropriate constant. 
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I have denoted the relevant exponent by d because, as we shall now see, it has 
always that value. Indeed, we know from (3.14) that 



The large x behaviour of is given, as is well-known from general theorems 
on Laplace transforms, by that of the nearest singularity of Cr(C, t). This near- 
est singularity is on the real axis, and we have assumed it to be a power-law 
singularity. It therefore certainly always describes the value of the d exponent. 

Does it also yield the value of the r exponent? Clearly, this depends somewhat 
on the meaning we give to the phrase "higher order terms" in (3.60). However, 
generally speaking, if we have assumed only one singularity dominating the 
whole picture, it is hard to see how it could be otherwise. One therefore deduces 
from the hypothesis just described that 



But this, we know from subsection 3.5.2 is simply wrong in the large majority 
of cases. Thus, when /i > 0, r = 1 + A, which cannot be equal to i?, which is 
A. The only cases in which this possibility exists are the loci fi = 0, for which 
r is non-universal, and v — 1 for which d is non-universal. It is certainly 
striking that all the exactly solved models known to this day find themselves 
on these loci, and that they indeed satisfy (3.62). This should therefore serve 
as a warning: some features of exactly solved models are highly atypical. 

In subsection 4.10 we shall see that the constant kernel with production term 
provides a counterexample to the dominant scaling hypothesis. In this 
will be shown later, one has r = 3/2 and ■& — 0. The generating function has 
infinitely many simple poles tending simultaneoulsy to the origin as t — > oo. 
The nearest such pole dominates the large mass behaviour and leads to $ — 0, 
whereas the combined effect of all the poles together upon values of ( just 
below zero leads to a —3/2 singularity. 

3.5.4 The Gelling Case: Scaling Region 

In the gelling case, we define scaling using the approach to scaling of second 
moments. This is appropriate as long as r < 3. In order to have no difficulty 
on that score, we limit ourselves in the following to the following range of 
parameters 




(3.61) 



r = ■&. 



(3.62) 



1 < A < 2 



v < 1. 



(3.63) 
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The condition A > 1 is necessary to ensure gelation, whereas the condition 
v < 1 is required to avoid pathologies such as instantaneous gelation or non- 
existence of solutions. The condition which we really only put in for conve- 
nience is A < 2. 

Our definition of scaling is therefore the following, there exist a typical size 
s(t) diverging at the gel time t c and a scaling function <£>(a;) such that 



lim 

tyt c 



1 



M 2 (t) 



m 2 c(m, t)f 



m 



At), 



dm 



oo 

W J x 2 <$>(x)f\x). (3.64) 



One now proceeds as in the non-gelling case to derive an integral equation for 
The details are carried out in Appendix F. If we now define 

s(t) = const. ■ (t c - ty 1/a , (3.65) 

It then follows from the work done in Appendix F that 

a = 1 + A - r. (3.66) 

In other words, it is no more possible to fix the exponent describing the growth 
of the typical size in terms of the exponents describing the kernel. Instead, we 
have a scaling relation between a and r. Finally, for the integral equation 
determining $(x), one finds 

oo 

J dx dyxK(x, y)&(x)&(y) [(x + y)f(x + y) — xf{x)\ = 
o 

oo 

J dxx 2 <S>(x) [(3 - r)f(x) + xf'(x)} (3.67) 
o 

The novel feature here again is, of course, the appearance of r in the equation. 
At the basic level, presumably, it does not change much: one can, in principle, 
solve (3.67) for each value of r and then see for which values of r (if any) 
the corresponding solution, which we may call 3> T (x), behaves as x~ T near the 
origin. 

The difficulty of a theoretical treatment, however, is thereby considerably 
increased. Let us first ask whether a formal solution of pure power-law form 
exists, akin to the £-( 1+A ) solution in the non-gelling case. Let x~ T be such 
a solution. In this case, it follows using partial integration that the r.h.s. of 
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(3.67) vanishes identically for arbitrary functions. Hence, putting f(x) = e x , 

oo 

J dx dy K(x, y)x x ~ T y~ T e~ x [(x + y) e - y - x] = 0. (3.68) 



It is a straightforward exercise in gamma functions to show that, for K(x,y) 
of the form x p y v + x"y^, the solution of (3.68) is 

r = ^ (3.69) 



Since, however, every kernel of homogeneity degree A can be expressed in 
terms of kernels of this form, it would appear to follow, at least on a formal 
level, that this value of r is the one that corresponds to the universal power- 
law solution for arbitrary K(x,y). As we shall see in the next subsection, this 
value is indeed highly significant: it is the only large-size behaviour compatible 
with a finite flow rate of mass to infinity. As such, it is the only possible large- 
mass exponent for times beyond the gel point. However, we are interested in 
the scaling region immediately before the gel point. There is, in principle, no 
reason to assume that this particular value of r plays any role in this region. 

Let us now attempt a scaling theory of the type we have achieved for the 
non-gelling case. Our results, unfortunately, will be much more fragmentary 
than in the regular case. As before, we use exponentials e~ px as functions f(x) 
and investigate the large p behaviour. Let us assume that <3>(x) is of the form 

= Ax- T \l + Bx A + o(x A )] (3.70) 



where we have taken next-to- leading behaviour of <&(x) into account, since 
this is what determines the behaviour of the r.h.s. of (3.67). One then obtains, 
under the hypothesis that the integral (3.68) is different from zero, after some 
tedious but straightforward calculations 

A = l + X-r = a. (3.71) 



On the other hand r cannot be determined by straightforward scaling argu- 
ments. Thus we have a general relation between the first correction to scaling 
and the growth of the typical size, but neither can be determined, both de- 
pending on the exponent r which must be determined by other considerations. 
It is not clear to me whether similar techniques to the ones that were success- 
ful in determining r to very high accuracy in Case II of the non-gelling case 
[18,19] may be generalized to this case. 

There are therefore two cases to be distinguished: 
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(1) The ordinary case: The integral (3.68) is different from zero. In this case 
we know nothing about the actual value of the exponent r, but we know 
how both o and the leading next order correction depend on it. The 
positivity of A, which follows from its definition, leads to the inequality 
for r 



(2) The extraordinary case: The integral (3.68) vanishes. In this case we have 
the "standard" value (3.69) for r. The computation of the correction to 
scaling exponents is somewhat more difficult, whereas the exponent a is 
given by 



The nomenclature used above is meant to suggest that the ordinary case is 
more common than the extraordinary. While I have no analytic arguments to 
show it, the following evidence exists: first, Lee [56] has studied a large number 
of gelling kernels of the form x^y v + a^y M , and has determined r numerically 
to high accuracy. In no instance except the case /i — v — 1 was r given by the 
"standard" value (3.69) but was always smaller. Second, even in the bilinear 
kernel, it can be shown (see Appendix L) that the case of power-law initial 
conditions leads to continuously variable values of r as well as an exponent 
A satisfying (3.71). Therefore, even the bilinear kernel sometimes yields the 
ordinary case. On the other hand, the extraordinary case really does occur for 
the bilinear kernel with initial conditions decaying rapidly with mass, since 
then A = 1, which contradicts (3.71) with A = 2 and r = 5/2. 

In principle, it would now be necessary to derive the value of r in the ordinary 
case from the overall structure of the reaction constants K(m,m'). As far 
as I can see, this problem is of considerable difficulty, and I have made no 
headway at all. Possibly the methods developed in [18,19] could be generalized 
to yield sharp upper and lower bounds on the r exponent. This is, of course, 
particularly important since all other exponents are expressed in terms of r 
and are not known explicitly. Lacking such a theory, we are led back to the 
numerical study either of the gelation transition in the original equation (2.2) 
or else of the scaling function $(x) in (3.67). 

Finally, let us note that the ordinary case can be characterized in an amusing 
way as follows: consider the average concentrations defined as follows 



2 < r < 1 + A 



(3.72) 



a = 



A - 1 
2 



(3.73) 




(3.74) 
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It is then a straightforward consequence of the definition of scaling and of the 
exponent r, that 

C(e) = const. ■ [es(t)]~ T (3.75) 



Now, it is easy to verify that the value of the correction to scaling exponent A 
given by (3.71) corresponds exactly to the fact that this limiting behaviour is 
reached linearly in time. This only holds in the ordinary case, so that it should 
be violated in the case of the bilinear kernel with initial conditions that decay 
rapidly as a function of mass, as indeed it is: there the average concentrations 
are readily seen to tend to their limit quadratically in time. 



3.5.5 The Gelling Case: After Gelation 

After gelation sets in, the power-law behaviour of the cluster size distribution 
remains for all times. Since, by definition, mass must continually decrease after 
gelation, the derivative must be a finite number, that is, 

d f 

< lim — / mc(m,t)dm < oo. (3.76) 
M->oo dt J 



But one finds 



M M 

d 



dt 

o 



J mc(m,t)dm — J dm\ J dm 2 Q(mi + m 2 — M) x 
o o 

xmiK(mi,m2)c(mi,t)c(m 2 ,t) (3.77) 



If one now inserts the power-law profile m Ts for c(m,t) in (3.77), one finds 
that 

d 1 

— J mc(m, t)dm = const. ■ M x+3 ~ 2Ts (3.78) 



from which follows that, in order for (3.76) to hold, t s must indeed have the 
"standard" value discussed in the previous subsection. 

Concerning the large-time behaviour of the concentrations, it can be shown 
rather generally that they all go as In the discrete case, this may be seen 
as follows: consider first the monomer equation 

oo 

c 1 = -c 1 J2K(hk)c k ~-c 1 M u (t). (3.79) 
fc=i 
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Since r — v > 1, the sum in M u (t) converges, so that this quantity behaves 
similarly to ci(i). Equation (3.79) then reduces to 



ci = —const. • c 1; 



(3.80) 



which has the stated behaviour. The general statement then follows by induc- 
tion, since it is straightforward to show that the production terms are of order 
t~ 2 . Since, however, the removal terms must dominate, they must be of order 
t~ 2 as well, implying the stated result. 

In fact, it is seen that the ansatz 



t + t 



(3.81) 



can be made to satisfy (2.2) if the ctj satisfy the following algebraic relations 



"3 2 



Lfc=l 



-1 .7-1 



J2 K (k,j ~ k)a k aj_ k . 
k=i 



(3.82) 



It turns out to be quite hard to know whether these equations indeed have 
a positive solution. For the case in which K(k,l) is of the form (kl) x ^ 2 , con- 
siderable progress can be made by reducing (3.82) to a recursion together 
with a self-consistency condition. It was shown [60] that in this case the co- 
exist and that for large j they display a power-law behaviour j~ Ts if A > 1. If 
< A < 1, these solutions presumably describe the behaviour of Cj(t) at fixed 
j and t large. They behave for large j as j _( - 1+A - ) , in agreement with the results 
derived previously for Case I. 



3.5.6 The N on- Gelling Case for A = 1 

As has been pointed out above, if A = 1 and v < 1, the ordinary scaling theory 
suggests r = 2. This, as first pointed out by van Dongen and Ernst [108], is a 
contradiction, since one then has a divergent total mass at the origin. The way 
out is to consider this as a borderline case and treat it via the convergence in 
second moment, in the same way as we have already done with gelation. 

The details are carried out in Appendix F and lead to the following relations 



s(t) = const. ■ exp {const. ■ y/tj 
«i(«) = (3.83) 
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Convergence to scaling then takes place as follows 



lim 



OO OO 

Vte'^ J m 2 c(m,t)f (me^) dm = j x 2 <$>(x)f \x)dx. (3.84) 



One may then absorb one factor of m into the function f(me ^*) and obtains 



lim 

t^oo 



oo oo 

\ft j mc(m,t)f |me"^j dm = J x<&(x)f(x)dx 



(3.85) 



The function <£>(x) must then satisfy the following integral equation derived 
in Appendix F 



dx dyxK(x, y)§(x)§(y) [(x + y)f(x + y) — xf(x)} 



oo 

J dxx 2 §(x) [f(x) +xf(x) 



(3.86) 



As in the gelling case, 3>(x) need not have finite mass. This means that, as in 
the gelling case, the total mass contained in the scaling region vanishes. 



3.6 Crossover 



Up to now, we have mainly considered the case in which the reaction rate 
K(m,m') are homogeneous in the masses. The justification of this hypothesis 
lies, as stated in the Introduction, in the fact that we are principally inter- 
ested in the behaviour of the system at large times and for large aggregate 
sizes. Arguably, in such circumstances, any simple physical mechanism for ag- 
gregation will have a simple behaviour with respect to scaling of the masses 
in this asymptotic regime. On the other hand, if two mechanisms act jointly, 
we may expect that the one that scales with the lesser homogeneity degree 
will be altogether negligible, so that we may, again, limit ourselves to the 
homogeneous case. These considerations already suggest, however, a possibly 
important exception: Assume a fast mechanism which only becomes appre- 
ciable at quite large sizes, superimposed on another mechanism, with a lesser 
degree of homogeneity, but effective at all sizes. The following illustrates such 
a possibility: 



K(m,m') = K[ Xl \m,m') + eK ( 2 A2> {m,m') 



(A 2 ), 



(3.87) 
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Here the Ai^ are two different degrees of homogeneity and e is a very small 
number. Under these circumstances, if A2 < Ai, the second term in the sum can 
obviously be neglected, but in the opposite case, it will eventually dominate. 
From the point of view of pure scaling thery, it is in principle sufficient to 
say this, and one has, indeed, a valid prediction for "sufficiently long" times. 
For practical purposes, however, this is often an excessive limitation, and one 
would like a scaling theory which gives correct results when times are large 
with respect to the scale determined by K± 1 as opposed to e _1 . To this end, 
we need a theory of crossover. 

To avoid unnecessary difficulties, we first fix the scale of the parameter e: 
indeed, it is obviously irrelevant whether we use e or /(e) as a crossover pa- 
rameter, where / represents an arbitrary monotonic function vanishing at the 
origin. In order to fix the parametrization, I therefore make the following con- 
vention 

K(am, am'] ej a) = a x K(m, m; e) (3.88) 



At a naive level, this means that e has a dimension of inverse mass, that is, the 
natural scale on which e varies is given by s(t) _1 . This requirement does not 
determine A uniquely, though: indeed, replacing K by K' = e a K changes A to 
\ + a while respecting (3.88). We therefore impose the additional requirement 
that K(l, 1) always remain of order one. Thus, for example, the kernel given 
in (3.87) is rewritten as 

K(m,m) = K[ Xl) (m,m) + e X2 - Xl K { 2 X2 \m,m) (3.89) 



Note that this only gives a meaningful result if Ai < A2, as was to be expected 
on intuitive grounds anyway. 

We now proceed as always: let us first define what we mean by a scaling limit. 
This is given by 

OO , X OO 

Urn / m'c(m',y/s(t);t)f I J -> W J x<5>{x,y)dx, (3.90) 



where c(m, e; t) denotes the value of the concentration of clusters of mass m 
at time t evolving under a dynamics using a kernel with a fixed value e of the 
crossover parameter. We now use the standard approach in order to derive the 
following equation for the typical size and the scaling function. For the typical 
size s(t) one obtains 

s(t) = Ws(t) x , (3.91) 
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where A is defined by (3.88). Note that, for the case described by (3.89), 
this means that the degree of homogeneity is Ai, that is the lesser degree of 
homnogeneity. From this follows that we will be looking at the size distribution 
on a scale defined altogether by the slow aggregation process defined by K[ Xl \ 
This may at first appear strange, but it should be noted that, in the scaling 
limit (3.90), the crossover parameter e goes to zero, so that the large-time 
behaviour of any specific system is not the one observed in the crossover limit. 

For the scaling function, we obtain the following equation in exactly the same 
way as we obtained (3.29): 



oo 

J dx x dx 2 x 1 K{x 1 , x 2 ; y)$(x u y)<S>(x 2 , y) x (3.92) 
o 

oo 

x [f( Xl + x 2 ) - /(xi)] = J dxxf(x) [y$ y (x,y) - x& x (x,y) -2$(x,y)], 



which must hold for all continuous function f(x) and for all y. From this we 
again derive an equation involving exponentials similar to (3.34): 



oo 

J dx! dx 2 x x e- pxx K{x x ,x 2 ; y)<&{x 2 , y) x (3.93) 



o 



X 



1 - e~ pX2 



dxxe px [2§(x,y) +x<5> x (x,y) -y%(x,y)}, 



If one now takes the limit of large p and makes exactly similar considerations as 
in the case without crossover, it is straightforward to show that the exponent 
r defined by 

lim $^) 



where y is taken to be fixed, is the same as for the kernel K(xi,x 2 ; 0), that is, 
it is identical to the value which would be established if the fast mechanism 
dominating at long times did not exist. In other words, the small size end of the 
cluster size distribution is not affected by the existence of the fast aggregation 
mechanism in the crossover limit as defined by (3.90). Again, this is not at 
variance with the prediction that the whole size distribution will eventually 
be determined by the fast process: as time goes on, at fixed e, the value of y 
diverges and the range of x for which the exponent r of the slow process yields 
a good description, tends to zero. 
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The general properties of (3.94) are not easy to work out, however. It is seen 
that there is a formal solution of the form 



<S>(x,y) = y x+1 X (xy), (3.95) 



where x( x ) is some function satisfying an equation which is readily obtained 
from (3.94). This equation, however, turns out to have no meaningful solutions, 
since setting p to zero in it leads to a contradiction. The solution given in (3.95) 
must therefore be viewed as a purely formal one, similar to the solution of the 
form $(x) = £ _ ( A+1 ) in the usual theory. 

Frequently, in practice, one is interested in the time evolution of the moments 
of the size distribution function, since these are often the only accessible quan- 
tities. In this case, if the system is at a fixed value of e, the resulting expression 
is given by 

oo 

M p (t; e) = W j x p $(x, es(t))dx, (3.96) 
o 

where s(t) is given by the solution to (3.91). Note, of course, that under these 
circumstances it is not possible to define s(t) via moment ratios, as is otherwise 
often done. 



3. 7 Extensions of the Model 



3. 7. 1 Monomer Production and Diffusion 

In many situations, it happens that monomers are continuously supplied to 
the system, so that total mass is not conserved any more. If this lasts only for 
a short time, the evolution of the system after the monomer injection can be 
viewed as ordinary aggregation starting from a somewhat polydisperse initial 
condition. This case is therefore not essentially new. If, on the other hand, 
the supply of monomer takes place in such a way that the total mass supplied 
diverges, new phenomena arise. 

Let us specifically look at the following equations also considered by various 
authors, see in particular [73,74] and references therein, as well as [20,57]: 



d t c(m,t) = J dmi dm 2 K(m 1 ,m 2 ) c(mi,t)c(m2,t) x 

x [5{mi + m 2 — m) — 5(mi — m) — 5(m 2 — m)} + 
+p{t + ty5{m-m ). 



(3.97) 
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Here the time shift to is introduced merely in order to avoid spurious diver- 
gences if u is negative. From (3.97) one immediately obtains for the total mass, 
assuming that gelation does not take place, that is, that A < 1, 

dmmc(m,t) = \ + J ■ ( 3 - 98 ) 



From this follows the following natural analogue to the scaling definition (3.12) 

on on 

1+u 



pm (t + ty+" 



mc(m,t)f[m/s(t)]dm — > W / x<b(x)f(x)dx (3.99) 





From this we may again, using methods entirely similar to those employed in 
Appendix A, obtain the following result for the typical size s(t): 

s(t) « (Wtyv+W 1 -*). (3.100) 



For the scaling function $(x) one obtains the following relation similar to 
(3.29): 



1+LO 



oo oo 

dxx®(x)f(x) + J x 2 ${x)f'(x)dx — J dxdy K(x,y)x$(x)$(y) x 



W 

x[f(x + y)-f(x)] + W- 2 f(0), (3.101) 





-2 . 



If one now substitutes f(x) by e px one obtain in the limit of large p the 
following condition for $(x): 

I(p) = W~\ (3.102) 



where I(p) is defined by (3.34). If we now use the results of Appendix C on 
the large p asymptotic behaviour of I(p), we find 



T = < 



(u<u+l) 

p + 2 (v>p+l) 



Assuming strong scaling we can state how the small clusters behave. If v < 
p, + 1 one finds 

Cj (t) » j-&+Wf»l* (3.104) 
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as stated in [20]. In particular, if oo = 0, that is, if monomers are injected into 
the system at a constant rate, one approaches a stationary distribution with a 
power-law exponent (A + 3)/2. That this value is ubiquitous when a stationary 
state is reached is well kmown (see e.g. [113]). On the other hand, if v > /j + 1, 
one finds 

Cj{t) « r (M+2) t (M+l-(-+lM/(l-A) (3 1Q5) 

which leads to decaying solutions when oo — 0. Therefore, in this case, con- 
stant monomer input does not lead to a stationary cluster size distribution, 
although no gelation takes place. A case where this in fact occurs is given, 
for instance, by the kernel mj~ 2 + 2 ■ In this case all concentrations decay 
as t -1 / 3 irrespective of the power law oo, that is, independently of the rate at 
which monomer is fed into the system. The systematic nature of the approach 
developed here is worth emphasizing: while I believe most of the results stated 
in this section to be well-known, I am not aware of a similarly straightforward 
unified derivation. 

3.7.2 Spatially inhomogeneous systems 

Let us now look at spatially inhomogeneous systems. As it frequently happens 
that monomer injection occurs at a well localized position, such systems are 
of genuine interest. At first, one might argue that they must lie outside the 
domain of validity of mean field theory. This is not generally the case, however: 
In order for the mean field approximation to hold, one requires that there be 
no correlations at the molecular level. Therefore, we must certainly exclude 
cases in which the density varies significantly on the scale of typical interag- 
gregate distances. It often happens, though, that one has homogeneity over 
such scales and yet a slow spatial variation giving rise to diffusive dynamics 
being superimposed upon the reaction process. This is the situation we shall 
consider. For closely related work, see e.g [17], but the literature on this kind 
of problems is considerable. 

We are therefore led to the following kinetic equations 

oo 

<9 t c(m,r) = | J dm\ dm 2 K[m\, m 2 )c(mi, r)c(m 2 , r) x 
o 

x [8{mi + m 2 — m) — 5{mi — m) — 5(m 2 — m)] + 
+D(m)A r c(m,r) + p5(r)5(m - m Q )(t + t) u} . (3.106) 

In the scaling approach we now introduce two time dependent growing quan- 
tities, namely s(t) which is the typical size of aggregates, and L(t), which the 
typical distance over which aggregate concentrations vary significantly. We 
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then require as a consistency condition on our approximations, that L(t) be 
much larger than the radius of a typical aggregate. As we have discussed in 
the Introduction, the relation betwen mass and radius is by no means a trivial 
issue and depends sensitively on the particular model. However, in a large 
variety of cases, the asymptotic relationship between the mass s(t) and the 
radius R(t) can be described by 

s(t)^R(t) D f, (3.107) 



where Df is the fractal dimension of the aggregate and lies between one and 
the space dimension d. The minimal requirement for the applicability of mean- 
field therefore becomes 

L{t) > s(ty/ D f (3.108) 



for al t. In particular, if the growth exponent of L(t) is less than z/Df, then 
mean-field will at best be applicable during a finite range of times, and non- 
mean-field behaviour will dominate the system at large times. In the opposite 
case, van Dongen has determined the limit of validity of mean-field theory by 
considering the size of the fluctuations and the efficiency of their transport. If 
one assumes that the diffusion constants D{m) decay with an exponent m™ 7 , 
where 7 > 0, it is found that the critical dimension d c above which mean-field 
holds, at least qualitatively, is given by 

d c = - \ (3.109) 

1 — A — 7 



if the r.h.s. is positive, and is infinite otherwise. 



3.7.3 Scaling Theory for Inhomogeneous Systems 

Let us develop a scaling theory for (3.106). We limit ourselves to the case 
in which the injection occurs at the origin and we transform to spherical 
coordinates in order to exploit the existing symmetry. As definition of the 
scaling limit we choose 

00 00 , . 00 

Jdmm J drr d - x c{m,r;t)f ^yj J dx x<Z>(x,y)f(x,y) (3.110) 



Here s(t) denotes, as usual, a typical size and L(t) a typical length. 
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We use the same procedure as always and find for the scaling function <3>(x, r) 
/i(p) = / 2 (p)+/ 3 (p) + l, (3.111) 



where the various Ij(p) are defined as 

oo 

Ii(p,a) = J dxidx2drr d ~ 1 K(xi.X2)xi&(xi.r)Q(x 2 .r) x 
(l - e~ pX2 ) 

I 2 (p } a) = a 2 J dxdrr d - 1 x 1 -"<§(x,r)e- e7r M - ^— ^ J (3.112) 
o \ ar / 

oo 

7 3 (p, <r) = y &drr' i - 1 2:$(i) [1 + u-zp- vo\ e - px - aT '. 
o 

Here z and z/ stand for the exponents with which the typical size and length 
grow with time respectively: 

s(t) = const. ■ t z L{t) = const. ■ t v (3.113) 



and their values are found as in the usual case during the derivation of (3.111) 
These are 



o 

x ( '" ! ' 

oo 



Z 



uj + 2 
1- A 



v = \{\- 1Z ). (3.114) 



The last relation is physically easy to inetrpret if one remarks that it is equiv- 
alent to 



L(t) ~ JD [s(t)\ t. (3.115) 



The fundamental length scale is therefore the one on which a typical cluster 
diffuses in time t. 

What can one deduce from (3.111)? These are still quite formidable equations, 
and it is not clear how to get results out of them. The following approach 
yields some results: consider an arbitrary exponent a > and look at the 
curve defined by a = p a , where both p and a go to infinity. If we further 
assume that 

lim $ [-,- J = Ms)p T{a \ (3.116) 
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we find straightforwardly the following orders of magnitude at large p and a 
for the Ij{p): 



I l (p,a) = p- x - 3+2T ^' da 

L,(p,a)=p^ d)a ^ +T{a) 
I 3 (p,a) = p- 2 - da+ ^) 



(3.117) 



Matching these and taking also the constant term into account one eventually 
finds the following relation between r(a) and a: first, if a > —7/2, then we 
define 



Otc. = 



1 - 2 7 - A 
A-d 



(3.118) 



Note that we shall always limit ourselves to the case d < 4, since otherwise it 
follows from well-known results that the aggregation process becomes irrele- 
vant and the whole system reduces to a non-interacting diffusion. 



One then finds for r(a) 



±(cfa + A + 3) 
1 + A + 2a + 7 



(a < a c ) 
(a > a c ) 



(3.119) 



and similarly if a < —7/2, we find a different value of a c 
1-A 



a, 



d 



(3.120) 



and the values of r(a) are given by 

A + 1 (a < a c ) 

^ \{da + A + 3) (a > a c ) 



r(a) = < 



(3.121) 



What do these results mean? A dependence r(a) suggests a singular behaviour 
at the origin of the type 



(f>(x, r) ~ exp 



In a; • r 



lnr N 
In a; 



(3.122) 



For the sake of comparison with known results, consider the case of the con- 
stamt kernel with constant diffusion. We are therefore always in the first case, 
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since a > 0. The singularity structure predicted by the theroy is therefore 

r) ~ x- 3/2 r-- d/2 + aTV -2 (3.123) 



and the critical line which separates the regions in which one or the other 
summand dominates is given by the relation r = x 1 ^ 4- **). The exact connection 
of these results with the related ones published in [17] is not entirely clear to 
me, but it is certainly very close. There it is claimed that in the stationary 
case, the scaling variable is m/r 4 ~ d and that the r exponent is given by (6 — 
d)/(A — d), which is exactly the value of r(a c ) in our formalism. 



3.7.4 Higher Order Reactions 

Let us now consider the case in which reactions of higher order than two may 
also be present. The reaction term of order r is given by 



— [ drrii ■ . . . ■ dm r K r (mi, . . . , m r )c(m 1 , t) ■ . . . ■ c(m r , t) x 
r\ J 

x 8(mi + . . . + m r — m) — 5(m r i — m) 



(3.124) 



Denoting (in this subsection only) the degree of homogeneity of the kernel 
K r by A r , we obtain by a formal reasoning entirely analogous to that used in 
Appendix A the following conditions on the typical size s(t): 



s(t) = W^sit) 



r-l „( + \X r -r+2 



(3.125) 



At this stage, however, I must emphasize that this derivation is not rigorous for 
reactions of order r, when r > 3, for reasons explained in detail in Appendix 
B. Let us nevertheless assume it for the time being. 

If therefore two different orders r\ and r 2 are relevant in the same physical 
process, we might expect that one will invariably dominate the other in the 
scaling limit, unless 

A n - K 2 — ri - r 2 . (3.126) 



As we shall see later, however, the situation is more complex. The crucial issue 
is (3.125). As stated above and set out in detail in Appendix B, this derivation 
is subject to certain objections which do not arise in the two-body case. It 
is perfectly possible for a scaling limit to exist, with a well-defined scaling 
function <&(x) satisfying some form of an equivalent scaling equation, and yet 
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for (3.125) not to hold. Under such circumstances, which we will discover in 
certain exactly solvable models, a more complicated situation arises, which 
we shall not discuss in the general framework of scaling theory because of 
its complexity. We shall fully discuss the exactly solvable case, however, in 
subsection 4.8. 

Let us now look more carefully at the pure case of three-body reactions. Let us 
first describe the various exponents describing a three-body kernel. We define 
A3, /ii and /i 2 as follows 

K(x,y,z) = z X3 k 2 (x/z,y/z) 

k 2 (y, z) = z x *h(y/z) (y, z « 1) (3.127) 

k 2 (y,z) = T 1 (y)z^ {z « 1). 



Note that, since the arguments of K(x, y, z) can always be written in ascend- 
ing order, it is enough to specify k 2 (x,y) on the triangle < x < y < 1. 
Conversely, k 2 (x,y) can be chosen arbitrarily on this domain, so that no lim- 
itations on the values of the above exponents hold. From this follows that 
(3.127) gives all the relevant information concerning the possible limiting be- 
haviours of K(x, y, z). In particular, one readily obtains the following limiting 
behaviour for T-\_(x): 

T 1 (x) = const. ■ x X2 ~» {x -> 0) (3.128) 



Note that (3.127) is not the most general behaviour imaginable: it would be 
possible, for example, to define k 2 (y, z) in such a way that its decay exponent 
ji as z — > depended on the value of y. We disregard such cases and leave 
their treatment to the interested reader. 

We now proceed quite similarly to the case of binary reactions and obtain an 
equation for the scaling function <&(x) of the following form 

00 

I 3 (p) = p J x 2 ®(x)e- px dx, (3.129) 



where Iz(p) is given by 

J 3 (p) = lj dxdydzK 3 (x,y,x)<Z>(x)<Z>(y)<Z>(z)xe- px [l - e - p{y+z) ] (3.130) 



A fundamental difference between this case and the two-body reaction case 
should again be observed: in the derivation of the scaling equation for 
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in the two-body case, we saw in appendix A that the behaviour of the small 
aggregates does not couple to that of the aggregates in the scaling regime: in 
other words, we did not need to invoke anything beyond scaling to justify the 
scaling equation (3.29). On the other hand, in the case of (3.130), a problem 
may arise from the three-body reactions involving two large particles and one 
small one. The technical issue involved is discussed in detail in Appendix B 
for the three body case. As we shall see in exactly solvable models, such cases 
actually occur and signal the necessity of choosing a different value for the 
growth exponent z than the one predicted by (3.125). This is a qualitatively 
new phenomenon, which does not occur in the two-body case. 

The large-p behaviour can again be estimated as in appendix C. This leads to 



I 3 (p) = B lp -^ +4 ^(l/p) 3 + B 2 p-^ +3 ^(l/p) 2 + 

+B 3 p' 2 ~^(l/p). (3.131) 

Here Bi are three integrals given in appendix D. One now proceeds as in the 
binary reaction case. The following cases arise, excepting always the possibil- 
ity that a divergence at the origin occurs leading to a different value of the 
exponent z: 

(1) Case I: p > and A3 — 2\ 2 < 0: This essentially corresponds to Case I 
in the binary situation. One finds r = 1 + A3/2 and the corresponding 
value of 1/2 for the decay exponent w for the small aggregates, which is 
the lowest possible value for w in the purely ternary case. For this to be 
consistent we additionally need the condition A 3 < 2p. If this is violated, 
we have a divergence at the lower end of the integral and a different value 
for the exponent z. 

(2) Case II: p > and A3 — 2A2 > 0: This case is new. One finds r = 1 + A2, 
and correspondingly 

w = ±^ (3.132) 

In this case, the second term in (C.4) dominates. Again, for consistency 
we need A2 < p, otherwise we again have divergence at the lower end of 
the integral and a modified value of the exponent z. 

(3) Case III: p = 0: This case corresponds to Case II for binary reactions. In 
this case, r may take any value within a certain range: in order for there 
to be no divergence at the origin, we require r < 1, and on the other 
hand, in order for the third term in (C.4) to dominate one needs 

r < min(l + A 3 /2,l + A 2 ). (3.133) 

(4) Case IV: p < 0: In this case, we cannot have a finite value of r. If we 
have normal behaviour, then the function $(:r) vanishes at the origin 
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faster than any power law. 



We shall find the above theory indeed to be violated in the following exactly 
solved cases: 



Ki has A3 equal to 1 and A2 and jj, equal to zero, but the z exponent is equal 
to 2 and r = 3/2 13 . For K 2 we have A3 equal to 2, A2 equal to one and \x 
equal to zero. In this case we find gelation with a divergence exponent z of 
—2 and r = 2.5, that is, behaviour quite similar to the product kernel in the 
binary case. These results will be shown in detail in subsection 4.8. 

3. 7. 5 Multicomponent Aggregation 

An interesting generalization of the aggregation model consists in allowing the 
mass parameter m to become a vector [69,51]. The meaning of the components 
could be, for instance, the concentration of various species present in the 
aggregates. Similarly, we can think of one species as being charged and the 
other neutral, so that the reactivity of an aggregate depends not only on its 
total mass, but on the total number of monomers of the charged species. The 
reaction rates K(mi,m 2 ) then depend on two additive parameters and are 
therefore an instance of the possibility of multicomponent aggregation. 

Apart from its relevance in realistic aplications (for an interesting use of re- 
lated concepts concerning the coagulation of charged aggregates, see for ex- 
ample [39]), the real interest in studying this case comes from the possibility 
of looking at a different kind of scaling behaviour: Indeed, so far, we have 
looked only at the average size of the aggregates. Here, as I shall show shortly, 
under quite general circumstances the limiting distribution consists of aggre- 
gates which have a fixed composition determined by the masses of the various 
species initially present. A natural question then follows: on what scale does 
the distribution of compositions vary? It is less than the typical size, so we 
obtain a scaling in which two different sizes play an important role. In fact, 
we show that this distribution is quite generally a Levy stable distribution 14 . 
Further, it can be shown that the index of the Levy distribution is given by the 
ratio of the exponents characterizing the growth of composition fluctuations 
to the exponent z for the growth of the typical cluster size. These results are 
in fact to be expected: indeed, we may consider the following stochastic model 

13 As if we were in fact dealing with Case I. Is this more than a coincidence? I make 
no claim to know 

14 possibly an ordinary Gaussian. 




(3.134) 
(3.135) 



K 2 (mi,m2, m 3 ) = m^j + m 2 m 3 + mim 3 
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for aggregation. Take a finite but large set of aggregates of various composi- 
tions. Then choose two at random according to the rates i^(m 1 ,m 2 ) and join 
them. Let this process be repeated indefinitely. It is then quite likely that an 
appropriate generalization of the law of large numbers will apply, so that all 
aggregates tend to a fixed composition. If one then wishes to analyze further 
on what scale the composition varies, one needs some version of the central 
limit theorem. (Similar random coagulation processes have been looked into 
in [82], but not, it seems to me, to the extent of proving such claims.) Such 
a theorem would state that the resulting distribution is always Levy and that 
the typical size of the fluctuations is connected to the index of the Levy law 
involved [8]. Note that the constant kernel corrresponding to such a case was 
already studied in [69,51]. Defining 



\m\ 



J2 m * ( 3 - 136 ) 



the scaling assumption now reads 

hm J dm\m\c(rh,t)f = J dx\x\&(x)f(x). (3.137) 



A*), 





Using exactly the same techniques as in the one component case, one derives 
the following equivalent of the scaling equation (3.29): 



J dxx-Vf{x)x — j dxidx2K(xi : X2)Q(xi)Q(x2) x 

x[/(f 1 +f 2 )-/(f 1 )] (3.138) 

In order to simplify notation, I limit myself to the case of two components, 
though the general case presents no difficulties of principle. A moment's thought 
shows that the following is a solution of (3.138): 

${x) = $ 1 (x 1 )6(x 1 -x 2 ), (3.139) 



where Xi denotes the i-th component of x and $i(x) is the scaling function cor- 
responding to the rates K(xi,Xi, x 2 , x 2 )- There exist further solutions differing 
from (3.139) by trivial normalizations, but since we can always fix units so as 
to have the same mass of both components initially present, this is irrelevant. 
I do not know whether further solutions exist, but the plausible considerations 
made at the beginning of this subsection concerning the law of large numbers 
appear to speak against such a possibility. 

Let us now look at the variations in composition and the way in which they 
scale. To fix ideas, we start from a system in which the mass initially contained 
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in both components is the same, and we assume that the scaling function given 
in (3.139) is indeed the relevant one. This permits us to neglect the dependence 
of the reaction rates on composition, since a fixed composition is reached at 
the scale defined buy the typical size s(t). . We then define a scaling function 
^(cr, 5) for the sum and difference between x\ and x 2 . We define it as follows: 



oo 

lim / dmidm 2 (mi — m 2 )c(m, t)f 

t-+oo J 



17li + 177-2 77li — WL2 

. Slit) ' s 2 (t) , 



oo oo 



J da J d66V(a,6)f(a,6). 



(3.140) 



-oo 



Here si(t) corresponds to the typical size of aggregates, whereas s 2 (t) gives 
the scale for the variations in composition. In Appendix E. we derive a scaling 
equation (E.l) for ty(a,5) and show that 

M,( a} 5)=a- a <S> l (a) Xa (^) (3.141) 



is a solution. Here <3>i(x) is the one-component scaling function of (3.139), 
Xa(x) is the symmetric Levy stable distribution of index 1/a given by the 
Fourier transform of exp(— |g| 1/,Q! ) with respect to q and a is given by 

a = lim (3.142) 

t-oo ^(t) S 2 (t) 



In order to determine a unambiguously, we argue as follows: All Levy distri- 
butions except for the Gaussian have power-law tails. Further, if a > 1, it is 
known that the first moment of the distribution diverges. This can certainly 
not happen if, as we have always supposed, we start from initial conditions 
having finite mass. It is also known that the Xa are only positive if a > 1/2. 
We therefore have 1/2 < a < 1. The value 1/2 corresponds to a gaussian dis- 
tribution, whereas all other allowable values correspond to a power-law decay 
of x~ x ~ x l a . Thus a value of a different from 1/2 implies that the differences in 
composition have power-law tails extending to infinity in the large time limit. 
Since no mechanisms are known to produce these unless they be present in the 
initial conditions, we may state that if the initial conditions do not have power- 
law tails, the final distribution must be Gaussian 15 . If the initial distribution 
did contain a power-law tail, then it is quite likely that the corresponding 
Levy stable law will be the relevant one in the scaling regime. We therefore 

15 This is maybe a little less obvious than it seems: as we shall see later in detail, 
the large x asymptotics in the scaling regime has no simple relationship to the large 
m asymptotics at fixed times. It is, however, a very reasonable assumption. 
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see that initial conditions are relevant in this case to determine the asymptotic 
behaviour of the typical size of the variations in the aggregate composition. 
We shall later find similar behaviour in ballistic aggregation, which is, in fact, 
a closely related problem. 

3.8 Non-scaling Asymptotics 

The above is a subject which, from the peculiarly biased viewpoint I have been 
taking in this paper, will naturally receive rather short shrift. Nevertheless, it 
is rather obvious that the following questions are relevant: 

(1) How does c(m,t) behave at large times for fixed ml 

(2) Conversely, how does c(m,t) behave for large m at fixed tl Obviously, 
the decay should, in some sense be exponential, or at least fast enough, 
but more detailed information concerning power-law corrections to the 
exponential decay is also important. 

(3) What is the intermediate asymptotics of the c(m, t)? By this I mean the 
following: frequently, there exists a time dependence f{t) such that 

c(m,t) = g(m)f(t)[l + o(l)] (t - oo) (3.143) 

for fixed m. In this case, one may ask what is the large m behaviour of 
g(m). In the absence of such a function f(t), the whole issue becomes 
rather ambiguous. 

(4) Similarly, the problem of the large time behaviour of the moments of the 
cluster-size distribution is of considerable physical relevance, since these 
are usually its most easily accessible characteristics. The problem splits 
into two disjoint issues. 

First, the large time behaviour of moments M p (t) of order p such that 

oo 

J x p $(x) = oo, (3.144) 



falls outside the purview of standard scaling theory. The problem, of 
course, lies in the fact that these low-order moments are dominated by 
the low-m part of the cluster-size distribution, whereas the higher-order 
moments are determined by the distribution as a whole. 

On the other hand, high-order moments may also cause difficulties, 
though perhaps of a more trivial kind: indeed, if the initial condition 
contains a power-law contribution, then some sufficiently high-order mo- 
ment will diverge from the very beginning., which is a clear sign that we 
may not blindly apply the scaling approach in evaluating moments. This 
is seen, for example, in the exactly solved model of the constant kernel, 
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for which I show that scaling is approached whenever the initial distribu- 
tion has finite first moment. Nevertheless, it is obvious that the second 
moment has no well-defined large-time asymptotics if its initial value is 
divergent. One therefore sees that stronger assumptions than merely the 
existence of the scaling limit are required in order to justify the usual 
scaling expressions even for high-order moments. 

In the literature, such issues have often been (mistakenly) discussed as so- 
called "violations of scaling" (see e.g. [63,64] for a characteristic example of 
such confusion. The issue of scaling in parity dependent kernels will be taken 
up again later in this paper.) In the following, we shall always refer to our 
standard definition of scaling, and will only consider as non-scaling such sys- 
tems as do not approach a limiting form in the sense defined by (3.12). In 
fact, as we shall see, all of the quantities described above may exhibit un- 
expected behaviour while the scaling hypothesis remains fulfilled. We shall, 
in our terminology, call such deviations from expected behaviour, violations 
of strong scaling. Since I have never been too explicit about what is meant 
by strong scaling, it is clear enough that anything peculiar enough to attract 
some attention might be brought under this heading. In the following, I shall 
largely follow the approach sketched in [62] 16 . 

Let us first turn to the large time behaviour of c(to, t) at fixed to. In this case, 
production of clusters of mass to can be neglected, and I may approximate 
(2.2) by 

oo 

c(m,t) — —c(m,t) J K(m,m')c(m' ,t)dm' . (3.145) 
o 



Using the definitions of /i and v given in (3.4) and (3.5) as well as the fact 
that the relevant masses to' in (3.145) are much larger than to, one eventually 
obtains the following approximate differential equation for c(to, t) 

c(to, t) = const. ■ m»M v {t). (3.146) 



In order to understand the behaviour of this system, we must first know the 
behaviour of M u (t). Assume strong scaling, after which one may check for 
self-consistency. In this case, one finds for kernels of type III 

M„(t) = const. ■ f("-i)/(i-A) > (3.147) 



6 The reader should be warned that this paper thoroughly confuses the two concepts 
of scaling regime and non-scaling asymptotics. 
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from which follows that c(m, t) decays as a stretched exponential of the form 



c(m, t) = const. ■ exp —const. ■ mH A - ) 



(3.148) 



Note that this is in perfect agreement with the predictions of strong scaling, 
so that no violations are to be expected in Case III. Of course, this might not 
be viewed as quite conclusive, since we need a hypothesis to determine the 
large-time behaviour of M u (t). As we shall see, however, the main hypothesis 
is that the expression of the moment as an integral of the scaling function 
should have no singularities at the origin, which is trivially satisfied in Case 
III. Another source of problems might arise from power. law initial conditions, 
but again, as long as v < 1 no divergence is possible. 

In case I, the above approach fails, since it is readily seen that M„(t) cannot 
be described by a convergent integral of the scaling function. The following 
observation, on the other hand, yields interesting results [67]: if one introduces 
the following new variables in the discrete representation, which turns out here 
to be more convenient, 



<f> j (e)= Cj (t)/c 1 (t) 

d6 = Cl (t)dt, (3.149) 
the original equations (2.2) become 



dfa ;_J 



J 3 



\ ]T K(k, j - k)fafa- k - fa ]T [K(k, I) - K(l, I)] fa. (3.150) 

k=l k=l 



These equations have a non-trivial equilibrium solution aj which can some- 
times be determined recursively. It can be worked out for kernels of the form 
K(j,k) = (j/c) A / 2 ,(see [61]), as well as more generally [103] that these con- 
stants behave as j~ T with r = 1 + A, which is therefore in good agreement 
with the qualitative behaviour predicted by strong scaling. This also settles 
the issue of intermediate asymptotics as defined above, showing that it indeed 
exists, and that it coincides with the asymptotic behaviour of small clusters 
in the scaling limit. 

One further finds using (2.2) that Ci(t) goes as from which follows that all 
Cj(t) do. This demonstration has been put on a firmer footing by a rigorous 
inductive proof by van Dongen and Ernst [103]. On the other hand, there 
remain discrepancies between the behaviour of the aj at large j and that of the 
scaling function 3>(x) at small x: in particular, at least under the assumption 
of regular behaviour for <&(#), the prefactors of both differ. 
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Finally, for case II, no such arguments are available. In particular, neither an 
approach based on (3.146) nor one based on recursion relations can be justified. 
Further, a large number of counterexamples are known, which precisely belong 
to this case. In fact, both non-gelling exactly solved classical models, the 
constant and the sum kernels, belong to case II, as also do some variations 
on the constant kernel for which the strangest deviations from strong scaling 
have been reported. 

It remains to discuss shortly the issue of the large m behaviour at fixed time 
[25]. In this case one may safely neglect the reaction of particles of mass m 
with larger aggregates, so that we have the following approximate equations, 
for the discrete case: 

j'-i 

d 3 = lYl K ( k iJ ~ k )°kCj-k- (3.151) 

k=l 



Again, for monodisperse initial conditions, this can be solved in a straightfor- 
ward way by recursion, using the ansatz: 

Cj(t) = ajt j -\ (3.152) 



One finds (for details see [25]) that Oj goes as j~ A for all kernels K(k, I) with 
v < 1. A comparison with the results obtained for the large x end of $(x) 
shows that a difference between the two can, and in fact does occur for the 
case v — 1 [106]. It therefore follows from these results that the large mass end 
of the cluster size distribution at fixed times always decays exponentially, a 
fact which we shall find to be of some importance in assesssing the possibility 
to describe certain models by mean-field equations. 

Let us finally consider the scaling prediction for the moments. Overall, the 
strong scaling prediction for the asymptotic behaviour of M p (t) is 

M p {t) = Hpsitf- 1 (t -> oo) (3.153) 



The immediate case in which this cannot hold is, of course, for values of p such 
that M p does not exist (diverges). Clearly, this can only occur when <&(#) has 
a power-law singularity at the origin, since the exponential decay of <&(#) at 
infinity shown in subsection 3.5.2 precludes any divergence there. Thus there 
are no problems in Case III, but whenever r is finite, the asymptotic behaviour 
(3.153) is invalid in the range 

p<-l + T. (3.154) 
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However, as we shall see later in exactly solved examples, this behaviour can 
also fail outside the range stated in (3.154). Indeed, in various kernels of type 
II one finds that for large t at fixed m 



c(m,t) =r u, ' (m) [l + o(l)] (t 



oo). 



(3.155) 



Here w'(m) varies continuously with m and bears no necessary relation to 
either z or r. In this case, it may well happen that (p — l)z be less than 
max m w/(m), so that the decay of M p (t) is dominated by that of a specific 
type of aggregate. We shall in fact see such examples in the next section. This 
definitely cannot happen, however, for moments which grow in time, that is 
whenever p > 1. 

An altogether different violation of the behaviour stated in (3.153) occurs 
whenever the initial conditions have a power-law tail. In this case, as one sees 
in many exactly solved cases, the large mass behaviour always maintains this 
tail, which cannot be eliminated by any aggregation process except gelation. In 
this case, clearly, whenever p is large enough to cause the divergence of M p (0), 
the moment M p (t) is identically infinite, so that (3.153) becomes meaningless. 
Nevertheless, scaling often holds also under these conditions. In such cases, 
we must reinterpret (3.153) to mean 

As(t) A 



where A is an arbitrary constant. The time necessary for (3.156) to become 
valid goes to infinity as A does. 



3.9 The Moment Equations 

The moments M p (t) of the cluster size distribution function satisfy some gen- 
eral exact equations, which are in many cases of considerable use, though we 
shall not make much explicit use of them in this paper. Indeed, from (2.2) 
follows immediately 




(3.156) 



o 



o 








x [(mi + m 2 ) 



v 



p 




(3.157) 
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For integer values of p, these take a particularly striking form: 



dM p 




dmidm2rn\m2 k K(mi, m 2 )c(mi, t)c(m 2 , t). (3.158) 



The two principal uses of these relations are the following: first, in the case 
of the bilinear kernel (4.68), they yield closed equations, which can be solved 
recursively. Their solution for low values of p is one of the quickest ways to 
obtain qualitative information concerning the solution. In particular, it readily 
yields the exponent z as well the gel time in gelling systems. 

On the other hand, they can also be used for estimating the solution from 
above or below. Thus White's proof of global existence for the solutions of 
(2.2), see [112], rests essentially on such an estimate which follows from the 
estimate 



The possibility of such estimates is a result of the positivity of the right-hand 
side of (3.158). One then obtains a closed set of differential inequalities, which 
can then, because of their recursive character, be solved to yield upper bounds 
on the c(m,t). Similarly, we used inequalities of the same nature in subsec- 
tion 3.5.2 to estimate the high moments of the scaling function $(x). Note 
carefully that the moment equations (3.158) are not at all equivalent to the 
corresponding equations (3.36) for the moments of the latter only de- 

scribe the moments if the system approaches scaling, whereas the former hold 
in every case. In fact, the knowledge of all moments is equivalent, in principle, 
to knowing the full cluster size distribution, so that the solution of (3.158) is 
in principle of the same degree of difficulty as the full solution of the original 
Smoluchowski equations (2.2). The solution of (3.158) therefore also contains 
indications of all corrections to scaling and non-scaling asymptotics, which, 
as we have seen, often differ considerably from the scaling behaviour. These 
two equations cannot, therefore, be identified. Nevertheless, it is certainly true 
that the leading behaviour of M p (t) for p sufficiently large only depends on 
the scaling function $(#). It might therefore be possible to prove the scaling 
hypothesis by a careful study of the leading terms for the large-time behaviour 



K(m, m) < C(m + m) 



(3.159) 



of (3.158). 
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4 Exactly Solved Models 



In the following, I shall discuss a large variety of exactly solved models. One 
can justify interest in such cases in several ways: they are in themselves a 
source of pleasure, they are an ideal test ground for numerical work and they 
allow us better to understand the structure of solutions in general. While I 
sympathize to some extent with all of these motives, in this paper I shall be 
primarily concerned, as I have been throughout, with the issue of the valid- 
ity of scaling theory. Exact models are, of course, invaluable in this respect, 
since we have ordinarily a complete overview of the solutions behaviour. It is 
therefore straightforward to say in which way the system does ineed satisfy 
the scaling hypothesis, and to what extent it may occasionally deviate from 
what is expected. As we shall see, the scaling hypothesis as I have defined it, 
see (3.12), holds essentially for all systems we shall have occasion to look at 17 . 
On the other hand, as I shall show in an extensive set of examples, the set of 
assumptions I have denoted by the general name of "strong scaling", are in 
general not fulfilled. To be more accurate, it is frequently possible to find more 
or less natural exactly solved cases for which these assumptions fail in several 
ways. The detailed illustration of such failure I consider to be instructive, since 
it shows with great detail exactly how much we may, and more importantly 
what we may not, deduce from the existence of a scaling limit. 

The exactly solved kernels fall in two categories. First, one might mention the 
so-called "classical" kernels 



which are known as the constant, sum and product kernels respectively. Their 
solution can be generalized to the case of an arbitrary linear combination of 
the three, see (4.68), which is known as the general bilinear kernel. The three 
kernels described in (4.1-4.3) are all quite different from each other, as should 
be expected from the general scaling theory exposed in the previous section: 
the first has A = 0, the second has A = 1 and the third has A = 2 and is 
accordingly a gelling kernel. In the first two cases, on the other hand, we have 
regular growth, proportional to time in the first case, exponentially fast in the 
second. All these predictions are indeed borne out by the exact solutions. 



The solitary exception I have in mind is the crossover between the sum kernel and 
the product kernel, for which I have not been able to find a scaling description. How- 
ever, this may still be a soluble problem. In any case, this only concerns a crossover 
phenomenon, for true scaling behaviour I am not aware of any counterexamples 





K 3 (m, m!) — mm! , 



(4.3) 
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As for the scaling theory, these kernels in general show a reasonable behaviour, 
that is, they usually satisfy scaling in an extremely strong form. However, 
even in these cases, it can be shown that many issues depend rather sensi- 
tively on initial conditions. We shall see, in particular, that initial conditions 
with power-law tails at large masses can significantly modify the non-scaling 
behaviour, even when the scaling function is unaffected. 

I shall also discuss some natural extensions of these kernels to ternary reac- 
tions. This will show the way in which the scaling approach may actually break 
down altogether if collisions involving particles outside the scaling regime dom- 
inate the physics. 

Further, since not only the three kernels above can be solved, but also the 
general bilinear kernel (4.68), we can clearly study various types of crossover. 
In fact, as I shall show, two of the three possibilities can be analyzed rather 
extensively and one arrives eventually at a crossover function as defined in 
subsection 3.6. 

Finally, there are non-classical kernels, such as the g-sum or the parity-dependent 
kernels. These are characterized by the fact that they are not exactly homoge- 
neous. As we shall see, they display a considerable variety of unexpected be- 
haviours in their non-scaling asymptotics, but all satisfy scaling in the sense 
specified in this paper. Furthermore, as we shall see, even these satisfy the 
Dominant Singularity Hypothesis, that is, the value for the small-size expo- 
nent r is the same as that of the large size exponent ■&. Nevertheless, the array 
of non-scaling anomalies shown in these cases is instructive, as it shows us 
that certain features one takes for granted are, in fact, peculiar features of the 
simpler exactly solved kernel and do not generalize to more general cases. 



4-1 The Constant Kernel 
The constant kernel defined by 



was first solved by von Smoluchowski [89] and its solution has been often 
rederived. The approach we follow is an elementary, yet instructive, use of the 
generating function approach: define 



K(m 1 , m 2 ) = 1 



(4.4) 




(4.5) 
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Here the integral must be replaced by a sum whenever c(m, t) is a sum of delta 
functions. The fundamental equation (2.2) now takes the form 



F t = ±F((,t) 2 -F((,t)F(0,t). (4.6) 



After making the replacement 

G(C,t) = F(C,t)-F(0,t), (4.7) 

(4.5) becomes 

G t = -G 2 
2 

oo 

G((, t = 0)= g(() := J c(m', 0) (e< m> - l) dm'. (4.8) 

o 

From this one readily obtains the following explicit form for G((,t) in terms 
of g((), which in turn is explicitly determined by the initial concentrations. 
The connection between the actual concentrations at time t and the initial 
values is, as in many other cases, quite opaque. We shall see, however, that 
for many relevant issues the specific nature of this connection is not needed. 
This expression for G((,t) is 



If, using the notations introduced in (2.3) and (2.4), I wish to compute Cj(t) 
for monodisperse initial conditions, it is a straightforward excercise in power 
series manipulation to show that 

4 / t 

<*(*) = ttt™ (rr^J • ( 4 - 10 ) 



(t + 2) 2 \t + 2 



This is readily seen to approach a scaling form: if we namely set s(t) = t, we 
obtain using the definition (3.12) that (4.10) yields as scaling function 

$(ar) = Ae~ 2x . (4.11) 



It therefore follows that r = 0, z — 1 and w — 2. This manifestly satisfies the 
scaling relation (3.25) derived in [110]. To tidy up various details, we point 
out that (4.11) is indeed a solution of (3.34): all functions of the form ae~ x l 2 
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satisfy it for K(x,y) = 1 and only (4.11) satisfies the normalization condition 
(3.38). Furthermore, setting s(t) equal to t, as we did, corresponds because 
of (A. 5) to setting W equal to one. Everything therefore confirms the scaling 
calculations down to the smallest detail. Note that, in the following, we shall 
usually not verify the normalization issues in great detail. Rather, we shall 
usually implicitly assume W — 1, which involves a certain choice of time scale 
such that s(t) satisfies excatly s = s x , with no further prefactors. 

Now we still need to show that this scaling form holds for arbitrary initial 
conditions. In itself, this is non-trivial: estimating the asymptotic behaviour 
of the formal expressions for c(m, t) for arbitrary initial conditions is quite 
cumbersome, even in this most elementary case. We therefore rely on an ob- 
servation made at the end of subsection 3.2: namely, in order to show scaling, 
it suffices to show convergence as defined by (3.14). In our case, it amounts to 
showing that 



exists. But this is quite straightforward using (4.9), if we may assume that 
g(() has finite first derivative at £ = 0, or, what amounts to the same, that 
the c(m, 0) have finite first moment. Under this hypothesis, one finds that the 
limit defined in (4.12) is given by 



which is, according to (3.14), the Laplace transform of x&(x) evaluated at — s. 

This system provides various interesting examples for several of the phenom- 
ena discussed in Section 3: Let us, for example, consider initial conditions for 
which Cj(0) decays as Bj~ a , with 2 < a < 3. In this case, the total mass is 
finite, but none of the higher integer moments converge. G^((, t) has therefore 
two singularities in ( at large times t: one at, and the other near, the origin. 
This last is the one giving rise to the scaling contribution, since we have just 
shown that convergence to a scaling form does take place under these circum- 
stances. On the other hand, the second and higher moments do not exist, due 
to the presence of the former singularity. It is therefore rather misleading to 
say that c(m,t) is "well approximated" by £ _2 e _J// '. A more correct statement 
is that at large times, one has 



lim G c (C,t)| 

t— >oo 



z=s/t 



(4.12) 



4 



(4.13) 




Cj (t) ~ t'^ij/t) + 4r 2 e- J '/( 2 *), 



(4.14) 
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where <&i(x) is a function which can, in principle be computed explicitly. In 
fact, one has, using a similar approach as in the derivation of (3.14): 

c[xMx)] = ^1-^12-^(3-.)^ (415) 

[z + s) 



where C denotes the Laplace transform. It behaves as x 2 ~ a as x — > 0, but 
becomes asymptotically equal to x~ a as x — > oo. Since a > 2, it follows that 
the second term of (4. 14) always dominates when j/t is kept fixed and t — > oo. 
On the other hand, the first term eventually dominates as j ; — > oo for fixed 
t. Choosing an appropriate definition of r' and w', however, one finds that 
t' — t and w' = w, so that strong scaling still holds in this particular case. 

Quite generally speaking, one should not confuse the fact that the distribution 
approaches a scaling limit, which is a global property of the full distribution 
at large times, with good approximations for individual values of Cj(t). The 
exact solution for monodisperse initial conditions could easily lead to wrong 
impressions on that subject: indeed, for (4.10), one has 

Cj (t)=t- 2 exp{j/t) [l + 0(jt- 2 ) + 0(t' 1 )] (4.16) 



for all j and t. This does not generalize to arbitrary initial conditions, however. 
In particular, for fixed values of j the constant 



,(») 



lim 

t^oo 



t 2 Ci 



(t) 



(4.17) 



is in general not equal to 4, and depends on the initial data, even when these 
decay exponentially, as is readily verified in specific instances [108]. If the 
initial condition does decay exponentially, however, then cf°^ also approaches 
4 exponentially fast [65] , as a straightforward estimate of the contour integrals 
yielding Cj(t) from (4.9) will show. 

4.2 The Sum Kernel 

The sum kernel is given by 

K(m,m) =m + m. (4.18) 



It has been solved exactly first by [79,87], see also [23] for further references 
on the subject. I will show two approaches, which can both be used to obtain 
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the full solution: the first involves a transformation, introduced by Lushnikov 
[67], which can always be performed when the kernel is of the form 

K(m, m) = f(m) + f(rri), (4.19) 



where / is an arbitrary function. On the other hand, the usual approach 
using generating functions works quite well also. The reason for showing both 
approaches is, as stated in the Introduction, to give a broad overview of the 
different ways to attack a given problem. 

Let us therefore first consider the discrete case with monodisperse initial con- 
ditions. We may then define, for arbitrary kernels of the form (4.19), 

r(t) 00 

L k =l c k{t) k= i 



where the Cj(t) are defined from the c(m, t) via (2.4). One then sees that (2.2) 
becomes 

-Jf = E fWMj-k - (4-21) 
fc=i 



This set of equations is recursive and can sometimes be solved explicitly. In 
the case at hand, when f(j) = j, the additional substitution 

= <t>je je (4.22) 
does the job, leading to the system of equations 
(Jib- 

^ = E Wrfi-k (4-23) 

lfc(0) = *,■,! (4.24) 
These are solved via the ansatz: 

^(0) = aj e j -\ (4.25) 

The coefficients dj are now determined by the following recurrence 



3-1 

(j - l)a,j = E ka kO-j-k 
k=i 

ai = l. (4.26) 
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This recursion is readily solved, for example using generating functions, see 
Appendix Q, and yields finally 

a, = jf^ y (4.27) 



The sum over all concentrations obeys a closed equation (if we assume that 
mass is conserved for all times, which will turn out to be justified here). This 
leads to 

oo oo 

Y / c k (t) = S e- t S = ]Tc fc (0). (4.28) 

k=l k=l 

Note that now S — 1. Combining all this eventually yields for the concentra- 
tions Cj(t) 



Note here an amusing connection: if one considers (4.21) for f(j) = j and 
make the substitution Xj = <f>j/j, one obtains 



x 



j 



\ k U - k)x k Xj- k - jxj (4.30) 



k=i 



the Xj therefore satisfy the Smoluchowski equations (2.2) for the kernel K(k, I) = 
kl, at least as long as the total mass Mi(t) is conserved, that is, before the 
gelation time. This is particularly useful to evaluate the exponent r, which 
does not depend on the time scale. A generalization of this remark was shown 
in [118]. 

One finds that the appropriate scaling variable is je~ 2t and the scaling function 
is given by 

-3/2-2/2 

* W = ,4.3!) 



To show that scaling occurs for all initial conditions, however, it is easier to 
use the generating function approach without further transformations. The 
equations for c(m, t) can be cast in the form 



dm'm'c(m' , t)c(m — m', t)dm' — 
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-S e~*mc(ra, t) - c{m, t), (4.32) 

where both mass conservation and (4.28) have been used to simplify the equa- 
tions. If we define F((,t) as 

oo 

F(C,t) = J c(m,t)e Cm dm, (4.33) 



one obtains using (4.32): 



8 t F = (F — S e~*)d c F - F (4.34) 
F((,t = 0) = f((). (4.35) 

This can be solved by the method of characteristics, with the following result, 
shown in appendix G 



F((,t) = f(( )e- t (4.36) 
C = Co-[/(Co)-/(0)](l-e- t ). (4.37) 

From this it follows after some algebra, that, if we assume that the initial 
cluster size distribution c(m,t) has finite second moment, then 

lim [d c F(C*)]L =se _ 2< = , L— ■ (4.38) 



1 - 2/"(0)s 



which gives for the scaling function $(rr): 

-3/2 -x/[2M 2 (0)] 

= . (4.39) 

v /27rM 2 (0) 



where we have used the fact that the second derivative of f(() is the sec- 
ond moment of the initial distribution. It is certainly both remarkable and 
somewhat anomalous that the scaling function depends (albeit in a somewhat 
trivial fashion) on the initial condition. 

This dependence on the initial condition, as well as the appearance of the 
second derivative of the generating function in the expression for the solution 
make it of particular importance to study the effect of highly polydisperse 
initial conditions, that is, such initial conditions as have c(m, 0) behaving as 
CmT a as m — > oo, where 2 < a < 3, as we saw in the case of the constant 
kernel. To avoid tedious repetitions, we use the remark made in (4.30) in order 
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to map this problem onto the corresponding problem for the product kernel, 
which will be treated in Appendix L. There we shall find that the typical size 
does not grow as e 2t any more, but rather as exp[(a — l)t/ (a — 2)]. This is the 
first sign of non-universal behaviour in this system. One also finds 



Note how, for a = 3, the traditional value 3/2 is recovered. For a > 3 the 
above computations cease to be valid. On the other hand, for a — > 2, one sees 
that r approaches the limiting (and not allowable) value of two, for which the 
total mass diverges. 



4-3 Parity- dependent Kernels 



We now consider a variation of the constant kernel which is specific to the 
discrete case, namely a kernel in which the reaction rates depend on the parity 
of the sizes of the reactants but are otherwise constant. We shall therefore use 
the discrete notation throughout this subsection. This kernel is defined as 
follows: 



K(k,l) 



K if k and I are odd 
L if k and I are even 
M otherwise 



(4.41) 



This model was introduced in [63,64]. The physical justification given there 
was that such parity effects might arise in the diffusion-limited aggregation of 
alternating linear copolymers. While this seems at best questionable, there is 
no doubt the model provides an important testing ground for the universality 
assumption implicit in the use of the scaling approach. Indeed, such an alter- 
nation, within constant bounds, of reactivity according to parity is a perfect 
example of the sort of detail that ought to be irrelevant if we want to be able 
to apply scaling theory in concrete situations. It does not matter that this 
particular instance of deviation from constancy may not appear reasonable 
physically: there is every reason to believe that similar deviations may pro- 
duce at least similar effects, so that the exact solutions we can find in this 
case are important because they teach us, to what extent the scaling theory 
may be expected to hold, and in what aspects it is likely to fail, once we leave 
the familiar ground of exactly solved systems. The discussion of the solution 
falls therefore squarely within the principal purpose of this review. 
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It turns out to be impossible to give a full and explicit solution for the kernel 
(4.41). What can be done, however, is the following: on the one hand, one 
can prove in full generality that the kernel (4.41) satisfies scaling as defined in 
(3.12), with a scaling function that is, however, different from the one given in 
(4.11) for the constant kernel. It is also possible to give some of the non-scaling 
asymptotics exactly for the general kernel (4.41). On the other hand, it is also 
possible fully to solve some particular instances. I shall restrict myself to the 
case in which L = AM, for which simple and completely explicit formulae are 
obtainable, as shown in [14]. For these, then, everything is known and can be 
discussed in full detail. 

Let us first consider the scaling limit as defined in (3.12). As stated above, 
the scaling limit differs from the one obtained for the constant kernel. This 
invalidates to some extent the usual claims of universality of the scaling the- 
ory. However, it does not contradict what was stated in Appendix A: indeed, 
there, we had shown that if a kernel behaves asymptotically as a homogeneous 
function, then the scaling function should be given by a solution of equation 
(3.29). But the kernel (4.41), which oscillates for ever around the constant ker- 
nel, does not satisfy this hypothesis. Still, this change in the scaling function is 
quite unexpected. The scaling function $(x) for the kernel (4.41) is evaluated 
in detail in Appendix H. As we shall see, however, the set of odd and even 
concentrations obey completely different asymptotics. How, then, can we have 
scaling for the full distribution, without having to divide it into odd and even 
clusters? The answer lies, again, in the very particular nature of the way in 
which the scaling limit is defined, see (3.12). Indeed, what is considered there 
is the average of a function which varies on a macroscopic scale, namely the 
typical size s(t). Such an average is clearly unaffected by any even-odd oscil- 
lations, since they are effectively averaged over. One might at first then argue 
that this is an artefact of an inappropriate definition. However, if one thinks 
about it from the purely physical point of view, the very reverse is the case: 
any reasonable physical measurement, say light scattering, will measure such 
an average and will be quite insensitive to the oscillations. In fact, the only 
measurement likely to feel such effects is a complete measurement of the clus- 
ter size distribution function at the level of accuracy of the monomer size unit. 
This is only very rarely practicable 18 , so that one usually has to be satisfied 
with the coarser measures which do not perceive the oscillations. However, the 
mathematical solution certainly has these oscillations, so that, as we shall see, 
no statement can be made stating that the individual concentrations are well 
approximated by the scaling form. 

In Appendix H it is shown that both the odd and even cluster-size distribu- 
tions approach scaling in the sense of (3.12) separately. This is given precise 
meaning in (H.ll). The scaling functions §d{x) and & p (x) for the odd and even 



A determination of this nature has, however, been achieved in [9]. 
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distributions respectively cannot be evaluated analytically, but the r and d 
exponents can. Here, as throughout in this subsection, the letter d refers to the 
odd (dispari) clusters, whereas the letter p refers to the even (pari) clusters 19 . 
One first defines 



K-M 



Poo = 



1 + 1 1 + 



KL 



(k - My 



(4.42) 



The expressions for r djP and d p are then found to be: 



M 
~K 



Po 



r p = mm 



M L 



(4.43) 

(4.44) 
(4.45) 



whereas for dd we only know that it is positive. Note finally that from the 
definition of scaling given in (H.7) immediately follows that the full cluster 
size distribution tends to a scaling limit with the scaling function <£>(x) given 
by 



<&(*) = \ [<&«,(*) + %(x)} 



(4.46) 



which is quite an extreme illustration of the fact that the scaling function need 
not yield a good approximation to any of the concentrations Cj(t) at all. 

Let us now turn to the non-scaling asymptotics: as usual, define the exponents 

c 2j+1 (t) = const. ■ r w 'd (4.47) 
c 2j (t) = const. -r w 'p (4.48) 

It is then straightforward to check, as is also shortly discussed in Appendix 
H, that 



M , 

^ = 1 + ^(1+6x0 

/ 2M „ 

w' p = min (1 + —poo , 2w d - 1 ) (4.49) 



19 This somewhat singular notation arose first in [14]. It is unambiguous and easy 
to remember after some time, so I decided to stick with it here. 
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Since these coincide with the values of Wd and w p that are obtained from the 
values of (4.43) and (4.44) for r d and t p via the scaling relation (3.25), we may 
state that strong scaling holds. 

The unexpected feature in this case is rather the non-universal nature of the 
result. That is, in general, one does not expect a simple odd-even oscillation 
in the reaction rates to produce such massive changes in the scaling function. 
However, one must remember that the constant kernel is in case II, which is 
quite non-universal in itself. In fact, as is shortly pointed out in Appendix H, 
these peculiarities may be typical of case II. 



4-3.1 Two Fully Solvable Subcases 

We may cite two cases that have been solved more or less fully, which therefore 
allow for a test of the above computations. The first is given by 

L = AM (4.50) 



In this case a completely explicit solution has been found in [14] for the initial 
condition 

c j (0) = a6 jtl + ±(l-a)6 jt2 . (4.51) 



It is then given explicitly by the formulae 



2j\ a [.s(0]~ 
3 , 



(4.53) 



'W = mtW 1 + aKt) + ^OTT (1 + aKt) ' Lim - (454) 

The quantity s(t) defined in (4.54) is therefore the most natural definition of 
typical size in this system, since the forms which define C2j+\{t) and C2j(t) are 
already very close to scaling form in their exact expression. One immediately 
verifies from (4.52,4.53) that = 1/2 and r p = 0, all of which indeed confirms 
exactly the results obtained in the previous subsection. Similarly, one verifies 
the value of zero for The exponent $d on the other hand, which could not 
be evaluated in general, has in this case the value 1/2. It therefore coincides 
with r d in this particularly simple case, once again satisfying the Dominant 
Singularity Hypothesis. It is not clear in general whether this will always be 
the case for these kernels, however. 
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Another case which can to some extent be solved exactly is the one in which 
M = \(K + L) (4.55) 



In this case the kernel is of the form of a sum kernel and can be solved 
using the techniques mentioned in subsection 4.2. The solution obtained in 
[64] is extremely opaque, however, and does not yield such explicit formulae 
for the concentrations as are given for the case L = AM. However, after some 
complicated calculations one obtains values of and r p again in full agreement 
with those found in the previous subsection. 

4-3.2 The addition- aggregation kernel 

In a recent paper [80] a minimal model showing anomalies very similar to 
those found in the parity-dependent kernel was studied. One defines constant 
rates, with the only additional feature that monomers are singled out and 
react differently from other aggregates, more precisely 

K{l,l) = e 

K(l, k) = K(k, 1) = 1 (4.56) 
K{k,l)= 1 , 

where k, I > 2. For this model one obtains a fairly bewildering array of large 
time behaviours for monomers of fixed size, yet scaling in the sense in which 
we have defined it continues to hold, as the authors show in an Appendix. Fur- 
thermore, the scaling function is the same as for the constant kernel, thereby 
vindicating the usual claims of universality. 

4-4 The q-sum kernel 

In a recent paper [13], a model quite different from the classical ones was 
solved, namely 

K(k, I) = 2 - q k - q l = 2 - e - bk - e~ bl . (4.57) 

where q = e~ b is an arbitrary number between zero and one. We shall use q and 
b interchangeably in the following. Again the model is restricted to discrete 
values of the masses only. This model is interesting for at least two different 
reasons: 

(1) It is a non-trivial variant on the constant kernel, in which the reactivity 
of the small aggregates is diminished with respect to that of larger ones. 
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The difference is only significant for rather small aggregates, at least in 
the general case in which q is fixed and not too near to one, so that one 
expects no major deviations from the ordinary scaling picture. As we 
shall see, this is indeed the case if we interpret scaling in the sense of this 
paper, that is, as the convergence in first moment as defined in (3.29). 
On the other hand, for non-scaling asymptotics, quite unexpected results 
are encountered. 

(2) On the other hand, it is also possible to consider a simultaneous asymp- 
totic behaviour in which q — > 1 at the same time as t — > oo. In this case 
we obtain a highly peculiar crossover situation in which the sum kernel 
drives the aggregation of small clusters, whereas the large clusters obey 
a constant kernel kinetics. We shall not have much to say about it, as the 
evaluation of the crossover function is a daunting task which is left to the 
inspired reader. It can be shown, however, that the formalism developed 
in 3.6 indeed yields the correct scaling variables. 

In order to obtain the solution, we note that the kernel (4.57) is of the form 
f(k) + f(l), with f(k) given by 1 — q k . One therefore proceeds along the same 
lines as described in subection 4.2 for the sum kernel, namely we use the 
transformation given by (4.20). This leads to the equations 

rich- 

-§ = EC 1 - <t)Mj-k - (1 - </%■ (4-58) 

These equations can in principle be solved recursively. The solution becomes 
very soon quite intricate, however, and it is not possible to obtain much in- 
formation in this way. 

The solution makes heavy use of the q factorial and q exponential, defined as 
follows 



n-l 



(a;q) n =l[(l-aq l ) (4.59) 

1=0 

oo 

e q {x) = \{{l-xq 1 )- 1 



1=0 

oo 



= V-^-, (4.60) 

where the last identity is familiar from the theory of the q exponential [33]. 
The solution, which is derived in detail in [13] is given as follows: one defines 

oo 

H(C,0) = j:M8){e bK -l) (4-61) 
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This function satisfies a non-linear differential-difference equation, the solu- 
tion of which is shortly sketched in Appendix I. The final result is given by 



#(C,0) = -^Ml + SK,fl)] (4.62) 

From these results and the definition of 9, see (4.20) one readily obtains 

U-Lf ^'j-' (4.64) 
e q {q) ^ {q; q) r 



From these results essentially everything can be computed, though this can 
sometimes be difficult. 

Let us first concentrate on the case in which q is fixed. In this case, the only 
limits of interest are the large-time limit at fixed j and the scaling limit. It is 
quite straightforward to obtain from (4.58) the asymptotic order of magnitude 
of the Cj(t). With some additional work one obtains the more precise result 



It is therefore clear that the large time asymptotics is markedly different from 
that of the constant kernel. In particular, the exponent w' is not even defined, 
since the large-time decay law depends on the size j of the cluster. The only 
feature that remains is the fact that the exponent tends to the value 2 for 
large values of j. On the other hand, an elementary analysis shows that this 
behaviour is only attained when 

t > exp(q~ j ) (4.66) 



which increases extremely fast with j. We have then obviously no contradiction 
with the scaling statement, which asserts that the decay goes as t~ 2 as long as 
j is of order t. Along similar lines, we see that the moments M p (t) for p < q — 1 
are dominated by ci(t) and hence decay as t~ 2+q instead of going as t p_1 as 
predicted by scaling theory. This is the case even though, as long as p > —I, 
the scaling expression for the moments converges at the origin. 

The scaling limit can also be shown to exist. This is proved by showing that 
the function H((, 9) has only one singularity which approaches zero as t — > oo, 
which is a simple pole. Due to this the scaling characteristics of the constant 
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kernel are obtained. The detailed proof of this fact is quite technical and is 
given in [65]. 

Let us now turn to the issues involved in the crossover regime, when 6^0. 
This has been treated fairly extensively in [65] , to which the interested reader 
is referred for more details. As things stand, this appears to be an instance 
in which the formalism developed in subsection 3.6 could be applied if we 
substitute K(k,l) by K(k,l)/b, so as to satisfy the requirement that K(l, 1) 
be always of order one: indeed one has A = 1 with b playing the role of the 
small parameter e in (3.89). One therefore expects a typical size of order e* 
and scaling parameters je~ f and bj. I have not been able to carry out this 
evaluation successfully, however. 

The results which have been obtained are the following: 

(1) When bt is of order one, the sum kernel is a good approximation both for 
clusters of fixed size j and for clusters in the scaling region. On the other 
hand, considerable discrepancies always exist in the large-size limit, no 
matter how short times are. This is related to the fact that the nearest 
singularity of the g-sum kernel is a pole, so that the d exponent is zero. 
Similarly, we may say that moments of fixed order p are, in this range, well 
approximated by those of the sum kernel. However, large order moments 
will show deviations earlier than low order moments. 

(2) For sufficiently large times, all features of the constant kernel are recov- 
ered. This is an immediate consequence of the statements made about 
the existence of the scaling limit. However, the time needed to reach in 
all respects this scaling limit is surprisingly large: the necessary condition 
is 



(3) A more accurate description of the regime lying between these two ex- 
tremes would clearly be highly desirable, but I am not aware of a way to 
obtain such results. It can be shown that the crossover limit in the sense 
of subsection 3.6 exists with the scaling variables je~ 2t and bj. This is 
shown in Appendix J. 



4-5 The Bilinear Kernel 

The most general kernel which has been solved in classical works on the subject 
is the so-called bilinear kernel: 



* > exp(l/6). 



(4.67) 



K(m, m')=A + B(m + rri) + Cmm'. 



(4.68) 
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The history of the various solutions to this kernel is too long to be told in 
detail. The crucial issue concerns the case in which C > 0. As we have already 
stated in subsection 3.5.4, such a case, involving as it does a homogeneity 
degree A > 1, necessarily leads to violation of mass conservation in finite time. 
While it was originally believed that no solutions could be found beyond that 
time (see for example [76,77]), it was eventually realized [115,58] that such 
solutions existed and were physically meaningful. The decrease in total mass 
was interpreted as the flow of mass to an infinite cluster, in a way quite similar 
to that which is observed in systems such as percolation: as the size of a finite 
system is taken to infinity, the total mass comprised in those clusters which are 
of a size that diverges as the system size may or may not be zero. In the latter 
case, this mass does not appear in the final bookkeeping once the system has 
been taken as infinite. There is hence no contradiction involved. Furthermore, 
from the mathematical point of view, it can in fact be shown that if the second 
moment of the cluster size distribution diverges, mass conservation can be 
violated. We shall further show explicitly how this happens in the bilinear 
model (4.68). 

Let us first limit ourselves to the simpler model 

K(m,m) =mm (4.69) 



We shall prove in Appendix M that for purposes of scaling (4.68) and (4.69) 
are equivalent. The possibility of solving the general kernel (4.68) exactly 
does, however, open new vistas concerning crossover, as we show in the next 
subsection. 

In order to gain some insight into the behaviour of (4.69), we introduce the 
following generating function 

oo 

G(C,t) =Jm'c(m', t)e Cm '. (4.70) 
o 



This satisfies the following PDE 

G t — [G — Mi(t)} Gq = 0. (4.71) 



Here we have not assumed that Mi(t) is constant, since this hypothesis cannot 
be sustained over all times, though a solution to the Smoluchowski equations 
(2.2) with the kernel (4.69) does indeed exist for all times. Let us first con- 
sider the case in which M\(t) is equal to one. In this straightforward 
application of the characteristic method yields 
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G[C(*)J=0(Co) 

C(*) = Co-^(Co)-i) 



(4.72) 



where g(() is the generating function of the initial conditions. For each value 
of t, the function ((t) takes its maximum value at a value of Co such that g'(Co) 
equal \jt. Above this value G((,t) is not defined. But G((,t) must clearly be 
defined for all negative values of (. The maximum value must therefore be at 
least equal to zero. But it is easy to verify, using the convexity properties of 
g(() following from the positivity of the initial conditions, that if t > l/g'(0), 
the maximum value of ((t) is strictly negative, so that the function G((,t) 
is not defined for all negative ( for such a value of t. This would imply a 
cluster size distribution c(m, t) that is exponentially increasing with m, an 
utter absurdity. For lesser values of t, however, everything is well-defined and 
the solution has all the expected properties. As an example we may cite the 
instance of monodisperse initial conditions, for which the exact computation 
is an easy excercise in contour integrals performed in Appendix K. The result 
is 

Cj (t) = -J^-t^e'^. (4.73) 



This is indeed a well-defined solution, but it is obvious that for t > 1 the 
corresponding generating function stops being a solution to (4.71) with Mi(t) 
equal to one, since all functions Cj(t) decrease. In fact, these are the so-called 
Flory solutions of the Smoluchowski equations corresponding to (4.69). They 
are obtained by starting from the assumption that Mi(t) is always equal to 
one, and therefore eliminating it from the Smoluchowski equations. In other 
words, instead of solving the original problem (2.2) for the kernel (4.69), one 
solves instead 

oo 

c(m, t) = | J m'{m — m')c(m ,t)c(m — m',t)dm — mc(m,t). (4.74) 
o 



Here Mi(t) has been set equal to one by fiat. It is then straightforward to 
check that (4.73) is the solution of (4.74) for all times. However, one then 
finds for this solution 



(* < 1) 

Mx(t) = { ~ (4.75) 

(t>l), 
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where ((t) is given as the solution to the equation 



Clearly, (4.75) casts some doubt on the validity of (4.74) beyond t — 1. It has, 
however, been defended on the following grounds: beyond the gel point, it is no 
longer possible to identify the finite systems and its infinite limit as we have 
always implicitly done so far. Rather, in the finite system, for times larger 
than the gelation time, we have two cluster populations: on the one hand, we 
have the ones that are small with respect to system size, the concentrations 
of which tend to a well defined limit as it is made to go to infinity. These 
are known as the sol particles. On the other hand, the large clusters of the 
finite system have sizes that go to infinity with the system size, and contain a 
finite portion of the mass. These so-called gel particles do interact, in the finite 
system, with the sol. Since, in the finite case, mass is certainly conserved, we 
might expect some variation of (4.74) to be derivable for the finite system. If 
one then goes to the limit in this fashion, the solution (4.75) is obtained for 
the sol for all times, and (4.76) simply tells us how the mass contained in the 
finite clusters decreases with time. As we shall see later, however, there exist 
exactly solved finite systems which tend either to the Flory solution (4.74) or 
to the so-called Stockmayer solution, which is the post-gel solution of (4.71). 
It is therefore not straightforward to decide which of the two solutions is the 
"correct" one. The issue is, to some extent, simplified by the fact that, in 
any realistic situation, the mean-field assumptions made in the derivation of 
the Smoluchowski equations (2.2) certainly break down near the gel point, so 
that the whole issue becomes irrelevant. A detailed discussion of the difference 
between the Flory and the Stockmayer solutions is found in [115]. 

Let us now derive this latter solution. The central issue is a geometric one, 
discussed in greater detail in Appendix K: it is found that the function G((, t) 
covers the (£, t)-plane twice and has a singularity over some curve on the 
(C,t) plane, on which the two branches of the surface join smoothly. This 
structure is a universal feature of all solutions of (4.71), independently of the 
shape of M\(i). However, since it is manifest that only the sheet of G((,t) 
that connects to ( = — oo without crossing the singularity, has meaning as 
a generating function of a solution of (4.71), it is clear that we must choose 
M 1 (t) in such a way that the singular curve always lie above ( — 0. On the 
other hand, after t — 1, it follows from the previous dicussion that we cannot 
make Mi(t) to remain constant. We must therefore allow a non-trivial time- 
dependence of Mi(t). But it is clear that such a dependence can only obtain 
if the second moment of the cluster size distribution diverges, that is, if the 
singular curve coincides with the ( = axis. This computation is carried out 
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in detail in Appendix K for monodisperse initial conditions, yielding 



Cj{t) = < 



(4.77) 



Note that both the values and the first derivatives of the Cj(t) are continuous, 
so that they really are (strong) solutions of the set of first order ordinary 
differential equations (2.2). 

Finally let us make some remarks on the applicability of the scaling approach. 
As always, one must first define the sense in which we want scaling to hold. 
Since most of the mass is contained in clusters of size of order one, even as 
the gel time is approached, the usual definition of scaling given by (3.12) does 
not hold. As we have already stated before, the definition of second moment 
scaling (see (3.15) for the precise definition) is the most appropriate one. The 
scaling function as t approaches t c is computed in Appendix K, with the result 



$(x) = —^x- 5 / 2 e- x/(2a) 

a=0^l (4 78) 

Two remarks are in order at this stage; First, the ordinary scaling we have just 
derived rests upon the finiteness of the third moment of the initial conditions. 
If this diverges and yet the second moment is finite, we are led to an interesting 
situation: the scaling behaviour is then given by different exponents, but the 
gel time remains finite. After this gel time, however, the asymptotic behaviour 
must necessarily have the exponent —5/2, since that is the only one that allows 
a finite rate of mass transport. In Appendix L we therefore study in detail what 
happens under these circumstances, as well as what happens in the case of a 
divergent second moment, in which case instantaneous gelation occurs. 



4-6 Crossovers in the Bilinear Kernel 

The bilinear kernel (4.68) has three variants which display crossover behaviour. 
Using the type of parametrization defined in (3.88), these are 

K(m,m') = l + e(m + rri) (4.79) 
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K(m, m') = l + e 2 mm (4.80) 
K(m,m') = m + m! + emm' (4-81) 

Let me first recall the definition stated in (3.90) of the crossover limit and the 
corresponding scaling function 

OO , x oo 

lim J m'c(m',y/s(t);t)f - W J x<£>(x,y)dx. (4.82) 



Here, as in subsection 3.6, the function c(m, e; t) denotes the solution of the 
Smoluchowski equations (2.2) with a kernel characterized by a given value of 
e. 

For the first case, the scaling function is evaluated exactly in Appendix N. It 
is given by the following expression 

M X y) _ V(2sy)**exp [-2xy/(l - e^)] 



This expression reduces to the scaling function (4.11) for the constant kernel 
when y — > 0, as of course it should. Here s(£) is taken to be t. The limit of 
large y is somewhat more subtle: taking the limit x — > oo and y — > oo with 
xye~ 2y = z constant, one obtains the correct limit apart from a divergent 
prefactor due to issues of normalization. This illustrates nicely the difficulties 
involved in studying the large y behaviour of such functions. However, the 
use of (4.83) is not, of course, in the large y regime, where everything is 
best described by the scaling function (4.39) of the sum kernel, but rather at 
intermediate times, for which the aggregates have already become large but 
not large with respect to e v . 

On the other hand, for any finite value of y, $(x, y) tends to a constant as 
x — > 0, as stated in subsection 3.6. Figure 1 shows the scaling function (4.83) 
as a function of x for various values of y. Note how the initial values at low 
x are always constant, corresponding to a r = regime, whereas a x~ 3 ^ 2 
behaviour arises over a very extended range for large values of y. 

It is also possible to go part of the way toward evaluating the crossover func- 
tions for one of the two other cases arising in the full bilinear kernel, namely 
(4.80). For the crossover from the constant to the product kernel given by 
(4.80), one finds the following result, the derivation of which is sketched in 
Appendix O: we first define 



<?(s,2/) = 7— - -2 T2 T7i : 7 ( 484 ) 

2 + Sq — [Sq cos y + 2s sin y) 
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Fig. 1. The crossover scaling function is plotted as a function of x for various 
values of y, namely 0.1, 0.5, 1, 2 and 10. Note how a long x~ 3 / 2 regime is eventually 
established. It requires, however, the very largest value of y. 



y 

s = s -t + J dy'q(s,y'). (4.85) 



Here (4.85) must first be inverted and the result substituted into (4.84) to yield 
the desired result. One then has for the crossover scaling function $(x, y) 



q(-s,y) = J x$(x,y)e- sx dx (4.86) 



The analytical (or even numerical) evaluation of §(x, y) from these equations 
appears to be a non-trivial task, which is left to future work. It is easy to 
verify, however, that this function has a singularity at y equal to tc/2, which 
corresponds to the gelling singularity. Indeed, as will be remembered from the 
definition of the crossover limit, the parameter y in this case is nothing else 
than a rescaled time et. But it is easy to verify, using for example the moment 
equations, that the gel time for the reaction rates (4.80) is of order 7r/(2e) as 
e -> 0. 

On the other hand, computing the scaling function in the case (4.81) of a 
crossover betweeen the sum and the product kernel presents considerable dif- 
ficulties. A naive approach does not appear to lead to true scaling behaviour 
in a straightforward way. While this is in itself a matter of some interest, it 
should be left as a subject for future work, as the issue is not settled. 



79 



4-7 Finite Systems 



It has been shown by Lushnikov and Piskunov [70] that the following finite 
versions of the classical reaction kernels Ki, K 2 &ndK 3 as defined in (4.1, 4.2, 
4.3) can be solved exactly at least in the discrete case: 



These have all been solved exactly in [70] . We shall first discuss in some detail 
the case of K[ N \k,l), showing how the exact solution can be obtained and 
discussing its behaviour also from the point of view of the crossover theory 
devfeloped in the previous section. We shortly state the exact solution for 
K2 N \k, I) and go on to analyze the limiting behaviour of K^ N \k, I). This last 
is of particular interest, as it tends to the kernel K 3 (k,l) as N — > oo, which 
displays gelation. On the other hand, it is obvious that K^ N \k,l) cannot 
display such a phenomenon. As we shall see, one can show that K^ N \k,l) 
leads to the Stockmayer kinetics as defined in subsection (4.5), that is, the 
exact solution to the equations defined by the kernel K% N \k,l) tend to the 
solution of the full Smoluchowski equations, not to those for which one has 
performed the substitution of the first moment by the constant value of the 
total mass (Flory solution). On the other hand, as can also easily be shown 
(see [4]) the following finite version of K 3 (k, I) will lead to the Flory solution 



We therefore have the very remarkable phenomenon, that the limiting be- 
haviour of the finite systems depends crucially on the nature of the cutoff. 
The reason for this has been pointed out already by Ziff et al. [117,24], where 
it is shown that the Stockmayer solution corresponds to a situation in which 
large aggregates do not react with small ones, whereas the Flory solution cor- 
repsonds to the case in which they do. Here, by "large", I mean comparable 
to the size of the cutoff. In this case, it is understandable that should 
tend to the Stockmayer solution, since reactions between large particles are 

— (AO 

strongly restricted, whereas in K 3 this is not the case. 

At this stage it is important to clear up a possible misunderstanding: with 
the above kernels one aims to study the gelation transition by viewing it as 
the limit of a set of non-gelling kernels, which represent systems with only a 
finite number of different reacting species. This should be sharply contrasted 




(min(ife, /) < N) 
(mm(k, I) > N) 



(4.87) 




max(fc, I) < N 
max(k, I) > N 



(4.88) 
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with the regularization induced by considering only a finite number of reactive 
aggregates in the first place. The effects of such an approximation are much 
more complex, since one must then take into account the qualitatively new 
phenomenon of fluctuations in the aggregate numbers. This has been done 
by van Dongen and Ernst in [104], which presents the whole issue with great 
clarity. This, however, belongs to those interesting subjects, the treatment of 
which falls outside of this review. 



4-7.1 The Finite Constant Kernel 

We shall here treat the case of the kernel K± N \k, I) as defined by (4.87). The 
exact solution can only be found in a straightforward way for the case of 
monodisperse initial conditions, which is why we have limited ourselves to the 
discrete case. The solution in this case is given by 



Cj (t) = ^- 1 exp 
e 

t — 2 J dy exp 



N nk 



fc=o 



k 



N k 
. fc=0 K 



(4.89) 
(4.90) 



For completeness' sake the derivation of this solution is given in Appendix P. 
It should, by the way, be emphasized that the solution here given is, in a sense, 
incomplete: Indeed, the functions Cj(t) for N < j < 2N are different from zero 
and are not given by the expressions (4.89). These are inactive clusters formed 
by the aggregation of clusters of size less than N. Eventually, of course, all 
the mass will come to be concentrated in those clusters, whereas the mass 
contained in the active clusters goes to zero. 

Qualitatively, one sees that the solution behaves as follows: as long as 9 remains 
significantly less than one, one finds that the solution given by (4.89) is very 
close to the exact solution of the constant kernel given in (4.10). On the other 
hand, as soon as 6 becomes significantly larger than one, the concentration 
profile becomes exponentially increasing and a large amount of the total mass 
finds itself in the inactive clusters. It also follows from (4.90). that the time 
at which this change in behaviour takes place is of the order of the cutoff N . 

It is, of course, quite obvious that such a kernel in the limit of large N offers a 
good example of a crossover situation, as discussed in a similar case in [111]. 
From the exact solution we may obtain the crossover behaviour as follows: 
first define a scaling variable £ given by 



t = N{l-6). 



(4.91) 



81 



We now take note of the following identity valid in the scaling regime in which 
£ is of order one, also shown in Appendix P: 




-£i(0-hi£ 

Ei(|£|)-ln|£| 



(£>0) 
(£<0) 



(4.92) 



where the exponential integrals Ei(x) and Ei(x) are defined as in [1]. If one 
now substitutes these expressions in (4.89) and (4.90) one eventually obtains 



where x is the scaling variable j/t. We see therefore that t 2 Cj(t) depend only 
on the scaling variables x and t/N, or equivalently £. The scaling variables 
required are indeed those predicted by the crossover theory developed in sub- 
section 3.6: the exponent A is equal to zero and e is equal to 1/N, as is readily 
seen if we write as follows 



As to the qualitative behaviour of the crossover function described by (4.93) 
and (4.94), it must be admitted that is is not wholly transparent. However, 
it does describe the transition between an exponentially decreasing cluster 
size distribution and an exponentially growing one as £ changes sign, which it 
clearly must do, in view of (4.94): as soon as t/N is larger than the (convergent) 
integral from £ to oo given on the right-hand side, £ must become negative, so 
that the concentration profile becomes exponentially increasing. Note further 
that the scaling picture shown here implies that for times of order N, that is 
in the scaling regime, the total mass in the active range is still of order one, 
and can be given in terms of the scaling variables by the following expression 




(4.93) 



(4.94) 




(4.95) 



where X[o,i]( x ) is the characteristic function of the unit interval. 



M 1 (t) 



exp[2£!(0] 
exp[-2Ei(|£|)] 



(*>0) 
(£<0) 



(4.96) 
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4-7.2 The Finite Sum Kernel 



For the kernel K < ^ J \k ) l) ) the results are entirely similar. I therefore limit 
myself to stating the exact solution, which is given by 



9 — \ exp 
A 



oo 

— / dw 



A 



Ef =1 q^- (fc+1) 

1 - EfcLi a k w~ k 
dff 



l-Ef = iO fe A(^)- fe 

N \ 

l-]Ta fe A(0)- fe 
k=i ) 



(4.97) 




(k-l)\ 



The derivation of this solution proceeds on exactly the same lines as that 
of the constant kernel and is given in Appendix Q. The detailed analysis of 
the crossover is left as an excercise to the enthusiastic reader. Concerning the 
asymptotic behaviour, however, one can make the following observations: as 
follows from the computations in Appendix Q, for N = oo, EfcLi < l kW~ h has 
a singularity at w — 1, where it has a value of one. There is hence for every 
finite N a value w c (N) < 1 such that 

TV 

J2a k w c (Ny k = l, (4.98) 

k=i 



which is the limiting value to which A tends as 9 — > oo. Since the singular- 
ity in the integral representing 6 is a simple pole at A = w c (N), 9 diverges 
logarithmically with t, in contradistinction to the N = oo case, for which it 
saturates at one. As for the constant kernel, we must distinguish between times 
for which the simple pole at w c (N) dominates and earlier times., for which we 
are far enough from the singularity for the finite sum to be well approximated 
by the infinite sum. 



4.7.3 The Finite Product Kernel 

The exact solution for K^ N \k,l) is surprisingly simple to obtain using the 
same approach as in the previous two cases. The final result is given by: 



83 



- A 

t = A exp J 
.0 



i-Eti^ 



(4.99) 




Using these formulae, one sees that A grows monotonically with time starting 
from zero. If N is infinite then, upon reaching the critical value e -1 , A will 
yield a finite value of the time, but will not be able to grow beyond this value, 
since the integrand then has an integrable singularity. This critical time is 
indeed the gel time and is, in this case, equal to one. Beyond this time, the 
above equations do not quite state unambiguously what happens, but a closer 
look at the derivation shows that A remains at its critical value of e -1 . We 
have therefore rederived the Stockmayer solution in this way. 

If, on the other hand, we take N to be finite, the behaviour becomes qualita- 
tively quite different: indeed, the nearest singularity of the integrand on the 
positive real axis is now a simple pole, since the denominator is a polynomial. 
Since this polynomial approximates the corresponding analytic function quite 
well, we see that this singularity lies beyond the value e _1 by an amount of 
order 1/JV. THerefore, as time goes on, A will approach ever nearer to this 
pole without ever reaching it. It then follows that the finite aggregates (those 
having j <C N) are given for all times very nearly by the Stockmayer solution. 
This therefore shows our claim, that the kernel K^k,l) yields in the limit 
N — > oo the Stockmayer kinetics for the finite aggregates. Note, on the other 
hand, that for times larger than t c , the aggregates of size of order n will be 
significantly different from the value predicted by the infinite model, since A 
is distant by a quantity of order 1/N from its theoretical value, thus implying 
a change for those aggregates by a factor of order one. 

Let us now look at the kernel K^\k,l). Since it is bounded from above by 
const. ■ N(i + j) it cannot display gelation, so that, if one considers all aggre- 
gates, mass must be conserved. From this follows that the kinetic equations 
(2.2) for K^\k, I) for the aggregates of size 1 < j < N read 



which are manifestly the same as the Flory equations (4.74) as pointed out by 
Bak and Heilman in [4]. They therefore have the same solution for the active 
clusters, so that the solutions of (4.100) indeed straightforwardly tend to the 
Flory solution. 

We are hence in the highly peculiar situation that the limiting dynamics de- 
pends on the sequence by which the infinite system is approached. It is intu- 



j'-i 

\ k U - k)c k Cj- k - jc 



(4.100) 



fc=i 
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itively clear that this can only happen in a gelling system. The issue at hand 
is that of the interaction between finite aggregates (having a size j <C N) and 
"gel aggregates", of size of order N. In fact, it follows from a careful study of 
the global existence results proved in [58] that such a situation can never arise 
in the situation studied there, namely whenever the kernel is bounded from 
above as follows: 

K(k,l)<r k n (r k = o{k)). (4.101) 

and, in fact, more generally whenever v < 1 with v defined as in (3.5). The 
peculiar phenomena described above can therefore only arise for a kernel with 
v = 1 and for which gelation occurs. (It has been shown by van Dongen [107] 
that kernels with v > 1 have instantaneous gelation, and in fact their solutions 
may in fact perhaps not exist at all). The intuitive reason for this is as follows: 
when v < 1, the reaction rates between a j'-mer and an iV-mer, with N ^> j 
is of order N v . If we now average this reaction rate over all gel aggregates 
and use the fact that the gel mass is finite, we see that the average reaction 
rate tends to zero. It is therefore essential that v — 1 for these pathologies to 
occur. This may perhaps diminish the sense of urgency in finding a "solution" 
to these problems. The above results do show, however, that no resolution can 
be expected from the consideration of finite systems when v — 1. 

4-8 Pure Three-body Reactions 

It is fairly straightforward to generalize the methods described in the preceding 
subsections to higher order reactions. Let us first describe what has already 
been done in the literature. In [42,43] Jiang et al. described the behaviour of 
the following n-particle kernel: 

n 

K( mi ,...m n ) = H(A mi + B) (4.102) 
i=i 

These show gelation with exponents z = — 2 and r = 2.5. Furthermore, the 
case of monodisperse initial conditions can be solved exactly in this case as 
well. These results are obtained using an extension of the techniques described 
in [58]. This agrees quite well with the overall picture given in Section 3 (see 
in particular subsection 3.7), that as far as scaling properties are concerned, 
many-body reaactions do not involve major changes from binary, apart from 
replacing A by A — (n — 2). It can also straightforwardly be verified that the 
constant kernel reaches scaling from arbitrary initial conditions with r = 
and z — l/(n — 1). In these examples we therefore do not see any variation 
from the expected picture. 
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We now turn to the two other examples given in (3.134) and (3.135). These, 
as we stated before, show some significant peculiarities. K 1 , in which the rates 
are given by m\ + m 2 + m 3 , can be straightforwardly solved by an extension 
of the techniques described for the sum kernel, see [50]. As usual, let us define 
the generating function: 

oo 

F((,t)= Jc(m,t)e mi , (4.103) 
o 



one finds the following PDE for F((,t) 



d t F = \[F 2 - S(t) 2 ]F c - S(t)F (4.104) 
S(t) = F(C = 0;t) (4.105) 
/(C) = F(C;t = 0) (4.106) 

These equations can again be solved via the characteristic function technique, 
as discussed in appendix R, with the following result 



C- 1 + t/(o) "^" (0)c °- (4108) 

From this, the solution for monodisperse initial conditions can be obtained 
explicitly and is given by 



, * (2k + If -1 1 / t \ k ( 2k + 1 t 
C2fc+i(*) = T^T, h-TT ex P 



2 k k\ l + t\l + tJ V 2 1 + t J 
C2*(t) = 0. (4.109) 

One finds that the typical size s(t) grows as t 2 . This leads immediately to the 
following scaling function: 

x -3/2 e -x/4 

= 2 ^ ■ (4-110) 



The remarkable result, however, is not in the scaling function, but rather in 
the rate of growth of the typical size. Ordinarily, as we have seen in many 
occasions, this quantity is not difficult to find out and its value rests on very 
simple-minded, and hence very robust, scaling arguments. Here, however, this 
is not the case: a perfectly straightforward scaling theory predicts that z — 1 
for this kernel, but one finds through an exact solution that z — 2. It is 
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also easy to see, using the same kind of techniques that were developed for 
other systems, that this scaling function is in fact valid for arbitrary initial 
conditions. We see, therefore, that this is an example of the phenomenon 
discussed earlier, namely in subsection 3.7: one indeed readily verifies that 
the condition r — /i < 1, which is necessary for the validity of the scaling 
equation (3.129) for three-body reactions is here violated. Since the derivation 
of (3.129) fails because of a divergence of the integrals involved, we see that 
the expression for the typical size will also be given incorrectly: in fact, the 
real typical size is larger than the one predicted by (3.125). The following 
rough argument may help to understand what is happening: if we cut the 
integrals appearing in (3.129) at a lower limit of l/s(t), these become time 
dependent quantities, diverging as s(t) 1 / 2 , thus leading to a modified version 
of the equation (3.125) for s(t), namely 

s(t) oc s^)*- 1 ^) 1 / 2 oc s(t) l '\ (4.111) 
which yields the correct large-time behaviour for s{t). 

It is also straightforward to verify that, for strongly polydisperse initial condi- 
tions, that is, if c(m, 0) goes as m~ a with 2 < a < 3, then both z and r vary 
continuously with a (see Appendix R). This is reminiscent of similar proper- 
ties for the (two-body) sum kernel. Nevertheless, the peculiar behaviour of z is 
much more striking in this three-body case, since we are dealing with continu- 
ously variable exponents, not a continuously variable exponential growth rate. 
Of course, since the various exponents are undetermined by our general con- 
siderations, there is no a priori reason why they should not vary continuously 
with the nature of the initial condition. 

Let us now consider instead the kernel K 2 given by 

K 2 (k,l,m) = kl + km + lm. (4.112) 

In order to avoid tedious repetitions, let us use a simpler though less complete 
method of solution: we write down the equations for the various moments M p 
as follows 

M v = \ £ '-^M P1+1 M P2+1 M P3 , (4.113) 

Pl+P2+P3=P 

where the primes indicates that the sum runs only over those indices for which 
no two of the pi simultaneously vanish. From this follows in particular 

M = -|M 
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M 2 = MfM + 2M 2 . 



(4.114) 



These equations can be solved exactly. However, it is almost immediately 
apparent, without any further work, that M 2 diverges at finite time. This is 
confirmed by the exact solution, which reads, for the case of monodisperse 
initial conditions: 



We therefore find that this kernel gels, although the straightforward formula 
for z leads us to expect an exponential growth of the typical size. From this 
follows, incidentally, that a theorem such as that proved in [112] cannot be 
generalized in the expected way to the pure three-body case. 

Summarizing, we have given two examples for which the naive scaling argu- 
ments give a value for the growth exponent z different from 1/ (2 — A3) , which 
is the result of the straightforward scaling approach. We may now ask: What 
happens in the more realistic case, in which the three-body term perturbs a 
given two-body aggregation process. In the next subsection, we shall consider 
the cases in which the kernel K\ perturbs the constant kernel, and the kernel 
K 2 perturbs the sum kernel. Thus in both cases we consider the case in which 
the two-body reaction rate have the same growth exponent z as 1/(2 — A 3 ) 
for the three-body kernel. In other words, we ask whether the naive value for 
the growth exponent may not play some role when the three-body dynamics 
is taken together with a corresponding two-body aggregation. As we shall see, 
the answer is (in part) affirmative. 



4-9 Three-body Reactions Perturbing Two-body Aggregation 

Let us now consider the case in which both two and three-body reactions 
occur. We therefore consider the system of equations 



M 2 (t) 



2-e 



(4.115) 



oo 




dmidm 2 K^ '(mi, m2)c{mi, t)c(m 2 , i) x 







[8{m\ + m 2 — m) — 5(mi — m) — 5(m 2 — m)} + 



(4.116) 



00 



+— / dmidm2dmsK ( '(mi, m 2 , m3)c(mi, t)c(m 2 , t))c(m 3 , t) x 






[5(mi + m 2 + m 3 — m) — 5(mi — m) — 5(m 2 — m) — 5(m 3 — m)} 
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Here we shall discuss the two following cases 



K( 2 \mi,m 2 ) = 1 K^\mi,m 2 ,m 3 ) = mi + m 2 + to 3 

(4.117) 

K^(mi, m 2 ) =mi + to 2 K®(m 1) m 2 , m 3 ) = m 1 m 2 + mim 2 + to 2 to 3 



In Appendix S we show in some detail using the method of characteristics that 
the first case does indeed satisfy scaling with the following exponents 

(a < \) 

z={ Aa (4.118) 

(« > 2) 




l + 2a 



This can also be seen more easily using the method of moments. Indeed, one 
finds 



Mo = -(| + «)M 2 

M 2 = l + a + 2aM M 2 . (4.119) 

From this the values of z stated in (4.118) follow straightforwardly if we assume 
scaling instead of proving it. This shows therefore that when the sum kernel 
perturbs a two-body kernel, it behaves very much as if it had the value of z 
predicted by a naive scaling theory: that is, it does not affect the growth expo- 
nent of a two-body kernel which itself has z equal to one unless the strength of 
the perturbation reaches a certain threshold value. Above this threshold value, 
on the other hand, it does affect the growth exponent, but only to the extent 
of making it continuously variable, rather than imposing its own value. This 
last only happens in the limit in which the three-body perturbation becomes 
infinitely strong. Yet another peculiar feature of this model is the following: 
it is an instance in which adding to a given aggregation mechanism, namely 
the three-body kernel, a very strong second mechanism, namely the constant 
kernel, actually leads to a slowing down of the first mechanism. This should 
be compared with a conjecture framed in [60] in which it was hypothesized 
that if 

K 1 (m,m')> K 2 (to, to') (4. 120) 



for all to, to', then the typical size of the system corresponding to K 1 would 
grow faster than the one corresponding to K 2 . Counterexamples were given in 
[12], but they involved channeling particles towards non-reactive states. The 
above example, although not strictly to the point since it involves ternary 
reactions, certainly shows the need for great care in making such claims. 
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From the results of Appendix S we also obtain the scaling function $(x), albeit 
in a very implicit form. However, this allows to estimate its small x behaviour 
and hence the exponent r. One obtains the result 



{a < 1/2) 
(a > 1/2) 



(4.121) 



Note that from this and the scaling relation (3.25) follows that, for a < 1/2, w 
is equal to one. However, as is easy to see, for monomers (and more generally 
speaking, whenever j is held fixed and t goes to infinity) the decay exponent 
w' is given by 



w = 



2 + 2a 
T+2a' 



(4.122) 



We therefore see that in this model strong scaling is violated. It is arguable 
whether this is a more "natural" model than the q model as an example of this 
phenomenon. In the q model, we require a non-homogeneous kernel, whereas 
here we have a homogeneous two-body kernel perturbed by a three body kernel 
of a degree of homogeneity which, according to a naive scaling argument, 
should yield the same growth rate. The type of violation involved is also, in a 
sense, less drastic, since the w' exponent is at least defined independently of 

3- 

On the other hand, writing down the moment equations for the sum kernel 
perturbed by the three-body sum-of-products kernels yields a gelling system: 
this is easily seen by examining the equations for M 2 (t) and M (t), which 
only involve one another and can be solved elementarily. Again, it is easy to 
verify that the solution for the second moment diverges no matter how small 
a is. We see that the behaviour is sometimes the one that would have been 
expected on naive grounds, but that one cannot count on it. 



4-10 Constant Kernel with Production Term 



The constant kernel with a time-independent production term can also be 
solved exactly. The results are quite instructive, so we present them here in 
some detail. The equations are 



d t c(m, t) — \ j c{m 1 )c{m 2 ) [5{m - m 1 - m 2 ) - 
o 

S(m — mi) — 5{m — m 2 )]dmidm 2 + P5(m — m ). (4.123) 
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For simplicity, we assume that no aggregates are present at time zero. From 
this follows that only aggregates with mass an integer multiple of m are 
produced, which allows to use the notation for discrete initial conditions. 

Under these circumstances a straightforward application of the generating 
function technique yields the following solution: define 

oo 

G(C,*) = E C *(*) (e*-l) 
fc=i 

One then obtains after a straightforward calculation 

G = -y / 2P(l -eC)tanh (\t^2P{\ - e<)) . (4.125) 



(4.124) 



Note how this function is meromorphic with a countable set of simple poles 
on the positive real axis which crowd toward the origin as t — > oo. Since all 
singularities are simple poles, one finds that $ = 0, whereas, as we have seen 
before, r = 3/2, thereby violating the Dominant Singularity Hypothesis. 

We now use the identity 

OO J 

tanhx = X | E o l2 + (2; + 1)%V4 . (4.126) 



in order to find the solution 

J2 oo 



Cfc(t) 



7T 



{21+ 1) 2 7T 2 ' 

2Pt 2 



Pt 3 



i + 



i=— oo 



-(fc+1) 



(4.127) 



which yields the various results stated in the previous section in a rather 
straightforward way: indeed, if we let t — > oo at fixed values of k, one obtains 
a constant value, since (4.127) is then a Riemann sum which can be replaced 
by an integral, namely 



lim Ck (t) = 1 X ^—rdx = \ — T [ k 7 ^ . (4.128) 

*-oo feW i i (l + x 2 ) fc+1 V vr + v ; 



This confirms the general results that such models tend to a stationary value, 
as well as the specific value for the limiting concentrations. Indeed these are 
well-known to be given by the second expression in (4.128). 

However, from (4.127) one also finds the general scaling form for large sizes 
and large times. This goes beyond the well-known results for the stationary 
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distribution, giving a fairly complete overview of the dynamical behaviour. 
The result is readily obtained by noting that, when k is of the order of t 2 , only 
the first few terms of the infinite sum in (4.127) contribute, and that these can 
be approximated accordingly. However, a simpler approach is the following: 
remembering the expression (3.14) for the scaling function <&(#), we obtain 
the following expression for the Laplace transform of $(x): 

oo 

Jdx^(x) (e- px - l) = -v^tanh^l (4.129) 
o 

from which it is straightforward to obtain the scaling function. 



4-11 The constant kernel with a localized source term 



Let us now consider the case in which one has the constant kernel with a 
localized source term and a diffusion term, that is 



oo 

9,c(m,f; t) = I j K ( mi , mMmM^Mm - - m 2 ) - 



— 5(m — mi) — 5{m — m'i)\dm\dm2 + 
+P5(m - m )5(r) + DAc. (4.130) 

We essentially follow the treatment given in [17]. In this case, it is an obvious 
observation that the moments obey a closed set of equations, namely 



p-i 



M p (f,t) = \ I P I M q M p _ q + Pm p 5(r) + DAM p . (4.131) 



These are diffusion equations with a recursively given inhomogeneity, so they 
can clearly be solved iteratively. Again, however, the explicit expressions rapidly 
become so unwieldy that they do not yield any understanding of the system. 
It is nevertheless possible to gain some insight by looking for scaling solutions 
of (4.131) as follows: define 

M p (t) = t a ^ p . (4.132) 

Matching powers on both sides of (4.131) and noting that 

J M 1 (f,t)df= Pm Q t, (4.133) 
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which sets the value of a± to 1 — d/2, one finally obtains 
a p = (2 - d/2)p - 1 



(4.134) 



The fact that the a p increase as p is crucial for the approximations to be 
consistent: indeed, the production term in (4.131) contributes to the total 
mass M 1 (f, t), and is neglected for higher values of p. This is only valid for 
d < 4, so that the upper crittical dimension of this system is 4. Indeed, 
for d > 4, it can be verified by a similar power counting argument that the 
aggregation terms can be discarded and the whole problem reduces to diffusion 
with a constant source. On the other hand, as we shall see in Section 5, for 
dimensions one and two the hypotheses of mean-field theory are bound to fail 
at some stage. This does not mean that the theory developed here is never 
applicable: if the reaction probability upon contact is very small, there will 
be a large intermediate regime for which the theory just developed will apply. 
However, eventually, it will fail, and the issues discussed in Section 5 become 
of central importance. 

A considerable simplification occurs in (4.130) if the time-dependence of the 
l.h.s. is neglected or, in other words, if we consider the stationary approxima- 
tion. This is in fact what was done in [17]: it is shown there that the c(m,r) 
have a scaling form defined as follows: 



The reader is referred to the original paper for further details. In particular, 
an exact solution can be obtained for the one-dimensional mean-field 
well as for the diffusion-limited case, for which techniques analogous to those 
to be discussed in the next section are necessary. 



5 Beyond Mean-field 

As has been stated many times, the mean-field equations (2.2) are only a 
valid description if no spatial correlations between the aggregates arise. The 
reason for this requirement is clear enough. Implicitly, we have always had 
the following picture in mind: aggregates, in order to grow, must first be 




(4.135) 



where the exponents r and z are given by 



z = A- d 
d-6 



(4.136) 



T d-A 
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brought into contact by some transport mechanism, and then react. Both the 
efficiency of the former and that of the latter may depend in an arbitrary man- 
ner on the masses m and m! of the two reacting aggregates. The reaction rate 
matrix K(m, m') is therefore a product of a transport-related collision cross- 
section k(m,m') and a reaction efficiency upon contact given by some matrix 
a(m,m f ). In this picture, however, we take it for granted that those pairs of 
aggregates which the transport mechanism is likely to bring together repre- 
sent an unbiased sample of all aggregates, in other words, that no correlations 
between aggregates exist on the scale determining their transport. 

This is, as we shall see in the following, a non-trivial assumption. More ex- 
actly, we shall see that in the majority of physically relevant cases, it is likely 
that this hypothesis cannot be sustained for large times. On the other hand, 
in very many systems, there is a large regime of intermediate times for which 
the mean-field theory described in the first part of this paper is valid. This 
is due to the fact that, for many systems of practical relevance, the limiting 
factor in the reaction is not the transport mechanism itself, but rather the 
efficiency of reaction upon contact. In that case, it is straightforward to see 
that, over the times necessary to produce one single aggregation event, the 
transport mechanism will have randomized the environment of the reacting 
pair so efficiently that the mean-field assumption is unproblematic. If, on the 
other hand, the transport mechanism is the rate-limiting factor, then matters 
are much more difficult: as we shall see, the transport and the reaction mecha- 
nisms may then collaborate to create correlations between particles, leading to 
a situation in which mean-field theory fails. In the following, we shall therefore 
not be concerned with solutions of (2.2) any more. Rather, we shall consider 
specific models for irreversible aggregation, in which a given transport mecha- 
nism is specified, as well as a reaction mechanism. We then proceed to analyze 
the cluster size ditributions generated by these models, comparing them with 
possible solutions of (2.2) and using the scaling techniques developed in the 
rest of this paper to characterize the behaviour of these distributions. 



5.1 Diffusion-limited Cluster- Cluster Aggregation 

The simplest model for irreversible aggregation, inspired from the physics of 
colloids, for which the transport mechanism is Brownian diffusion, is known 
as the Particle Coalescence Model (PCM) and goes as follows: consider point 
particles diffusing on a lattice and carrying an integer which we call their 
mass. When two particles, of mass m and m' respectively, alight on the same 
site, they combine at a given rate k to form a particle of mass m + m'. The 
rate k is often taken to be infinite, that is, the formation of the new particle 
is instantaneous. The new particle is again pointlike. On the other hand, the 
rate at which it diffuses may depend on the mass m as an arbitrary function 
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D{m), which goes asymptotically asm 7 . 

Note in passing that this model already contains a considerable amount of 
structure. Its main unrealistic fetaure is that the particles remain pointlike. 
The great advantage of this assumption is, that it frees us from considering 
issues of morphology. Indeed, if we assume that clusters connect rigidly when 
they approach within one lattice spacing of each other, the model becomes 
much more complex: in particular, it then generates fractal structures, and 
the radius of the cluster will depend in a highly non-trivial way on the mass. 
Such a model, simulated in three dimensions, already deserves to be called a 
realistic model of colloid formation in the Brownian regime. However, since 
nearly nothing can be said about it at the analytical level, I shall not treat it 
any further, but refer the interested reader to the extensive literature (see e.g. 
[47,78] and references therein). 

The PCM, introduced by [48] can be treated rigorously as long as the diffusion 
constant is mass-independent. In this case, see [66] it has been shown that for 
space dimension d > 3 the model behaves qualitatively like the mean-field 
model with a constant kernel. While this may seem encouraging, we shall in 
fact see that this feature is intimately linked to the pointlike nature of the 
particles. On the other hand, in two space dimensions, logarithmic corrections 
appear, whereas in one dimension an explicit exact solution is available [93,94], 
which displays a number of highly anomalous features. 

For arbitrary diffusion constants, on the other hand, no such rigorous treat- 
ment is available. Let us shortly argue what mean-field theory we should 
expect to hold: as is well known, the current which flows into an immobile 
absorbing sphere of radius R surrounded by a sea of particles diffusing with 
a diffusion constant D is DR d ~ 2 for d > 2. From this we may plausibly argue 
that the rate at which one diffusing particle encounters another is propor- 
tional to D\ + D2 as well as to + R2) d ~ 2 ■ Here Ri are the radii of the two 
particles and Di their diffusion constants. Thus, if Ri is always the same, we 
may disregard it and we are eventually left with 

K(m, rri) = D(m) + D(mf). (5.1) 

If the space dimension d < 2, matters are more complex, since any two par- 
ticles are sure to react eventually. As we shall see, however, for d < 2, mean 
field theory is not expected to work at all, so the precise form of the reaction 
coefficients is a rather academic issue. 

At this point, it may be worthwhile to point out that it is occasionally possible 
to mimic the results produced by correlations through an appropriate choice of 
reaction term, different from quadratic. Here I shall systematically disregard 
such possibilities: when a reaction involves two bodies, I shall take the attitude 
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that only a quadratic expression can be viewed as an adequate reaction term 
describing it. 



5.1.1 Simple Scaling Arguments 

In the following, I shall show how very straightforward arguments of the kind 
that have been used throughout this paper, can be used to predict the be- 
haviour of the PCM. In the following subsection, I shall briefly present the 
exact one-dimensional solution due to Spouge [93,94], as well as shortly men- 
tion the corresponding rigorous results concerning the two-dimensional model, 
thus partly confirming the results to be derived now. 

First let us derive the growth of the average cluster size. Let us assume that a 
"typical" particle starts out by being a monomer and grows in such a way as 
to be typical, that is, of size s(t) all the time. If we consider the probability 
P(x, i) of finding the particle at x at time t given that it was at the origin at 
time t — 0, one has 



— = D[s(t)]AP = s(t)^AP 

P(x,0) = 5(x). (5.2) 

From this follows that the probability of finding the particle at a distance of 
order L from the origin is of order one if 

L = O (y/D [s(t)\ tj = 0(t (1 -* 7)/2 ). (5.3) 



Therefore we expect the particle to have covered a space of order L in time t. 
Since originally the total mass in this interval was of order L, and it did not 
change except through diffusion at the boundaries, which can be neglected, 
we expect the typical size to grow as L, or in other words 



Now let us turn to the values of r and w. At this (low) level of sophistica- 
tion, we shall assume that strong scaling holds throughout and shall have no 
qualms in identifying the decay exponent for monomers with w. We may pic- 
ture the monomer decay problem as follows: each monomer is surrounded by 
two particles of typical size. These grow normally by coalescing with other 
particles outside the interval which connects them and which contains the 
monomer. These particles' growth in no way affects the monomer's survival. 
The problem therefore reduces to a three particle problem, where the two 
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outer particles have a time-dependent diffusion constant, whereas the middle 
particle diffuses normally. The question is then to determine the decay law for 
the survival probability of the middle particle. This can then be identified with 
the monomer concentration decay at large times. This problem is somewhat 
technical and was studied in [35], with the following result: if 7 < 0, that is if 
the aggregate diffusion constant grows with size, there is a fair probability that 
the two large particles will diffuse quite far away, letting the middle particle 
survive for quite long times. The result then is given by 

ci(t) = canst. ■ r 2 ' {2+ ^ (t -> 00) (5.5) 



For 7 = 0, on the other hand, a straightforward solution is known [48,27] 
which yields 

Cl (t) = const. ■ r 3/2 (t -> 00). (5.6) 



If we now apply (3.25) to the values of w obtained from (5.5) and (5.6), one 
obtains 



T 





-1 



( 7 <0) 
(7 = 0) 



(5.7) 



If we now compare this with the strict mean-field prediction for which we use 
the kernel (5.1), which is really the only one with a more or less fundamental 
significance for the PCM, we obtain for the mean- field growth exponent zmf 

1 . . 

zmf = — — • (5.8) 
1+7 



With respect to r MF our knowledge is less complete, since kernel (5.1) is of 
type II, for which no closed expression for r is available. Nevertheless we know 
from (3.47) that it satisfies the inequality 

tmf < 1 - 7- (5-9) 



Both (5.8) and (5.9) are, of course, incompatible with our scaling theory. 

At this stage, one might nonetheless wonder whether some appropriate modi- 
fication of the equations could possibly describe the system. From the value of 
z one deduces that the degree of homogeneity A should be equal to —7/2. The 
fact that w 7^ 1 excludes a type I kernel, so we are left with the possibility of 
type II. From the inequality (3.47) we see that 7 > r, which is indeed always 
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satisfied. But in fact, there is a far more serious objection to such attempts: 
let us compute the probability of finding a particularly large cluster. This is 
clearly related to the possibility that any given particle might go particularly 
far in time t. This probability decays as a Gaussian, and not as an exponential. 
But we have seen (see subsection 3.5.2), that under very general hypotheses, 
the large size behaviour of a cluster size distribution is always given by an ex- 
ponential. We see therefore that even in this apparently exceptionally simple 
case, no simple adaptation of mean-field will give an adequate description of 
the system's behaviour at large times. 

Let us now turn back to the case in which 7 > 0. Then the two particles 
surrounding the monomer diffuse apart at an ever decreasing rate, leading 
thereby to very rapid trapping of the central monomer. More specifically, one 
has 

ci(t) = const. ■ exp —const. ■ r (t — > 00) (5.10) 



The determination of f3 involves severe analytical difficulties. Hellen et al. [35] 
suggest the following approximate value obtained from a Lifshitz tail argument 

P=P^- (5-11) 

' 4 + 27 y ' 



Again this is in qualitative agreement with the Smoluchowski picture, but in 
clear quantitative disagreement: indeed, for 7 > 0, kernel (5.1) is of type III, 
so that we expect small monomers to decay as stated in (5.10), with (3 given 
by 

= 7, (5-12) 



in clear contrast to (5.11). The authors in [35] also state unambiguously that 
even fairly coarse numerical measurements invalidate the mean-field prediction 
(5.12) quite clearly. 



5.1.2 An Exact Solution 



In the one-dimensional case, if the diffusion constants of the aggregates are 
size-independent, that is, if 7 = in the notation of the previous subsection 
and if, additionally, reaction between aggregates is instantaneous upon con- 
tact, an exact solution was found by Spouge in [93,94]. Here I wish to give an 
outline of the way in which this solution can be obtained, as well as to point 
out that it does, in fact, confirm the results suggested by the previous, ad- 
mittedly coarse, scaling analysis. We shall, however, follow a slightly different 
approach in the spirit of that pioneered by Doering and ben-Avraham [21,22]. 
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It is clear, by the definition of the model described above, that we are dealing 
with a Markov process. This implies that one can write down a master equation 
for the joint probability distribution function at time t for finding the first 
particle at x 1 , the second at x 2 , and so on to the iV-th, where the Xk represent 
the lattice position at time t of the particle k. This equation, however, is far 
too complex to be of any use, and also contains more information than is 
really needed. The standard approach consists, as is well-known, in taking 
correlation functions involving one or more particles. In this case, this leads 
to an intractable set of equations, in which each correlation function couples 
to one of higher order. If, on the other hand, we consider a different kind of 
reduced description, which takes the particular one-dimensional nature of the 
system fully into account, one can indeed find a closed set of equations. This 
goes as follows: define E(k,x) as the probability of finding a total mass k 
contained in an interval of length x. The crucial simplifying feature about this 
particular choice of correlation function is, that it remains unchanged by all 
reactions that take place either inside or outside the relevant interval. In other 
words, it is only affected by events in which particles move inside or outside 
the given interval via diffusion. And this, understandably, can be solved. 

Let us translate the above verbal arguments into equations. First define F(k, x) 
to be the probability of having a total mass k inside a given interval of length 
x, and additionally that the leftmost site of this interval be occupied, irrespec- 
tively of the state of the rightmost site. In this case, one has 



where a is the lattice spacing, which also indicates the length which particles 
jump in one diffusion step. 

For simplicity we now restrict ourselves to a continuum limit, in which the 
particles are far apart on the scale a. In this case, we may discard the possibility 
of both ends of the interval being occupied, so we have for the time variation 
oiE(k,x) 



where D denotes the macroscopic diffusion constant arising from the contin- 
uum limit. In the derivation of (5.14) we have made essential use of the fact 
that sites cannot be doubly occupied, which follows from the instantaneous 
reaction rate. We now need initial and boundary conditions for (5.14). One 
finds, after some careful considerations involving the original discrete model, 
that 



F(k, x) = E(k, x) - E(k, x - a), 



(5.13) 



d t E(k, x) = F(k, x + a) - F(k, x) 
= DAE(k,x), 



(5.14) 
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E(k, 0) = 
d x E(k,x)\ x=0 = c k {t) 



(5.15) 



for all times t. Furthermore, if the initial particle density is p , one has 

E(k, x;t = 0) = ^^e- pox , (5.16) 



at least, if the monomers are initially distributed at random 20 . It is now 
straightforward to solve for (5.14) with initial conditions (5.16) and the ab- 
sorbing boundary condition at the origin (5.15). Using the method of images 
yields 



E (k,x) = P ° / y k e'P°y e -{-yr/(^t) _ e -(x +y y/(4Dt) (51?) 



2k\V^Di 



(x-yf/(4Dt) _ -{x+yf/{ADt) 







and hence 



-2 




Analytically it is not quite straightforward to evaluate this expression, or even 
to show that it scales. In terms of generating functions, however, this is easy 
and is shown in Appendix T. The resulting scaling function <&(x) is given by 

$(x) = (5.19) 
v 7r 



and the typical size indeed grows as VDt. We have therefore confirmed that 
the decay at large x of the scaling function is faster than exponential, thereby 
precluding any identification with a mean-field model of the type we have been 
discussing in the body of this paper. We also show rigorously that the r and 
w values are in fact given by the values predicted in the previous subsection, 
namely —1 and 3/2. Finally, we may note that (5.18) shows that strong scaling 
holds in this system: indeed, if t — > oo in the r.h.s. of (5.18) while k remains 
fixed, the integral becomes time-independent and one finds 



lim 



(Dt) 3/2 c k (t)} = P °~ {k R t 1] (5.20) 



20 All the above restrictions can, and have, been lifted. Again, for greater details see 
[93,94]. 
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What have we learned from this example? Primarily we see that in one di- 
mension the mean-field approach fails, essentially because of the buildup of 
correlations between particles due to the interplay between reaction and dif- 
fusion. In particular, one finds that the probability of finding two unreacted 
aggregates nearby is considerably smaller than expected on the assumption 
that the reactants are uncorrelated: thus, if we ask for the probability of 
finding exactly one monomer in each half of an interval of length x, that is 
E(2,x) — 2E(2,x/2), one finds it to decrease as x 3 as x — > 0, whereas the 
square of the probability of finding a monomer in such an interval goes as 
x 2 . There is, therefore, a strong repulsion between nearby particles, due to 
the efficiency with which diffusion in one dimension brings nearby particles to 
react. 

It should be pointed out that we have assumed instantaneous reaction rate 
whenever particles meet. Without this assumption, the model cannot be solved 
exactly as far as is known. However, it is easy to realize that in the limit where 
the reaction rate k becomes very small, the transport mechanism becomes 
irrelevant and mean-field theory holds. However, if particles have an initial 
separation of order L, they will on average meet L times in a time of order 
L 2 , that is, in a time sufficient for them to come close at all. This follows 
from well-knnown properties of one-dimensional diffusion. Therefore, if kL is 
of order one, the problem reduces to the one in which reaction is essentially 
certain. Since the typical interparticle distance L(t) is growing with time, we 
should say that when L(t c )k is of order one, the results obtained above should 
start to hold, whereas before the description using the constant kernel should 
be valid. Numerical work does indeed confirm this expectation. The time t c 
can be evaluated under the assumption that for t = t c . L(t) can be evaluated 
both using the formulae of the constant kernel as those of the diffusive model. 
One then has 



which implies that t c goes as D/k 2 . One sees therefore that it may be possible 
to apply the mean-field theory quite meaningfully, even in a system for which 
we know that it cannot be applied in the true asymptotic regime. 

A similar reasoning leads one to conjecture that the mean-field description 
for the PCM with 7 = becomes accurate in dimensions 3 and higher. This 
is due to the fact that random walks are transient in this case, that is, they 
visit each point only a finite number of times. One does not expect, therefore, 
that the correlations created by the diffusion and the reaction process should 
build up indefinitely in time. Indeed, these expectations have been confirmed: 




(5.21) 
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specifically, it has been shown (see e.g. [66]) that the cluster size distribution 
can in this case be scaled onto an exponential distribution, just as for the 
constant kernel, which is the one we expect to model this particular system. 

Finally, I believe this system also shows very nicely the extent to which the 
concepts originally developed in the context of the mean-field equations can 
be carried over to more general situations: thus the basic scaling relationship 
(3.25) connecting r, w and z follows from mass conservation alone and holds 
in all cases. Similarly, the distinction between power-law polydisperse and 
bell shaped scaling functions also generalizes, as we have seen in the previous 
subsection. 

5.2 Ballistic Aggregation 

Another simple model of aggregation is the following: consider particles of a 
given radius a originally spread at random in 1Z d with random velocities. These 
then move freely, and stick irreversibly upon any two coming closer to each 
other than the capture radius a. In order to specify the model completely, we 
require two additional pieces of information: first, how does the radius grow 
upon irreversible sticking? Here many solutions are in principle possible. We 
will focus entirely on the two following: either the radius does not grow, or 
else it grows deterministically in such a way as to maintain the total aggregate 
volume constant. The second is of course more realistic, but considerably more 
difficult to discuss, as well as more controversial. Second, we need to know how 
the velocity of the compound particle is obtained. Since we are physically as- 
suming that this represents sticking of freely moving particles, we shall always 
assume that momentum conservation determines the final velocity, which it 
does in a unique way. Since the two particles combine to one inelastically, we 
have, of course, no conservation of energy. 

Let us first consider the one-dimensional case. The two growth rules for the 
radius are then equivalent. This model was first considered by Carnevale et al. 
[15] and later investigated further by Jiang and Leyvraz [45] using numerical 
simulations together with qualitative scaling arguments. We present these first, 
and later discuss shortly an elegant exact solution [75,29], which confirms the 
results anticipated in [45] in a very satisfactory way. 

5.2.1 Simple Scaling Arguments 

As we have seen in the case of the diffusive one-dimensional model, it is often 
a good strategy to search for the long time behaviour of the typical size and 
the long-time decay of very small clusters, since this contains essentially most 
of the relevant information. In [15] it was shown that the typical size in one- 
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dimension grows as t 2 / 3 using the following argument: an aggregate of mass m 
has arisen, in the discrete picture, out of m monomers. These all had momenta 
of order one, which averaged, as vectors, to zero. If we therefore assume that 
they are essentially independent, we are led to assume that the momentum 
pirn) of the m-mer is of order m 1 / 2 , which leads to a velocity v(m) of the 
order of mT 1 ! 2 . This leads to a typical aggregate of size s(t) crossing in a time 
t a length of order t>[s(t)]t, which is of order t/^Js(t). Since the number of 
particles initially contained in this interval will be of the same order as s(t), 
we wind up with 

s(t) = const. ■ ip t) 2/3 , (5.22) 



where p again denotes the initial particle concentration. 

This reasoning is quite robust, but does make an assumption concerning the 
initial velocity distribution function (VDF), namely that it has no large ve- 
locity power-law tails. Indeed, the argument which invokes the central limit 
theorem fails, and s(t) grows in a different way when such tails are present. 
Indeed, one finds [45] that if the VDF has a distribution with a behaviour of 
the type 

p( v ) = const. ■ v'"- 1 (1 < a < 2) (5.23) 



then one has using an obvious extension of the above reasoning, that the 
exponent z is given by 

2 - 5^y. (5.24) 



The assumption a > 1 is necessary, since otherwise the typical velocities would 
grow on average, making the whole approach questionable. No meaningful 
numerical results have been obtained for this case, whereas (5.24) has been 
well confirmed numerically for the range 1 < a < 2 [45]. 

To complete our scaling argument, we should find the asymptotic decay for 
very small clusters. A completely general lower bound on Ci(t) can be derived, 
quite independently of the initial VDF, in the following way: consider the 
aggregates at time t. It is possible, of course, to identify the positions at time 
zero of the monomers conforming them. These fall in consecutive intervals of 
typical length s(t) separated by intervals of length one. If both the following 
events occur, the initial particle will survive until time t, 

(1) The particle should first have been placed on one of the empty regions 
between two of the intervals constituting the two neighbouring typical 
clusters. This has probability l/s(t). 
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(2) The particle should have been launched with a velocity small enough that 
it will not reach either of the two neighbouring clusters in time t. Since 
the two clusters are eventually separated by a distance s(t), this limiting 
velocity is of the order s(t)/t. If the VDF has a finite probability density 
for velocity zero, this leads to a probability of the order s(t)/t as well. 

Combining the two probabilities one obtains 

d(t) > const. • t" 1 . (5.25) 



It further appears quite unlikely that monomers could survive under any other 
circumstances than the ones stated above, so that (5.25) is presumably the 
exact order of magnitude of the monomer decay. Note in passing that we have 
assumed that the VDF is finite at zero velocities. What happens, for example, 
if all initial velocities are taken with values ±1? At first sight, this changes a 
lot, since odd clusters always have a lower (mass-dependent) bound on their 
velocities, whereas even ones always have a finite probability of having veloc- 
ity exactly zero. Therefore, for a fixed even value of the mass one has C2j(t) 
going asymptotically as l/s(t), which differs from (5.25). Similarly, for fixed 
odd values of the mass, we get exponential decay of C2j+i(t) for large times. 
However, it is easy to see that in the scaling limit the continuum approxima- 
tions can in fact be made, so that the scaling exponent w is 1 in this case 
also. 



From this evaluation of w follows via (3.25) that 
l/a (1 < a < 2) 



T = < 



(a > 2), 



(5.26) 



where the usual case is contained in the second. Numerically, these claims have 
been verified 21 in [45] . 

Now let us compare with mean-field theory. To this end, we need to develop a 
plausible expression for the reaction rates K(m, ml). A common and frequently 
used one for related problems (see e,g, [23]), is the following: 

K(m, m') = \v(m) - v(m') | [R(m) + Rim')} 11 ' 1 , (5.27) 



Persistent findings stating that r = in the case of velocities with no power- 
law tail have repeatedly appeared in the literature. Since they are now found to 
be contrary to exact results [75,29], it may be enough to say here that they are 
erroneous. 
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where v(m) is the velocity of a typical aggregate of mass to and R(m) its 
radius. Using the standard approach outlined in this paper, one finds for the 
growth exponent z in the case where v(m) goes as mr 1 ! 2 and R{m) constant, 
that z is in fact 2/3 as expected. On the other hand, in this case, the kernel 
(5.18) is of type III, so that the monomer concentration ci(t) should go as a 
stretched exponential, specifically 

ci(t) = const. ■ exp (-const. ■ t 1/3 ) (5.28) 

It is, of course, quite evident what is going wrong here: in (5.27) we are as- 
suming that monomers always react as if they were moving at a velocity of 
order one, whereas it is clear that the velocity of the surviving monomers is 
a decreasing function of time, due to fast monomers being selectively elimi- 
nated by the reaction. In fact, it is seen numerically that the velocity profile 
v(m) does not decay with mass, but is on the contrary rather flat, decaying 
as s(t) 1 / 2 with time. There is also (somewhat inconclusive) evidence [45] that 
(5.28) actually holds for very small reaction rates, so one may in fact argue 
that one is using the "correct" mean-field theory. 

A more sophisticated approach involves considering the joint distribution of 
mass and momentum. In this case we may view the vector (to, p) as a vector 
valued mass, in a way similar to the one in which we treated multicomponent 
aggregation. One then finds for the reaction kernel 

[R(m) - R(m')] d - 1 . (5.29) 

Since the p are on average zero, we are in the same situation as when discussing 
the scaling of composition. However, it is easy to see that, at least within the 
treatment of composition scaling that we have given, the behaviour of small 
clusters is essentially the same. The gaussian behaviour of the momenta is 
justified in detail by our approach. On the other hand, it should be admitted 
that we have not shown that these are the only possible solutions to the 
multicomponent scaling equations (3.140). 

5.2.2 Exact Results 

Here I shortly review the exact solution derived in [75,28,29,30]. The precise 
form of the solution is quite formidable and its derivation is an impressive feat 
indeed. I shall not, however, go into the details, which are rather technical, 
and for which the interested reader is referred to the original literature. 

In order to state the result in a readily understandable form, let me define 
some notation. I denote by T(xi.yi, X2, y-i) the following probability: Let y(x) 



K(m,p] to' ,p) 



TO 
p 



TO 

7 
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be a random function of x chosen according to the Wiener measure, or, in other 
words, a Brownian motion directed along the x axis. T is then the probability 
that it satisfy 



vM = yi 
y(x 2 ) = y 2 
y(x) < x 2 



(xi < x < x 2 ) 



(5.30) 



In other words, it is the probability that a directed random walk start at 
(xi, yi) and reach (x 2 , y 2 ) without having ever touched the "absorbing parabola" 
y = x 2 . This can also be descibed as the solution to the following partial dif- 
ferential equation 



dT(x u yi,x,y) = d 2 T(x u yi, x, y) 
dx dy 2 



(5.31) 



with the boundary conditions 

T(x 1 ,y 1 ;x 1 ,y) = 5(y-y 1 ) 
T(x 1 ,y 1 ;x,x 2 ) = 



(5.32) 
(5.33) 



From this one now defines a few auxiliary functions: first I introduce Zx>(^i, yi), 
which is the probability of reaching infinite values of x without having been 
absorbed by the parabola. This is given by 



x 

%c(xi,yi) = Jim J T(x 1 ,y 1 ]x,y')dy'. 



(5.34) 



Finally, we shall be interested in the relative probabilities of escape of particles 
that start very near to the parabola. To this end we define 



I(X!,X 2 ) = 



d 2 



dyidy. 



■T(x 1 ,y 1 ;x 2 ,y 2 ) 



yi=xf;y2=x^ 



d 



J(xi) = -^—%o(xi,yi] 
oyi 



(5.35) 
(5.36) 



yi=xf 



The final result can now be stated as follows: the concentrations c(m,v,t) 
obey the following scaling law 



c(m,v,t)=t 1 $(a;, r/). 



(5.37) 
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Here x and 77 represent the following scaling variables 
m 

X= ^t) 2 / 3 

ri=(2at) 1/3 v, (5.38) 
where a is the r.m.s. velocity in the initial VDF. One then finds for rf) 




(5.39) 



Several points should be noted about this result: first, the mass and velocity 
dependence do not factorize, nor do they show the simple kind of interdepen- 
dence we have found in multicomponent aggregation. In fact, it is explicitly 
shown in [28] that the collision probability between two aggregates does not 
factorize. In this way, it is clear that any attempt to build a mean-field the- 
ory is ill-founded, and indeed such attempts have failed to reproduce known 
results. 

The derivation of this exact solution is very interesting in its basic ideas, but 
quite laborious in the details. In Appendix U, I therefore attempt to present 
the concepts to the extent that a sufficiently enthusiastic reader should be 
able to fill in the missing algebra. 

One can then use (5.31) to evaluate T, and hence all other functions defined 
above, in terms of Airy functions. This opens the way to an asymptotic analysis 
of $(rr, rf), and of its integral over 77. This is a laborious undertaking, which 
was performed in the above references, where they confirm the scaling results 
stated in [45]. Additionally, an unexpected large x asymptotic is obtained, 
namely 

00 

J ®(x,ri)dr] = const. ■ exp (-x 3 /12^j (x > 1) (5.40) 

—00 

However, it is not clear whether this might not be an artefact of initial condi- 
tions in which all particles are put on a lattice: it seems incredible that normal 
density fluctuations when particles are put at random on the line should still 
satisfy (5.40). In any case, however, this large x behaviour once more shows 
how impossible it is to treat this model in mean-field terms. 

5.2.3 Higher Dimensions 

The case of higher dimensions is both far more complex and, to some ex- 
tent, controversial. First, it is obvious that it is important, in this case, to 
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specify how the particle radii grow upon aggregation. We shall consider two 
possibilities: 

(1) The particle radius remains fixed when the particles stick together. This 
is obviously unrealistic, but it has the merit of greater simplicity. In par- 
ticular, there is a real hope that mean-field theory might hold in d > 2. 

(2) The particle radius grows so as to maintain the total d dimensional aggre- 
gate volume constant. Since the total occupied volume fraction remains 
always finite, this is a really difficult model. Thus, even three-particle col- 
lisions cannot be discarded out of hand as being asymptotically irrelevant. 
I shall have little to say on that count. 

Consider first the case of fixed radius. We may then attempt to formulate 
once more the mean-field theory. In this case expression (5.27) for the reaction 
kernel still holds, and leads to 



independently of d. The monomer decay predicted by (5.28), namely a stretched 
exponential with exponent 1/3, still holds independently of d. This is in fact 
quite natural, since the reaction kernel defined by (5.27) only depends on d via 
the geometric cross-section, which we have assumed constant. In this sense, 
kernel (5.27) is much more general than the previous analysis indicated. 

Is this mean-field theory likely to be correct? Since the total number concen- 
tration decays as t~ 2 / 3 , the mean free path grows as t 2 ^ -1 )/ 3 , which is always 
more than the interparticle distance £ 2 /( 3(i ) when d > 1. This is therefore ex- 
tremely encouraging, since it shows that typically colliding particles will be 
distant on the scale of nearest neighbour separation. In related work it has 
been shown that the Boltzmann equation becomes exact at large times [85] 
for the model of d dimensional ballistic annihilation when d > 1. A numerical 
confirmation of all this would be of considerable interest, but has not, to the 
best of my knowledge, been performed to date. 

The model in which the radii grow as m 1 ^, on the other hand, is far more 
complex. The obvious mean-field theory has again a reaction kernel of the 
form (5.27), which has values of A and \i given by 



_d-2 
X -^2d~ 

A* 5 ( 5 - 42 ) 

from which one deduces in the ordinary way that 
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2d 



(5.43) 




ci(t) = const. ■ exp -const. ■ t d ^ d+i) 



(5.44) 



There is therefore again a stretched exponential decay of the monomers ac- 
cording to this theory. This is again due to the fact that we do not take into 
account the selective elimination of high velocity monomers in this model. But 
since, in this model, the mean free path of a monomer always remains of the 
order of the interparticle distance (due to the conservation of the occupied 
volume fraction), it appears very likely that such an elimination will in fact 
take place. 

Trizac and Hansen have performed some high quality numerical simulation 
on the two dimensional version of this model [98]. These lead apparently to 
predictions quite at variance with the above theory. However, it should be 
borne in mind that the data were never explicitly analyzed in terms of the 
theory described here, but rather compared with a rather peculiar "mean- 
field" theory [84], which includes, however, some correlation effects by means 
of various uncontrolled approximations. This theory's principal merit was to 
provide a rationale for the observation of a "pure exponential" cluster size 
distribution, together with r = 0, which was then thought to be the correct 
result in one dimension [15]. The agreement of the two-dimensional simulations 
performed in [98] with the theory presented here is therefore rather difficult 
to assess. 

Nevertheless, let us us shortly discuss the most salient differences between 
the simulation [98] and theory. The measurement of the z exponent 22 yields 
a result which is claimed to be significantly different from the theoretical 
value of one. The figure published does not, however, show a very convincing 
straight line behaviour over any range of sizes. The best that can be said is 
that strong crossover effects are certainly present, which make it difficult to 
say with certainty what the actual value of z is. The quoted value of 0.8 (for 
the dilute system, which is the one of greater theoretical interest, as well as 
the one for which deviations from theory are largest) is really quite far from 
the expected value, however, so I feel strongly that a more careful analysis of 
the data is required before any definite statement concerning the validity of 
mean-field theory can be made. 

No measurements either of the monomer decay or of the r exponent were made. 
It is therefore not possible to decide the crucial issue whether the stretched 
exponential decay predicted by (5.44) indeed takes place. The approximately 
exponential appearance of the cluster size distribution does, however, make it 
highly unlikely. 



called £ in [98] 
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On the other hand, the average kinetic energy per particle was measured and 
yielded a power law decay t~ s incompatible with the prediction of mean-field 
theory. Indeed, from the mean-field arguments given above one obtains 

5 = z, (5.45) 

whereas the quoted value for the dilute case is 1.12. This is certainly a serious 
disagreement with the measured value of z. If we take this result seriously, it 
suggests that the cancellation between velocities of colliding particles leads to 
velocities far less than what would be implied by the Central Limit Theorem. 
The intriguing suggestion that this might be due to the kinematics of collision 
was made in [99]. A considerable amount of work clearly remains to be done 
to settle these important issues. 



6 Conclusions and Outlook 

I have attempted to give a comprehensive overview of the way in which scaling 
concepts can be used in the context of the kinetics of irerversible aggregation. 
We have seen that applications are numerous and the method is quite powerful 
in predicting, among others, various qualitative behaviours both for the shape 
of the cluster size distribution and for the time dependence of the rescaling 
factors involved, such as the typical size or length. There remain nevertheless 
several shortcomings to the method, which might perhaps be overcome and 
lead to a technique that could be more quantitative even for fairly complicated 
systems. Here are some aspects which I believe to be of interest. 

• Once one has an integral equation for the scaling function $(x), difficulties 
are by no means over. It is then necessary, in practical applications, to solve 
it, generally by some numerical approach. This turns out to be a difficult 
problem. The various forms into which the basic equation (3.12) can be cast, 
due to the presence of an arbitrary function f(x), might possibly be of use. 
On the other hand, for numerical work, the original form of the equation as 
given in (3.35) is presumably optimal, since it contains the desired function 
in the r.h.s. of the equation, whereas the corresponding relations involving 
exponentials yield objects such as the Laplace transform of 3>(x), which is 
definitely less direct. 

The issue of computing scaling functions becomes particularly acute in 
cases such as crossover or inhomogeneous systems, for which the function 
<&(x, y) has two variables and the integral equation is complicated. In these 
cases, even the behaviour of the function at the origin is not straightforward. 
In this respect, it would be of interest to have a theory of comparable 
simplicity as that of regular variation for functions of two or more variables. 
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It is not clear that this is an attainable goal, however, as functions in two 
variables can exhibit a strikingly complex set of singular behaviours. 
I believe crossover may be of considerable use in practical applications, 
since a competition of two reaction mechanisms is a frequent occurence. 
The progress reported here, while possibly of interest, is certainly only a 
first step: more exactly solved models exhibiting crossover should be looked 
into and those for which I have not been able to give solutions, such as the 
q sum kernel in the limit q — > 1, should be viewed as relevant problems. 
Further, the integral equation has almost yielded no results of interest on 
the form of the two-variable scaling function To obtain a theory 

yielding the form of the singularity at least in some limiting cases is of 
crucial importance. It must be admitted, however, that the existing exact 
solutions do not hold out much hope for broad ranges of simple behaviour. 
Exact solutions for models with production and diffusion are, to the best 
of my knowledge, absent except for some work relying on the simplifying 
features arising from stationarity. I believe any such model would be of 
interest, since we have shown that the stationary limit in this case does not 
arise in an easy way from the non- stationary case. 

Are there cases in which the scaling function $(:r) can be obtained exactly, 
though the kernel itself be not solvable? This would increase considerably 
our understanding of aggregation, in particular if we could obtain rigorous 
results for some type of kernel of type III or type I, for which, so far, no 
analytical results are known. 

Can one go any further than I have done towards elucidating the singularity 
of the scaling function for the gelling case? It should be possible to 
obtain some kind of estimates on r, perhaps by using methods similar to 
those developed by Cueille and Sire for Case II kernels [18,19]. From the 
mathematical viewpoint this is certainly important. On the other hand, 
from the point of view of applications, it may well be enough to compute 
such exponents to sufficient accuracy. To this end, however, it would again be 
essential to have a good way to solve the equation, in spite of the appearance 
of the unknown exponent r in it. 

There also remains a certain mystery around the Case I non-gelling kernels. 
These are quite important, yet there remains some doubt concerning the 
very existence of a scaling function for this case. The numerical work of 
Lee [56] certainly makes a very strong case for existence, but numerical 
evidence, of course, is never entirely conclusive on its own. This work seems 
to point to a highly anomalous patterns of corrections to the leading power- 
law behaviour, that are oscillatory and might, in point of fact, violate the 
hypothesis of regular behaviour near the origin. A different mathematical 
framework should be designed in order to derive such results as are strongly 
suggested by the numerics. 

Finally, I believe corrections to scaling may be of interest: by this I do not 
mean the corrections to the leading power-law of the scaling function: this 
has already been done, at least for Case II in [108]. There is, however, an- 
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other issue, hinted at in [7], namely that scaling behaviour is not reached 
immediately. Thus, initial conditions may create broad transients, which 
should also be treated at least to some approximation. In [7] the conse- 
quences of having an initial condition containing both monomer and N-mer 
are worked out explicitly. More work along this direction would be desirable, 
as this is not at all explored until now. 

There also remain several mathematical problems related to the theory of the 
Smoluchowski equations. First and foremost, I would mention the necessity of 
obtaining optimal results for the system of equations describing coupled dif- 
fusion and aggregation. I believe present results, which only concern diffusion 
constants bounded from below for large masses, are probably not as strong 
as they could be made. Finally, several old problems are still there: proving 
the scaling hypothesis under fairly general hypotheses is one. An issue which 
may also be related concerns the relation between the stochastic models of 
aggregation and the deterministic Smoluchowski equations. In the framework 
of van Kampen's Q expansion, very remarkable results have in fact been ob- 
tained. A corresponding rigorous result was proved in [81], which may well 
pave the way to future progress. 
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A Derivation of the Scaling Equation for $(x) 

The fact that the c(m, t) approach a scaling form, as stated in (3.12), can be 
reformulated as follows: 




(A.l) 



We now rewrite the integrand using (2.2) as follows 




dmi dm,2 mim2c(mi,t)c(m2, t)K(mi, rri2)G [mi/ s(t),m 2 /s(t)} 
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S ^ f dmm^-f[m/s(t)]. (A.2) 



s(t) J ' s(t) 
Here G(x,y) is given by 

G(x, y) = — [(x + y)f(x + y) - xf(x) - yf(y)] (A.3) 
xy 

If we now assume that we can replace i^(m 1 ,m 2 ) by its continuous limit- 
ing form s{t) x K[mi/ s{t), m 2 /s(t)] — something we shall justify shortly — then 
(A.2) is of the same form as the l.h.s of (3.12), since G(x, y) is bounded when- 
ever f(x) is bounded and differentiate. It therefore follows that (A.2) can be 
recast as 



W 2 s(t) x ~ 1 

1 dx dyK{x, y)<5>{x)<5>{y) [{x + y)f(x + y) - xf(x) - yf(y)} 



o 

oo 



-W^\ [ dxx 2 f(x)$(x). (A.4) 
s {t){ 

If this, integrated over t, should give zero, as stated in (A.l), then surely one 
must have 

s(t) = W lS x , (A.5) 



for some constant W±. This allows, in principle, to determine s(t) completely, 
including the prefactor. We find 

s(t) = \(l - X)W 1 t + s(0) 1 - A l 1/(1_A) (A.6) 



Note that some difficulties arise when A = 1, and that the whole scheme be- 
comes inconsistent when A > 1. The former difficulties are somewhat technical 
and will be treated in another part of this paper. The second involve, as we 
shall see, convergence to a scaling form of the second moment rather than the 
first, and therefore require different considerations. 

The determination of W\ rests on issues which strongly depend on conventions 
and personal preferences, such as the exact way in which typical size is defined. 
We therefore assume in the following that W\ has been somehow determined, 
and proceed to set the constant W, which we had previously left free, equal to 
W\. One can then take W 2 out of the whole expression in (A.4). One further 
deduces, from (A.l), that (A.4) must vanish, yielding (3.29) after observing 
that 
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G(x, y) = - [f(x + y)- f(x)\ + - [f(x + y) - f(y)\ 

y *^ 

= G 1 (x,y) + G 1 (y,x). (A.7) 

Since the rest of the integral is symmetric, we may limit ourselves, up to a 
factor of two, to integrating over Gi(x,y). 

The one point that is left concerns the justification of substituting K (mi,m 2 ) 
by 

s(t)- x K[m 1 /s(t),m 2 /s(t)] (A.8) 

Let us first consider the case in which f(x) vanishes identically on some interval 
of the form [0, e]. In this case, the statement is obvious: as s(t) — > oo, the 
masses summed over in A. 2) become increasingly large, and the convergence 
of K to K as stated in (3.30) is enough to justify the substitution. 

We have therefore shown, that (3.29) must hold for all functions f(x) which 
vanish on an interval near the origin. The general statement is an immediate 
consequence of this together with the monotone convergence theorem. 



B Derivation of the Scaling Equation in the Three-Body Case 

From the point of view of the algebraic manipulations involved, the three- 
body case is wholly similar to the corresponding two-body computation, and 
the reader is referred to appendix A for notation and the basic computations. 
The issue we shall be addressing here concerns not so much the nature of the 
manipulations involved, but their admissibility. 

Up to (A.2,A.3), everything is strictly analogous to the case of appendix A. 
G is now a function of three variables given by 

G( X , y, z) = (* + V + ')n* + V + ')-*f(*)-vf(v)-'M. (B.1) 

The crucial difference between the two and three body cases now apppears in 
the following fact: G(x,y,z) is not a bounded function, no matter how f(x) 
is chosen, unless it be a constant. Indeed, choose x and y so that 

(x + y)f(x + y)?xf(x)+yf(y). (B.2) 
It then immediately follows that G(x, y, z) — > oo as z — > 0. The only possibility 
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left is, therefore, to assume that things will go well. If we then go through the 
computations that correspond to those given in appendix A, we obtain 

s(t) = W^ 3 ' 1 , (B.3) 



which is the result stated in (3.125). We may then, as in Appendix A, set W\ = 
W 2 and therefore obtain the result stated in (3.125). There is, however, no 
guarantee that this derivation goes through. It is, in particular, quite possible 
that the integral 

J dx dy dz xyz K(x, y, z) G(x, y, z) ${x)$(y)$(z) (B.4) 



diverges near the origin. This then indicates, from a physical viewpoint, that 
processes in which one very small particle coalesces with two particles of typ- 
ical (but not equal) sizes affect the scaling of the aggregation process in a sin- 
gular way. I shall not attempt to make a scaling theory of such cases, merely 
limiting myself to stating when the results obtained by the above theory are 
self-consistent, that is, when the above mentioned divergences near the origin 
do not occur. It should be clear that all the results discussed in this appendix 
carry over straightforwardly to the general case of a reaction of order r > 3. 

We shall see in subsection 4.8 that various exactly solved cases fall under this 
category. Further, we shall see that in this case, it is not always possible to 
neglect a two-body contribution which shoud be negligible on naive grounds. 
Thus, the kernel mi + m 2 + m 3 has a naive growth exponent z — 1. However, 
the real growth exponent in the pure three-body case is z = 2. Nevertheless, 
if we perturb the constant kernel, which also has z — 1, using this three-body 
kernel, the combined process has a growth exponent of z = 1, that is, the 
expected dominance of the model with the higher value of z does not take 
place. For larger values of the perturbation, on the other hand, the exponents 
vary continuously. 



C The Large-p Behaviour of I(p) 



In the following, we shall always assume that <&(x) has regular behaviour near 
the origin in the sense of (3.41) and that the corresponding power law is given 
by r. We now divide I{p) in two parts: 



C/p 

I 1= I dx / dyl((x,y)§(x)§(y)xe~ px 



1-e 



-py 
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oo oo 



h = J dx J dy K(x,y)&(x)&(y)xe 



-px 



1-e 



-py 



(C.l) 



C/p 

I now treat l\ and l 2 separately. For Ii, a change of variables yields 

/ , n 2 00 C 



h = p -(A+3) $ 



.9, 



dxx T e x dy K(x,y)y 



1 - e~ y . 



(C.2) 











Here the substitution of $(y/p)/<E>(l/p) by y~ T is justified by the fact that the 
range of integration of y is finite. A similar justification for x exists always, 
that is, both for I\ and I 2 . Should this yield divergent integrals, one would 
have to conclude that the regularity assumption was not self-consistent. 



For I 2 , we only rescale x by p and obtain using (3.2) and (3.3) 



1 



I 2 = - / dxx l ~ T e 



00 / I \ 

7 A » Xk U) ^ 



(C.3) 



Assuming further that fc(z) behaves regularly (in the sense of (3.41) again) 
with exponent \x (see (3.4)) one finds the final expression 



( \ / 1 \ 00 00 
-J k I - J J dxx l - T e- x J dyy x -^(y). 
^' C/p 



(C.4) 



Now let us consider the issue of convergence of these various integrals. It is 
obvious by the above calculation that the sum Ii + I 2 cannot depend on the 
constant C. One sees, moreover, that the integral appearing in the expression 
for Ji converges exactly when the one appearing in I 2 diverges. The limit 
separating these cases is 



l + A-//-r = 0. 



(C.5) 



Whichever integral diverges is therefore compensated by the correction to 
scaling generated by the finite cutoff in the converging integral. The expression 
containing the convergent integral must therefore always dominate. We may 
therefore write 

/ = B 1P -^® + B 2P -^ Q k Q , (C.6) 



where Bi and B 2 are given by 
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1 - e~ v . 



Bi = J dxx 1 T e x j dy K(x,y)y~ 
o o 

oo oo oo 

B 2 = J dxx l - T e~ x J dyy x ~^(y) = T(2-r) j dyy x -^(y). (C.7) 



(C.6) must be here understood in the sense that only the term for which the 
integral in the corresponding prefactor converges should be taken into account. 

Finally, let us consider the case in which $(#) decays faster than any power 
near the origin. In this case, one can neglect the contribution due to e~ py in 
the expression for I{p) given in (3.34). It can therefore be rewritten as 



oo oo 

j dxx 1+x <S>(x)e- px J dyk (£) $(y). 



(C.8) 



But one has y ^> x always, since x must be of order \ j p and y cannot be small, 
due to the behaviour of $(y) at the origin. One can thus use the following 
approximation for k(z) 



k(z) = Tz u [l + o(l)} (z^oo) 



(C.9) 



from which the result stated in the text readily follows. 



D The Large-p Behaviour for the Three-body Reaction 



One follows exactly the same strategy as in appendix C to evaluate the various 
integrals involved. The result is, with exactly the same notations, and the 
various homogeneity exponents defined as in (3.127): 



oo C 

B 1 = j ' dx J dydzx x - T e~ x y' T z~ T K{x,y,z) [l - e~ {y+z) ~ 



7 (°r \ ( 7 

B 9 = 2 I dxx x ~ T + M e~ x j dyy- T k x (yjx) / dz z x *- X2 ${z) 



o 

oo 



\0 



\C/P 



= J dxx x - T e~ x J dydzz^-T^y/zWyWz) 

C/p 



(D.1) 
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E The Distribution of Compositions in Multicomponent Aggrega- 
tion 



In this appendix we show how the compositions scale in the multicomponent 
case. To this end we require first a scaling equation of the same type as (3.29) 
or (3.140) for the function \I>(<7, 5). Again, we use the same approach as in 
Appendix A and obtain from (3.137) 



oo oo 

da 



-oo 
oo 



J d55[af a (a,5) + a5f s }^(a,5) = 

-oo 

o oo 

J da 1 da 2 J dS 1 d5 2 K(a u a 2 )^(a u S 1 )^(a 2 , S 2 ) x 

-oo 

x [(Si + 5 2 )f(a 1 + a 2 , 5i + 5 2 ) - 5 1 f(a 1 , 5 ± ) - 5 2 f(a 2 , 5 2 )] 
I now substitute for / 

/(^ 8) = — exp(-pa - iq5/a a ). 
s 



(E.l) 



(E.2) 



I now take (3.141) as an ansatz to be substituted into (E.l) in order to deter- 
mine the unknown functions $i and Xa- This leads to the equation 



oo oo 

-p J da <j a+1 <S> 1 {a)e- p(T J d5 Xa (5)e^ q5 = 

-oo 
oo oo 

da x da 2 J d5i d5 2 K(a u a 2 )^ 1 (ai)^(a 2 )xa(S 1 )xa(S 2 ) x 



x 



(<ti + a 2 ) a e- p{ai+r72) exp - 



iq- 



{a l + a 2 ) a / 



-a\ e 



(E.3) 



For this equation to hold, we need that the r.h.s. show a factorization in the 
dependence on q, since the l.h.s. does. For this to be the case, we need 

°r ( rj a + rr a \ °r 

J d5 1 d5 2Xa (5 1 )xa(5 2 )ex V i-tq 1 ^ 2 2 J = J d5xa(5)e'^ s . (E.4) 



As is readily seen, however, this equation is satisfied by the stable symmetric 
Levy distributions and by these only. Indeed, the Fourier transform Xa(o) °f 
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Xa{$) must satisfy the functional equation 

Xa{q) = Xc 



:Q Xc 



(01 + 2 



(E.5) 



for which the only solution is exp(— Igl 1 /"). 



F The Integral Equation for the Gelling Case 



In this Appendix we derive the basic equation to be satisfied by the scaling 
function <£>(x) when convergence to scaling occurs in the sense of the second 
moment only, as is, for example, characteristic of the situation in which gela- 
tion takes place. We follow exactly the same steps as in Appendix A and refer 
the reader to this Appendix for more detailed explanations. 

The definition we shall use of convergence to scaling is the one given in (3.64). 
In order to ontain the equation for $(x) we shall again start from an obvious 
consequence, namely 



lim 

tl — >oo 



*2 



o, 



(F.l) 



where we have defined 



'« = it 



00 

M 2 (ty 1 J dmc(m,t)f 



m 



;f. 2 ) 



Let us now evaluate I(t): 

I(t) = W 2 A(t)s(t) x - 2 M 2 (t) - WB(t)^Q - WC(t) S ^ 



M 2 (t) 



s(ty 



(F.3) 



Here the A(t), B(t) and C(t) are defined by 



A(t) = n f dmidm 2 m 2 rnl K 



At)'*®. 



mi \ , (mi 

At) 1 ' 



s(t) \s(t) 



At), 



mi+m 2 \ (mi+m 2 



s(t) 

s(t) 4 
J m\m\ 



s(t) 
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B(t)=Jm*c(m,t)f(^ (F.4) 

oo / \ 

C(t)= J m 2 c(m,t)f 



m \ m 



s{t) s (ty 



Because of the scaling hypothesis, the quantities A(t), B(t) and C(t) all ap- 
proach limits which can be expressed as follows in terms of the scaling function 
$(x): 



A = J dxdyxK(x,y)<&(x)<&(y) [(x + y)f(x + y) — xf(x)] 
o 

oo 

B = J dxx 2 $(x)f(x) (F.5) 



oo 

C = J dxx 3 ®(x)f(x). 



In order to go further, we must distinguish the gelling case and the non- 
gelling case with A = 1. Let us first consider the former. It is then reasonable 
to postulate that the various quantities involved such as M 2 (t) diverge as a 
power-law as the critical time t c is reached. One thus finds that 

m) _ ci m_ c 2 

M 2 (t) t c -t s(t) t c -t { ' ' 



Assuming that the remaining term of (F.3) (proportional to A) also has the 
same order of magnitude yields 

s(t)^M 2 (t) = (F.7) 



We now need to evaluate the different constants. An appropriate choice of 
W allows to eliminate everything, save the ratio Ci/C 2 , which enters via the 
relation 

M 2 (t) = const. ■ s(t) Cl/C2 . (F.8) 



We therefore need to evaluate this exponent. To this end we make the following 
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(not really rigorous) observation: from the scaling hypothesis follows 



lid ( \ 

J mc(m,t)dm = — ^ J x<S>(x). (F.9) 



M (t) 



The r.h.s of (F.9) varies as e 2 T M 2 (t)/ s(t). If, however, we allow e to go as 
l/s(t), then the l.h.s should be constant. From this we deduce 

^ = 3 - r. (F.10) 



For large times, (F.l) reduces to 

t i im oo/7r(^-(3-r)5-C)) = o (F.ll) 

*1 



from which follows 

A-(3-r) J B-C = (F.12) 



which is exactly (3.67) as stated in the text. Further, from (F.7), (F.8) and 
(F.10) follows straightforwardly 

a = l + X-T (F.13) 



Now let us go over to the case in which we have A = 1. In this case we know 
that r is in a sense 2, but this is not appropriate for convergence in first 
moment, since in that case we must always have 



oc 

/ 



x$(x)dx < oo. (F.14) 



So let us consider instead convergence in second moment. Up to (F.5) every- 
thing remains strictly as above. Now we must attempt to compare the different 
contributions. Let us first find a connection between M 2 (t) and s(t). Following 
the same line of reasoning as above, we obtain from (F.9), which is also valid 
in this case, that 

M 2 (t) = const. ■ (F.15) 
v ; lns(t) v ; 
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From this we may now deduce, by identifying the orders of magnitude of the 
coefficients multiplying A and C, that 



s(t) 

s(t) = const. ■ — -Vr, (F.16) 
w Inst v 1 



from which finally follows 

s(t) = const. ■ exp {const. ■ Vtj , (F.17) 



which is the result obtained by van Dongen [108] in a slightly different fashion. 
We may now choose W in such a way as to eliminate the various constants 
involved, so as finally to obtain 

A — B — C — 0. (F.18) 



For this function we only require that x 2 Q(x) be integrable at the origin, so 
that a r = 2 power-law poses no problems. If we wish to cast the definition 
of convergence to scaling in a form that is as similar as possible to the usual 
one of the non-gelling case, we may reduce it to 



lim 

t— >oo 



oo oo 

\Tt J mc(m,t)f (me"^] dm = j x$(x)f(x)dx. (F.19) 



This again is identical to the ansatz discovered by van Dongen in [108]. Note, 
however, that the above reasoning could require modifications in the presence 
of additional logarithmic corrections to the singular behaviour of near 
the origin. Whether this may haoppen or not is not clear to me at present. 



G The Sum Kernel: The Technique of Characteristics 



In this appendix we present the general solution of (4.34) with initial condi- 
tions given by (4.35) by the method of characteristics: if we define ((t) and 
F(t) through the equations 



F = -F (G.l) 
t=S e- t -F, (G.2) 
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then it immediately follows that, if F((,t) is any solution of (4.34), then 
F C(0'^ satisfies (G.l). It therefore suffices to match the initial conditions. 
This is done as follows; define 



_C(0) = Co (G.3) 
F(0) = f((a) (G.4) 

where /(C) is defined via (4.35). The characteristic equations (G.l), (G.2) can 
now straightforwardly be solved together with the initial conditions (G.3), 
(G.4) to yield the result given in (4.36), (4.37). Nevertheless, we present the 
details of a particularly straightforward way to perform this algebra, since this 
will be used extensively in the sequel for more complex cases. 

We first display an explicit solution of (G.l, G.2) in the scaling limit. This 
limit will turn out to be that in which t goes to infinity while the following 
variables 

-s = Coe* -s = (e 2t (G.5) 
remain constant. The solution of (G.l, G.2) is given by 

F = f(( )e- t = -s e- 2t (G.6) 
C = [/(Co - /(0)] e~* + Co - /(Co) + /(0) = -e' 2t [s + \f (0)**]. (G.7) 

where the first equality represents the exact solution and the second involves 
the scaling limit defined in (G.5). Now let us note the following fundamental 
relation: in the scaling limit one has, because of the defining relation for <&(#) 
(3.14): 

oo oo 

J dmmc(m,t)e m<: -> J x^{x)e~ ix dx. (G.8) 



But the following relation holds for the l.h.s. of (G.8): 

J dmmc{m,t)e mC = =^ = (G.9) 

n ^ Co ^ s o 



where the partial derivatives in the final expressions refer to the scaling ex- 
pressions derived in (G.5). This now readily yields 

£|rfWIW = tttW (ai0) 
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From (G.5,G.7) one now obtains the necessary relation between s an d s, 
namely 



s = s o + if"(0)sl (G.ll) 



from which the final result eventually follows, namely 

C[x$(x)](s) = —. 1 (G.12) 
f + 2/"(0) S 



Finally, the results stated in the text can be obtained in an elementary way 
from this equation and (G.ll). 



H The Parity Dependent Kernel: Scaling Theory 



In the following, we develop the technical details of the parity-dependent ker- 
nel. The kinetic equations for the concentrations read 



C 2j+1 = M c 2fc+lC 2 (j-fc) - Kc 2j+ 1 C 2fe+1 - Mc 2j+1 C 2fe 
fc=0 fc=0 k=l 

K j ~ 1 L j ~ 1 

% = 7T Yl C 2fe+l C 2(i-fc)-l + — C 2kC 2 (j-k) ~ (H.l) 
Z k=0 Z k=l 

oo oo 

-Mc 2j c 2fe+i - Lc 2 j ^ °2k- 

k=0 k=l 

If one introduces the generating functions 



£(c,*)=E%+iWc 2j+1 

3=0 

oo 

^(C,t) = E%WC 2j (H.2) 

3=1 

and the corresponding moments 



d(t)=D(l,t) 

p(t)=P(l,t), (H.3) 
one eventually finds the following equations: 
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D = MDP - [(Mp(t) + Kd{t)]D 

P = ^-D 2 + ^P 2 - P[Lp(t) + Md{t)\ . (H.4) 

Note that these are ODE, the only dependence in ( coming in through the 
initial conditions. Equations (H.4) then yield the closed relations for d(t) and 
pit): 



d = -Kd 2 

p=\d 2 -\p 2 -Mpd. (H.5) 

These can be solved exactly. For our purposes, we shall only need the following 
elementary facts: 



/ x / ,/ x ~ K-M 
hm p{t)/d{t) =Poo = j 



KL 

1 + 4/1 + 



(k - My 



(H.7) 



Here a is the initial value of d(t), which plays no significant role in the follow- 
ing, since we shall always be interested in the large-time limit. 

Let us now introduce 



D((,9) = D((,t)/d(t) 

P((,9) = P((,t)/d(t) (H.8) 
p = p(t)/d(t) 

6(t) = Kd(t) 
One then finds: 



K^ = MD(P-p) 

K< % = j 32 + ¥ 2 + {K - M - Lp)p - (H - 9) 

The initial conditions are very close to the solution given by D — 1 and 
P{9) = p{9). They therefore follow it for considerable time and eventually can 
be described as solutions of the following equations 



dD ~ ,~ 
K— = MD(P-p OQ ) 
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^ = jD 2 + ^P 2 + (K-M-L Poo )P, 



(H.10) 



which is a two-dimensional autonomous system. The functions D and P then 
start near the equilibrium D = 1 and P = p^. It is then fairly straightforward 
to convince oneself that letting ( tend to one and simultaneously letting t tend 
to infinity in such a way that t(l — () remains constant will give a non-trivial 
limiting behaviour for both D and P. This then is easily seen to be equivalent 
to the existence of a scaling limit. The evaluation of the scaling function by 
these means, however, seems to be a hopeless task 

What is possible, however, is the evaluation of the r and $ exponents for the 
odd and even concentrations. To this end, we assume that 



lim f)(2j + l)c 2j+1 (t)f (?^±±) = Jx* d {x)f(x), 

oo ■ 00 

lim ^2jc 2j (t)f (J^ = J x%(x)f(x). 



(H.11) 



We may then proceed once more exactly in Appendix A in order to obtain an 
equation of the type of (3.29) or (3.34). The final result is 



r Mr 

J dxx 2 <S> d (x)f'(x) = —J dxdyx<S> d (x)%(y)f(x + y) 



r r m 

J dxdyx$ d (x) $ d (y) + — %(y) 



00 00 ^ 

J dxx 2 %(x)f (x) = \ J dxdyx $ d (x)$ d (y) + — %(x)%(y) 



Jdxdyx%{x) —<5> d (y) + —$ p (y) 



(H.12) 

f(x + y)- 
(H.13) 



for arbitrary continuous functions f(x). In order to fix the normalizations 
properly, the simplest way is to note that, from (H.ll) using f(x) — 1/x one 
finds 



00 

/ $ d (x) = lim [Ktdit)} = 1 

J t^oo 


00 

/ $ p (x) = lim [Ktp(t)\ = po 



(H.14) 
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Now, putting f(x) equal to e~ px into (H.12,H.13) and using the standard 
approach, we find from (H.12), that all terms except the first on the r.h.s. are of 
the same order for large p, namely of order p 2 ~ Td . The first is manifestly always 
subdominant. In order to evaluate one then works out all the prefactors and 
finds in the end 

M 

T d =l-— Poo (H.15) 



The equation (H.13) is a bit more complex, since one must distinguish between 
the case in which the first term dominates or is subdominant. All other terms 
are of order p 2_Tp , whereas the first is of order p 3_2rd . Assuming first that 
the first term is again subdominant leads once more, using a comparison of 
coefficients to 

TV, = 2 - ^ + -IP*,. (H.16) 



On the other hand, if the first term dominates, a simple power-law comparison 
leads to 

r p = 2r d -l. (H.17) 



Finally, in order to evaluate the exponent d for large cluster sizes we proceed 
again in the standard manner: we develop (H.12, H.13) around a singularity 
at some p c < 0. This time, it is straightforward to check that the d p exponent 
is zero in all cases. As for the ^ exponent, one can not easily say anything, 
except that it must be positive. If this is so, however, then the $ exponent of 
the full distribution, given by (4.46) is equal to zero. 

As a final remark concerning the scaling theory, note that equations similar 
to (H.12, H.13) can easily be derived in the much more general case 



K(k, I) 



a 1 K Q {k,l) 
a 2 K (k,l) 
a 3 K Q (k,l) 



for k and I odd 
for k and / even , 
otherwise 



(H.18) 



where K (k,l) is a homogeneous kernel that is otherwise arbitrary. From a 
study of these equations, it might perhaps be shown in general under what 
circumstances such changes in the scaling function are to be expected. In 
particular, it might confirm a suspicion that this is characteristic of Type II 
kernels, the scaling theory of which is in any case non-universal. 
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Finally, let us evaluate the non-scaling asymptotic behaviour. The easiest ap- 
proach is to evaluate first the long-time behaviour of Ci(t) and C2(t) and then 
to convince oneself through induction that this is indeed the large-time be- 
haviour at fixed j for j arbitrary. From (H.l) one finds for monomers and 
dimers 



ci = - [Kd(t) + Mp(t)} ci 

C2 = \c\- [Md(t) +Lp(t)]c 2 (H. 19) 

From these and (H.7) the relations (4.49) stated in the text are readily verified. 



I The g-sum Solution: A Sketch 



Here we shortly sketch the derivation of the solution (4.62,4.63) for the g-sum 
kernel. From (4.58) the following equation for H((,9) follows, see (4.61): 

— (C, 0) = H(C, 6) [H(C, 6) - H({ - 1, 0)] (1.1) 



This equation is considerably simplified by rewriting H((, 0) in the following 
quotient form: 



One then finds for h((, 9) the linear equation 

§(C,0) = -MC-M) (1.3) 



This equation is linear and translation invariant with respect to 0, so that it 
can be solved by Fourier transform. After some work one verifies that 



oo am oo 

ku) = e« e - n (i - eb(c_0 )" ( L4 ) 



m=0 m - l=m 



From this it is straightforward, though somewhat tedious, to obtain the final 
result stated in the text. 
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J Crossover in the g-sum kernel 



In the following I wish to show that je~ 2t and bj are the correct scaling 
variables for the crossover transition for the modified g-sum kernel (2 — q k — 
q l )/b. I do not claim any expression, even formal, for the scaling function: the 
methods I shall employ will not be precise enough to fix any kind of functional 
dependence. They will, however, determine the order of magnitude of various 
quantities, and this will be sufficient. 

If the above is correct, then we expect from the definition (3.90) for the scaling 
in crossover, that 

1 °r 

tf(C,0;ye- 2t )j = J S(x, y) (e^ - l) dx, (J.l) 



lim < y 

t—>oo 



E<*(*) 



.fc=] 



where we have singled out the dependence of H((, 9; b) on the argument b. To 
show this, I shall prove the following: first I shall identify the crossover regime 
in time as given by finite values of the parameter 



In other words, we require b6 to be close to one and the difference of order 
6 1 / 2 . Then I show that in this case one has 



E<*(0 



= (Jib' 1 ' 2 ) 



(J.3) 



Ufc=l 



and finally I show that in the crossover regime as defined above one has for 
the generating function H: 

H((,6) = Vb<f ) ((,0- (J.4) 



Here and throughout this appendix 23 will denote an arbitrary function of 
the scaling variables ( and £, which may vary from equation to equation. 

Let us first start with the elementary identity: 

1 ^ q ra f(q r ) ^ f {m \0) f a , /j r\ 



e q (q a )^ (q;q) r ^ m\ 



23 but nowhere else in this paper! 
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valid for any entire function / for which one of the sides converges. This is 
obtained by developing / in a Taylor series, inverting the two sums and using 
the two definitions of the g-exponential given in (4.60). 

Using (J. 5) one obtains from 4.64) the following expression for t 



Note the extra factor of b due to the fact that we are working with a modified 
kernel with an additional factor b~ l in its definition. It is now easy to convince 
oneself that the quantity 



decays rapidly when m becomes significantly larger than Vb and is constant 
before. We may therefore replace, as far as orders of magnitude go, this pref- 
actor by an appropriate cutoff of the sum, leading to 

t = V ^ (J.8) 

where N(b) = b^ 1 ! 2 . Using identity (P.ll) one then finds 

t = lniV(6) + 0(O, (J.9) 

from which (J. 2) follows. An entirely similar evaluation yields (J. 3) from the 
expression 



E<*(*) 



1 1 00 r 9q r 

= ^rEf^- (J.10) 



Again, a similar transformation applied to S((, 6) yields 



00 



Once more, we are led to cut the sum off at N(b) leading to 
N(b) (a- a^) 

S((,9)=J2(b9r [ -?^-^. (J.12) 

m=l W' Q) m 
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But, for this range of values of m one has, as q — > 1, 



m~ Q . (J.13) 



If one inserts (J.13) into (J. 12) and reduces the sum to an integral, it follows 
that 

S(C,6) = N(b) 1 -^(C,0, (J-14) 



from which (J. 4) follows straightforwardly. All the results initially claimed are 
hence established, and the statements made in the text follow. 



K The Simple Product Kernel: Technicalities 



In the following, we work out the detailed properties of the simple product 
kernel, which has all the necessary ingredients to understand the phenomenon 
of gelation. It is worthwhile to do this, as the full bilinear kernel is very much 
more complex. 

We use the generating function G((,t) defined in (4.70), which satisfies the 
PDE (4.71). The characteristic equations of this last are given by 



dC 
dt 



-\G-M l {t)] = -{G-l)-h{t) 



(K.1) 



where the second equality of the first equation defines h(t). From this follows 
immediately that 



G 



Co - fo(Co) - 1]* - / h(t')dt',t 



g(Co). 



(K.2) 



Setting h(t) identically zero for the time being, one obtains by the inverse 
Laplace transform of (K.2) 

mc(m,t) = — § dCo[l-*ff'(Co)MCo)exp{-m[Co-^(Co)-l)t]}(K.3) 

—ioo 
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The case in which (7 (Co) is e"* , which corresponds to monodisperse initial con- 
ditions, is now reduced to an exercise, which the reader can do in order to 
verify the result (4.74) stated in the text. 

We now need to gain some understanding of the failure of the solution (K.3) 
and of the assumption that Mi(t) remains constant. Equations (K.l) specify 
the surface G((, t) by displaying a pencil of straight lines as the level curves of 
this surface. These straight lines intersect, however, meaning that the surface 
must be multivalued. It is important to understand this issue well: indeed, 
since we have derived the PDE (4.71) without giving much thought to ques- 
tionns of convergence, we must now take care that the solution we obtain has 
the necessary analyticity properties. We may therefore not take into consid- 
eration the second sheet of the solution, but only that sheet which connects 
smoothly to ( = —00, since there the convergence of the series is unproblem- 
atic. 

We therefore need to determine the locus of points on the (C, t) plane for which 
the function G((,t) becomes singular. This will define a curve, the position of 
which we must determine in such a way as to obtain a meaningful solution for 
the equations (2.2) with the kernel (4.69). The determination of the singular 
curve proceeds as follows: If one takes two infinitesimally close straight lines 
from the pencil defined by the first equation (K.l), their intersection belongs 
to the singular curve, as is intuitively clear and readily verified analytically 24 . 
Take Co and Ci two infinitely close initial values for the first of the equations 
(K.l). This yields for the time at which these intersect 

*c(Co) = W(Co) (K.4) 



and hence 

(t \ - W(Co) 

Cc(Co) = Co - 9 -^%r— - [ h(t>)dt>. (K.5) 
9 \So) J 

As long as t < £ c (0), therefore, there is no problem and the solution (K.3) is 
valid. This is due to the fact that the function G((,t) is analytic at C = 
for such times, so that the PDE (4.71) can be applied in order to show that 
G(0,t), which is the first moment, is indeed constant. This therefore singles 



Note that a similar pencil exists in the sum kernel and also has intersections at 
finite time. The reason this creates no problem is that, for the sum kernel, all the 
computations are carried out using a rescaled time 9, and all such intersections 
occur at values of 6 which do not correspond to any physical time. See Appendix G 
for the detailed definitions. 
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out t c (0), which is given by 

t c = t c (0) = l/( ? , (0) = M 2 (0)- 1 (K.6) 

as the limiting time beyond which (K.3) ceases to hold. It is known as the gel 
time, and its determination is in general not a trivial task; (K.6) solves it for 
arbitrary initial conditions in the case of the product kernel (4.69). 

Beyond t c , we require a condition that allows to determine the amount of mass 
carried to infinity by the run-away aggregation process. The following remark 
allows to do this: since G((,t) is a Laplace transform of a positive real func- 
tion, its singularity nearest to the origin must be on the real ( axis. We may 
therefore identify the singular curve (Cc(Co), ^c(Co)), as defined in (K.5,K.6), as 
giving the nearest singularity in ( of G((,t) for fixed t. Here Co is simply a 
parameter for this curve. But since c(m, t) can at most decay as a power-law, 
since otherwise Mi(t) would remain constant, and since on the other hand 
c(m, t) cannot grow exponentially, we are left with the necessary condition 
that the singular curve coincide with the ( = line on the ((,t) plane for 
t >t c . In other words 

Cc[Co(*)]=0, (K.7) 

where one defines ( (t) through the relation 

t = W(Co(0)- (K.8) 

h(t) is then given by 

Ht) = ^-~j t [«o(01 " 1)] • (K.9) 

In the particular case of monodisperse initial conditions, this yields 

h(t) = i-r\ (K.io) 

which corresponds to the results that follows from the solution (4.74) stated 
in the text. To derive the solution itself once h(t) is known, it is enough to 
proceed with the expression (K.2) in the same way as we obtained (K.3) from 
(K.2) under the assumption that h(t) was identically zero. The final result is 
given by 



^ ioo 

mc(m,t) = — § dCo[l-*</(Co)]<7(Co) 



x 
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x exp < — m 



z 

Co -(g(to)- l)t-Jh(t')dt' 



(K.ll) 



Again, it is straightforward to verify that the Stockmayer solution given in 
the text follows from (K.ll) in the case of monodisperse initial conditions. 

Let us now finish the subject by deriving the scaling theory for the product 
kernel. We shall here derive it in some detail for this simpler case and later, 
in Appendix M, show tersely that the general case brings nothing new. The 
reader who has well absorbed the following will fill in the necessary details. 

Since we are dealing with a gelling system, we need to use the definition (3.15) 
of scaling of the n-th moment, for n > 1. In the following, we will take n equal 
to two. In principle, one could look at two limits: the large time limit, and the 
behaviour near t c . In the following, we will only consider the behaviour near 
t c . In fact, it turns out that the large-time behaviour is considerably simpler. 

Using the monodisperse solution, the moment equations or otherwise, it is 
easy to see that the typical size diverges as (t c — t)~ 2 near t c . We therefore 
wish to establish this rigorously for all initial conditions, and to obtain the 
scaling function. To this end we first need to compute the second moment. 
Using the moment equations (3.157) we find for the times before t c 

m '« = t^& = ^ (K ' 12) 



The limiting statement predicted by the definition of scaling in second moment 
is therefore given by 

oo 

\im {(t c - t)- 1 [G(s(t c - t) 2 ,t) -l]} = J x<$>(x)(e sx - l)dx. (K.13) 



To evaluate the function G at this value of its arguments we must solve 

( -[g((o)-i]t = s(t c -t) 2 . (K.U) 



Remembering that t c is l/g'(0), developing everything to lowest order and 
introducing w for (t c — t)" 1 ^ we reduce (K.14) to 

g"(0)-2g'(0) 2 w + 2</(0)s = 0. (K.15) 
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This then leads to 



f x$(x)(e sx - l)dx = g'(0)w = - 1 - yl - 2(3s 
J pi". 



(K.16) 



where (5 is defined as 

a _ fl"(0) 
P ^'(0)3 



(K.17) 



and from this finally follows the expression (4.78) given in the text. 



L The Product Kernel: Power-law initial Conditions 



In the following we treat the effect of power-law initial conditons on the gelling 
transition. For simplicity, we limit ourselves to the case of the product ker- 
nel. We shall use all the notation of Appendix K. We shall start from the 
assumption that the large m behaviour of c(m, t) is given by 



c{m,t) = Am _1 - Q [l + o(l)], 



(L.1) 



which implies for the generating funcion g(Co) 

[1 + CoS'(O) + AT(1 - aX-Co)"- 1 ] [1 + o(l)] 

9 Co = { 

[1 - AF(1 - a)(-( o r-i] [1 + 0(1)] 



(2 < a < 3) 
(1< a < 2) 



(L.2) 



We shall treat the two cases separately: in the first (2 < a < 3) the third 
moment is infinite, so that the scaling approach developed in Appendix K 
fails, but the second moment is finite, so that the derivation of the gel time 
in (K.6) holds good. In the second case (1 < a < 2), on the other hand, 
the second moment diverges. From a generalization of (K.6) one is therefore 
tempted to think that instantaneous gelation takes place, an impression that 
shall indeed be confirmed in this Appendix. 



In the first case, we must solve the equation 

( -[g(( )-l]t = s(t c -t) 1 ^, 



(L.3) 
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where a is as yet unknown. Developing to leading order in ( and t c — t results 



in: 



M C V{1 - a)(-Co) a ^ - (tc - t)C ° + s(t c - = 0. (L.4) 



The three terms are of the same order if we introduce the rescaled variable w 
as Co(^ c — t) 1 /(«-2) and take 

a= -. (L.5) 

a — i 



The scaled version of (L.4) then reads 

At c T{\ - a){-w) a ~ l -- + s = 0. (L.6) 



From the definition of scaling, one obtains using the same reasoning as in the 
derivation of (K.16) 

oo 

J x$(x)(e- sx -l)dx = g'{0)w (L.7) 
o 



From this and (L.6) one could presumably obtain the exact form of <&(#). This 
task is left to the enthusiastic reader. On the other hand, it follows from (L.6) 
that — w goes as sV( Q-1 ) as s — > oo, so that the r exponent can be found 
easily, and is given by 

r = ^=i. (L.8) 
a — 1 



Note that these results are in agreement with the scaling law (3.25) connecting 
r and a in the gelling case, which in this case is given by 

3 - r = a. (L.9) 



In order to evaluate the correction to scaling exponent A, we remark that the 
subdominant behavior of — w as a function of s for s — > oo is linear in s, so 
that 

A = ~~r (L.10) 

ct — 1 
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From (L.10) follows that A and a are identical, as predicted by the scaling 
theory developed in the text. 

Now let us evaluate the behaviour of h(t) shortly beyond the gel time. To 
this end we use (K.9), where we first determine the behaviour of ( (t) in this 
region. Since we are near the gel time, (o(t) is small, so that we may use the 
asymptotic expressions (L.2) for g((o). Using these yields 

h(t) = const. ■ (t - t c ) 1/{a - 2) (L.ll) 



This implies that the flow of mass to infinity is slower than in the case where 
the initial conditions are short-range, which does seem rather unexpected. In 
particular, this is an instance in which the flow rate of sol to gel is zero at 
the gel time, therefore strongly suggesting that r > 5/2, as is indeed borne 
out by (L.8). On the other hand, it is of course readily verified that, past the 
gel time, the singularity of G((,t) at ( = is always of square root type, so 
that the large j behaviour of the cluster size distribution for fixed times larger 
than t c always has an asymptotic power law decay as j~ 5 / 2 . 

Let us now turn to the second case. There it follows from (L.2) that the 
generalization of (K.4) is 



From this it now follows along the same lines as those presented above or in 
Appendix K, that 

h(t) = const. ■ t 1/{2 - a) (L.13) 



for small values of t, and that t c is indeed zero. Again this implies that the 
flow rate of sol to gel is initially zero, in remarkable contrast to the case in 
which the initial conditions are short range. This underlines the arbitrariness 
of the assumption (often made in the literature, see in particular [58]) that 
the flow rate at t c is of order unity. 



M The Full Bilinear Kernel: A Sketch of the Scaling Theory 



In this Appendix, we give a very short argument showing that the full bilinear 
kernel scales in the same way as the product kernel. As the calculations become 
quite involved I shall be brief. 
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If we define as usual 

oo 

G((,t) = Jc(m,t) (eO-l), 



(M.l) 



one obtains the following PDE gor G((,t) 
G t = \G 2 - BG (G c - 1) - \CG\. 



(M.2) 



The characteristic equations for this PDE can be written as follows: 



| = -G- C (,-1) 



ApG + Bp(q - 1) 



dp 
~dl 

^ = AqG + Bq(q-l) 
dG 



(M.3) 



dt 



■■p-q[BG + C(q-l) 



The curves defined by the above system of ODE's has properties similar to 
those of the characteristic equations in the quasilinear case. Specifically one 
has 



G 



C(t),t =G(t) 



G t ((t),t 



Gt 



C(t),t 



: P (t) 



(MA) 



In order to prove that gelation occurs at finite time, it is sufficient to look at 
the moment equations (3.157). In order to show that the scaling theory is the 
same, it is enough to show that the surface described by q (which corresponds 
to the function G of Appendix K has the same geometrical structure as the 
function G described in Appendix K. 

In order to see this, one makes the following elementary remarks: for the usual 
initial conditions with ( negative but very small, both q and G are small and 
negative. It therefore follows that dq/dt is also small and negative. Since the 
gel time is finite, this means that q is still of order ( at the gel time. Of G 
we can at least say that it is no larger than of order ( . The behaviour of ((t) 
is therefore qualitatively similar to the one found for the product kernel in 
Appendix K, that is, the integral curves for ((t) will cross at some finite time, 
which becomes the gel time as Co 0. We will therefore observe a similar 
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backbending of the surface described by q and hence a similar square root 
singularity, which will therefore lead to the same scaling function. 



N The B ^ Crossover When C = 

In this appendix we study the crossover behaviour of the kernel 

K(m,m') = l + e(m + m'). (N.l) 

For this we use the generating function 

oo 

G(C, t; e) = J dm'c(m\ t) [e m ' c - l] , (N.2) 



where the e dependence as been explicitly put in, as we shall have occasion to 
single it out in the sequel. G((,t) then satisfies the following PDE 

G t - eG ( = ^ - eG (N.3) 



The characteristic equations are given by 



dt 

dG G . . 

Ht = ~2~ (N ' 4) 

These equations have, as always, the property that 

G(C(t),t;e) = G(t:e) (N.5) 

We wish to show that scaling holds with s(t) equal to t. If this is true, it 

follows, using the definition (3.90) of the scaling limit in crossover, that 

oo 

hm [tG(s/t; y/t)] = f dx $(x, y) (e sx - 1) (N.6) 



Note that this is a legitimate use of (3.90), since it involves using as a function 
f(x) the expression (e x — l)/x, which is indeed bounded at the origin, as 
required. 
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From (N.5) and (N.6) it follows that we need to evaluate tG(t) as t — > oo 
while et is kept equal to y and t((t) is equal to s. This means that the initial 
conditions Co an d Co are varied appropriately. On the other hand, we have, 
as always, 

G(C , : e) = G(0 : e) = J dmc(m, 0) (e<<>™ - l) = Co + O(Co) (N.7) 



We now solve (N.4) explicitly with the initial conditions given by (N.7), where 
we shall neglect the higher order terms: 



G(f e) - 2 <° e 
G( '' e) "C e--(Co_-2e) 

C(t; e) = Co + 2e In ~f ° ~ ^ (N.8) 

From these results one finds that in the scaling limit it is necessary to take 
Co of order 1/t (note here in passing the considerable difference to the similar 
computations when the ordinary scaling limit is taken in the sum kernel: in 
that case one finds a quadratic relationship between Co an d C(0- Again this 
shows that we are in two quite different sorts of asymptotic regimes). If we 
now define w as the value at which we keep tCo fixed, we obtain in the scaling 
limit 



t, x „ , / we 9 - (w - 2y) 

s= t((t) = w + 2y\n\ — 



2y 



2ywe~y 



t G(t; e) = — (N.9) 

v ; we~y -(w-2y) K J 

From (N.6) and (N.9) one finally obtains, by substracting from the r.h.s. of 
(N.6) its limit as s — > oo 

/ dx $(x, y)e sx = xr y — , ^ (N.10) 



This together with the first of the relations (N.9) yields the result stated in 
the text, by using the following inverse transform 

1 / 4y2 \ c+ ;°° ds e ~^> 

$(x,y) = J- [-^—) I dw^~ — (N.ll) 

v ' y; 2m\l-e-v) J dw2y-w(l-e-y) y ' 

c—ioo 
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This integral can be evaluated and yields the result (4.83) as stated in the 
text. 



O The C -> Crossover When B = 



Here we treat the crossover from constant to product kernel defined by the 
kernel (4.80). We start by defining the following generating function: 



G((,t) = J dmc(m,t) 
o 



e mt -l 



(0.1) 



For this function one finds, in the pre-gel regime, which is the only one we 
shall be concerned with here: 



G t = \\G 2 + e 2 {G c -l) 



(0.2) 



We may now write the characteristic equations for this system, which is now 
quasilinear, in contrast to the others treated thus far. One finds 



p = pG 
q = qG 

G~ = p-e 2 <p(l + <t ) ). 



(0.3) 
(0.4) 

(0.5) 

(0.6) 



Here p and q denote G t and Gq respectively. Solving these equations with the 
appropriate initial conditions yields the usual property for the characteristic 
curves namely that they generate a solution surface: 



G 



C(t),t =G(t) 



(0.7) 



We also introduce the important auxiliary variable <f> — q — 1, which satisfies 
= G(l + 0) (0.8) 



As a first step, we want to scale away the parameter e from the equations. 
This can be done via the following transformations 

y = et p=\ G=% f=- (0-9) 
e 2 e 1 e 
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After these substitutions, the rescaled quantities satisfy (0.3-0.6) with e equal 
to one. The initial conditions are now given by 



Go — Co 
Po = |Co 2 > 



(O.10) 



where \i2 denotes the second moment of the initial conditions. It readily follows 
from (0.3) and (0.4) that q and p are proportional. If one now substitutes 



d_ 

dt 



lnp = G 



(0.11) 



into (0.6) one obtains using the proportionality of p and q 



2 
so 



(0.12) 



This can be straightforwardly solved by noticing the analogy to a Newtonian 
equation of motion and using conservation of energy. This eventually leads to 



p = 



C 2 
so 



2 + Co 2 



(Co cos y + 2Co sin yj ' 



(0.13) 



from which the result stated in the text readily follows upon noticing that, 
from the definition of the crossover limit 



J dxx<&(x, y)e~ sx = q(—s,y), 



(0.14) 



where q((,y) is precisely the function defined above in the scaling described 
in (0.9). 



P The Solution of the Finite Constant Kernel 



The equations for the finite constant kernel may be written as follows: 

Cj(t) = \ E - Cj (t)S N (t) (1 < j < N) (P.l) 

k=i 
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N 

S W (t) = E^(i). (P.2) 
k=i 

The strategy, which in fact is the one that generalizes straightforwardly to 
all other K^\k,l) is the following: first one evaluates Cj(t) for an arbitrary 
function Sn^) using only (P.l). One then considers (P.2) as a self-consistency 
equation for Sjv(£). After having solved for it, one has an explicit expression 
for Suit) and hence for Cj(t). 

In the case of (PI) the first step is reasonably easy: one transforms to new 
dependent and independent variables as follows 



^■(0) = Cj (t) exp yj S N (t') dt' 

dQ = \ exp ^- J S N {t') dt^j dt. (P.3) 
Equations (P.l) then become 

^ = (P4) 

fc=0 

which, for monodisperse initial conditions, have the solution 

^(6) = V-\ (P.5) 

The self-consistency condition (P.2) can now be written as 
df) N 

Sn(0) = 2-Y:^-\ (P.6) 
One now notes that 

_ d . do (dey 1 d 2 e 



If we now substitute (P. 7) into (P.6) one obtains, after transforming to the 
variables p = d6/dt and 9 

fc=i 
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from which the formulae given in the text follow immediately. 
We now proceed to give the asymptotic expressions for the sum 

N nk 

x(o) = Y. T - lnN 
k=i K 



(P.9) 



in the scaling regime, that is, when £ defined as N(l — 9) is of order one and 
N — > oo. We first observe that 



9 k - 1 f y N - 1 



k=l 



k 



f y N -. 

J y-1 







1 + ^ 



(P. 10) 



r e y — 1 
= / ^ + 0(iV- 1 ), 



where in the last equality one uses the scaling limit. The identity stated in the 
text then follows straightforwardly. 



Q The Finite-sum Kernel 



The general equations for the finite sum kernel read as follows 

d 3 = E k °k c j-k - jCjM (t) - CjMi(t), 
k=l 



(Q-i) 



where, originally, M (t) and M 1 (t) are arbitrary functions. We shall then im- 
pose the self-consistency conditiions 



M (f) = Ec,-(f) 

Mi(*) = Eic i (o. 

J'=l 



(Q.2) 



The first step is accomplished as follows: introduce the new variables 
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<t>j{e) = Cj (t) exp { / dt' [jMo(t') + Mi(0 



= exp 



L 



y dt'M^t') 



(Q.3) 
(Q-4) 



One then finds, denoting the derivative with respect to 6* by a prime 

3-1 

k=l 

from which one immediately obtains the solution by induction 
OjiO) = «;,(>' ■ 

where the aj satisfy the relation 

3-i 



U ~ l ) a j = ka kaj-k- 



k=i 



If one defines 

oo 
3=1 

then one finds immediately that 
Fe~ F = e<, 



(Q.5) 



(Q.6) 



(Q-7) 



(Q.8) 



(Q.9) 



from which both the explicit expressions for the a,j cited in the text as well as 
the various other statements made there are easily derived. 

Let us now define 



X = exp|y" dt' [Mi(t')] 
F = exp|y dt' [M (t')\ 



(Q.10) 
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From this we derive the equations 



N 
3=1 

N 

Y' = Y,a j 9 J - 1 Y 1 - j . (Q.ll) 
3=1 

If we now introduce A — Y/9, we obtain the solution for A described in the text 
from the second equation (Q.ll). From the first equation (Q.ll) one obtains 



dlnX 

Ef =1 j%A-i 



dln9 = ^ N - — . (Q.12) 



Identifying this with the correponding expression of d\n9 in terms of A one 
eventually obtains 

X = - ., , (Q.13) 

1 - E,-=i ajA 3 



from which all results quoted in the text follow. 



R The Three-body Kernel mi + m 2 + m 3 



If we define G to be F((,t) — S(t), with the definitions given in the text, we 
find that the characteristic equations for the PDE satisfied by G are as follows 



C = -\G{G + 2S) 

t=-SG (R.l) 

s = -s 2 . 

Here S is defined as in (4.105). Let us introduce 9 through the equation 

d9 = S(t)dt. (R.2) 

From this the general solution follows straightforwardly: 

C = Co-| /(C ^ ( " 0) /(0)2 (l-e- e ) (R.3) 
G=[f(Co)-f(0)]e- 9 - (R.4) 
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In order to compute the exact solution for the monodisperse initial condition, 
one first sets /(Co) equal to e'' in (R.4) and then uses the following represen- 
tation for the concentrations Cj(t): 

c j (t) = ±fG(C,t)e-XdC, (R.5) 

c 



where the contour C runs from —a — ixi to —a + ni, where a is an arbitrary 
real number, which should be large enough for all sums to converge. We then 
use the basic properties of characteristic equations, that is, 

G[C(0,Co),0]=B(0,Co). ( R - 6 ) 



Here, both ((9, Co) and G(9, Co) are explicitly known functions of Co via (R.4). 
It is therefore possible to rewrite the integral (R.5) in the Co variables. After 
some fairly tedious but straightforward evaluations of contour integrals, the 
result (4.109) is obtained. This method is quite general and can be used to 
evaluate virtually all kernels for which the generating function is obtained 
via characteristics. It is, however, almost invariaby quite cumbersome, as this 
example shows. 

In order to prove scaling, we proceed entirely as in the sum kernel case (see 
appendix G for the details of the approach followed, as well as the notation.) 
The scaling limit is defined by 

s = -(e 2e s = -Coe 9 . (R7) 



In this limit (R.3) yields the following relation 

s = s + ^P-sl (R.8) 



By the usual techniques this yields the results stated in the text. 
If, on the other hand, we start with initial conditions of the type 

c(m, 0) = Am- a (2 < a < 3), (R.9) 



one finds for large t and hence large 9, that 

AT(1 - aX-Co)"- 1 + Coe- e = C (RIO) 



147 



Hence, in order for all the terms in (R.10) to scale similarly, one must introduce 
the scaling variables 

s = -Ce^-W*- 2 ) s = -Coe e/{a - 2) . (R.11) 

From this immediately follows that 
a — 1 a 

Z = T = . 

a — 2 a — 1 



(R.12) 



S The Kernel m\ + m2 + Perturbing the Constant Kernel 



We wish to show here the scaling behaviour of the mixed two- and three-body 
system described in (4.117) in the specific case given by the first equation 
(4.117). For convenience I limit myself to the case of monodisperse initial 
conditions, but the general case presents no difficulties of principle. In this 
case, if one defines 



oo 

s(*)=E<*(*) (s.i) 

G{C,t) = F(C,t)-S(t), 
one obtains the following set of equations 



C 2 

G t -\uG{G + 2S)G z = — -aSG (S.2) 
S=-(l + a )S 2 . (S.3) 

Applying the method of characteristics to (S.2) one gets 

t=-±aG(G + 2S) (S.4) 
— G 2 

G = —-aSG (S.5) 

If we now introduce the rescaled time 6 defined by 

d6 = S(t)dt (S.6) 
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and define the initial conditions for G and C to be G and Co respectively, one 
obtains for G 



G(7 0) = G ° 6e/2 e -(V2+a)0 

Co^ 2 

_-(l/2+a)0 



where in the second line, the smallness of Co has been used as well as the 
following relation 

G(Co, 0) = e& - 1 = Co + Ko + O(C 3 )- (S.8) 



We now want to obtain scaling from this solution of (S.2). This is not quite 
a trivial task, so we carry it out in some detail. As was stated in (3.14), it is 
sufficient to show that there exists a function s(t) such that 



We therefore wish to express the l.h.s. of (S.9) in terms of a function of a 
scaling quantity alone. To do this, we first relate it to the solutions of the 
characteristic equations (S.4, S.5). Since, because of the definition of the char- 
acteristic equations, we have 

G[((t),t]=G(t), (S.10) 



one verifies using L'Hospital's rule, that 



lim = GAC(t),t), 



(S.ll) 



at least whenever ((t) ^0 when Co — > for times in the scaling limit, which 
is a condition we shall verify shortly. 

We therefore need an expression for C(£). In terms of 9, the following expression 
follows from (S.4,S.7) after some routine algebra: 



C = Co 



-Cue 9 / 2 



a(-Co) 



2a 



X 



-2a 



-Co 



2x + 3 
{l + xf 



dx + e 



-a0 



(S.12) 



From now on, we restrict ourselves explicitly to the case a < 1/2. Then the 
integral in the r.h.s. of (S.12) converges at the lower end. From this further 
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follows that we may take Coe 0//2 as the scaling variable. It follows that 



£(0) e (V2+«)0 = Z(-Coe fl / 2 ), (S.13) 

where Z(x) is a function which is expressed in terms of quadratures as follows 

y- 2a (2y + 3) 



Z(x) = —x 



1 + ax 2a J dy V —- 



(S.14) 



Note that ((9) — > as Co — ► at fixed value of so = Co e ^ 2 - Let us now define 
s as £(#)e (1 / 2+a - )6 \ We then find that these two quantities are connected by 

s = Z(s ). (S.15) 



We may therefore use the relation (S.ll) to evaluate the r.h.s. of (S.9). This 
leads to 

£ [x*(x)] (s) = — = 7T~~ \ ' (S- 16 ) 

Z(so){l-s ) s(l-so) 



This gives the scaling function explicitly although a more detailed evaluation 
would require considerable algebra. Since this is not our purpose here, however, 
we merely conclude these remarks by evaluating the large-s asymptotics of 
which immediately yields the value of the exponent r. But, in the 
limit of large s, one also has sq — > oo, so that the asymptotics of the r.h.s. of 
(S.16) is of order s' 1 , which corresponds to the value of one for the exponent 
r. On the other hand, since the scaling variable is (t, as readily follows from 
(S.3,S.6), it follows that the typical size grows as t and hence the growth 
exponent is one, as stated in the text. 

We now turn to the case a > 1/2. Since the results are less important and can 
be obtained by a straightforward extension of the previous approach I shall 
be brief. In the limit Co —> the relation (S.12) reduces to 

Co 2 + Coe-^ = C (S.17) 



l-2a 



by extending the integral to infinity. If we therefore now define 

Co = s e- ae C = se- 2a \ (S.18) 
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we obtain the following relation 

3a 



1 -2a 



s 2 o + soe = s. (S.19) 



If we use this in connection with (S.7), for which we did not require the 
hypothesis a < 1/2, we obtain in the scaling limit 

C [x$(x)] (s) = ^fl, (S.20) 



where sq is computed from s via (S.19). From this follows the result that 
r = 3/2, as stated in the text. 



T The Scaling Function for the One-dimensional PCM 



Let us show how to prove that (5.19) holds in general. In this case, the first 
moment is equal to p , so that we should properly first divide by it to obtain the 
correctly normalized scaling function. Let us compute the generaing function 
G((,t) given by 

G((,t)= Po - 1 Y / kc(k,t)e k <. (T.l) 
k=i 

After some straightforward algebra one obtains 

G((,t) = -= / y 2 exp 2p v / Dt(e c - 1) e~ y dy. (T.2) 

If one now defines 

s(t) = 2p v / ^Dt o = (s(t) (T.3) 



and noting that, as follows from the definition of scaling through generating 
functions stated in (3.14), one has 



4 v 



-= J y 2 e- y2 e ay = J x<5>{x)e ax dx, (T.4) 



from which the result stated in the text follows immediately through inverse 
Laplace transform. 
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U Ballistic Aggregation: Sketch of Exact Solution 

The analysis starts by the following basic fact: once the initial condition is 
fully known, the cluster size distribution function is known for all times. In the 
model considered in [75,28,29,30], the particle positions are fixed (they occupy 
a lattice, which we take to have unit spacing), so that everything is reduced to 
averaging over the initial momenta. This is in turn considerably simplified by 
the assumption that the initial momenta are independent Gaussian random 
variables. 

One must therefore express the event that, say, particles 1 to m form one 
cluster at time t in terms of the initial momenta of all particles. After this, 
one evaluates the probability of this event using the independence of the initial 
momenta. In the following, we shall make heavy use of the fact that, if a cluster 
has formed at time t, the center of mass of the particles that constitute it 
has been moving on a straight line with the cluster's velocity ever since the 
beginning. 

There are three conditions (or rather, classes of conditions) which must all be 
satisfied for this event to occur: 

(1) A set of conditions guaranteeing that the particles 1 to m do in fact come 
to form a cluster. One sees that it is necessary and sufficient that, for all 
mi, the particle cluster from 1 to mi have sufficient velocity to collide 
with the corresponding cluster from mi + 1 to m before time t. This 
is easily seen to be equivalent to the following set of conditions on the 
velocities v (m 2 ) of the particles at site m 2 with respect to the center of 
mass of the cluster 



for all mi. Now we see a connection with the functions defined in the 
text: due to the independence of the v(m2), the partial sums on the l.h.s. 
of (U.l) perform a random walk starting from the origin and returning 
to it. Conditions (U.l) then state that this random walk must always 
remain below the parabolic profile given by the r.h.s. of (U.l). 
(2) A set of conditions guaranteeing that the particles to the right of the 
cluster will not themselves coalesce with the cluster. If we again measure 
all velocities with respect to that of the center of mass of the forming 
cluster, one obtains the following condition for all velocities of particles 
to the right of the cluster: 




m2=l 





(U.2) 



m,2=m+l 
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where v denotes the velocity of the forming cluster in the rest frame. 
These conditions state that no cluster formed by particles m + 1 to m\ 
can collide before time t with the center of mass of the forming cluster 
which, by assumption, is at rest in the frame we consider. Again, we see 
that a given random walk is constrained to stay below a certain parabolic 
profile. In this case, this must occur for all values of the upper limit of 
the sum, so we are dealing with events of the type described by T^. 
(3) An entirely similar set of conditions for particles to the left of the forming 
cluster. The equations are, of course, entirely similar to (U.2), with v 
replacd by — v. 

From these considerations the results stated in the text can be obtained more 
or less straightforwardly. The crucial remark is the following: in the scaling 
limit, that is, when the cluster masses m are large, and f > 1 as well, the 
random walks will tend, under an appropriate scaling of the variables, to a 
continuous diffusion. The fact that all random walks start on the absorbing 
parabola causes technical problems, which are excellently discussed in the 
various references given above. 

As a final excercise for the interested reader, consider the following problem: 
if the initial velocities are independent, but distributed according to some 
Levy stable law, what can we say about the solution? Everything up to (U.l) 
and (U.2) still holds. The real difficulty consists in building up a continuuum 
model, since there is, to my knowledge, no Levy analogue to Brownian motion. 
The discrete problem, on the other hand, is quite formidable. The numerical 
results found in [45] make it probable that a scaling picture does indeed ex- 
ist when the average velocity is still finite, whereas the other case is quite 
pathological. 
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